Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	28,116	T→A	6.2%	L275Q (CTG→CAG)	rihC →	ribonucleoside hydrolase 3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	28,116	0	T	A	6.2%	226.8 / 4.4	65	L275Q (CTG→CAG)	rihC	ribonucleoside hydrolase 3
Reads supporting (aligned to +/- strand):  ref base T (29/32);  new base A (2/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

GTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACA  >  NC_000913/27831‑28389                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                  gTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGACTTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTggg                                                                                                                                                                                                                                                                                >  1:35566/1‑289 (MQ=255)            gTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAgccagcca                                                                                                                                                                                                                                                                     <  2:416659/289‑1 (MQ=255)                       tttGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGGAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTg                                                                                                                                                                                                                                                          >  1:314588/1‑289 (MQ=255)                       tttGCATGTCACCAGTCAGGAAATATTAACTCCGGACTGTCTCTCTACCCTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGt                                                                                                                                                                                                                                                           <  1:112628/288‑1 (MQ=255)                          ggATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCa                                                                                                                                                                                                                                                       <  2:294247/289‑1 (MQ=255)                                cACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGAt                                                                                                                                                                                                                                                 <  2:35566/289‑1 (MQ=255)                                             atatTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAg                                                                                                                                                                                                                                    >  2:163430/1‑289 (MQ=255)                                                        tGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCAt                                                                                                                                                                                                                                                      <  2:220729/260‑1 (MQ=255)                                                        tGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCAt                                                                                                                                                                                                                                                      >  1:220729/1‑260 (MQ=255)                                                        tGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCcagcag                                                                                                                                                                                                                         >  2:390307/1‑289 (MQ=255)                                                           cTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAgtg                                                                                                                                                                                                                       <  2:274481/288‑1 (MQ=255)                                                               ctctctACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAgtgggtg                                                                                                                                                                                                                   >  2:15947/1‑288 (MQ=255)                                                                  tctACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTg                                                                                                                                                                                                               <  1:106852/289‑1 (MQ=255)                                                                       acTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTg                                                                                                                                                                                                          <  1:15947/289‑1 (MQ=255)                                                                            cGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGaa                                                                                                                                                                                                                                      >  1:146623/1‑256 (MQ=255)                                                                            cGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGaa                                                                                                                                                                                                                                      <  2:146623/256‑1 (MQ=255)                                                                            cGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTgg                                                                                                                                                                                                      >  2:68448/1‑288 (MQ=255)                                                                            cGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTgg                                                                                                                                                                                                      >  2:376327/1‑288 (MQ=255)                                                                            cGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCAGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTgg                                                                                                                                                                                                      >  2:458278/1‑288 (MQ=255)                                                                            cGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCAGGGCAAGCCAGCCAATGGACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTg                                                                                                                                                                                                       >  2:458274/1‑287 (MQ=255)                                                                                            cGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCACTCTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAAc                                                                                                                                                                                       <  2:159198/289‑1 (MQ=255)                                                                                                 aaaTGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTcac                                                                                                                                                                                >  2:163315/1‑289 (MQ=255)                                                                                                 aaaTGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTcac                                                                                                                                                                                <  2:390905/289‑1 (MQ=255)                                                                                                                   ttAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGt                                                                                                                                                              <  1:163315/289‑1 (MQ=255)                                                                                                                   ttAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGt                                                                                                                                                              <  1:163430/289‑1 (MQ=255)                                                                                                                      gCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGCGGTCGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGtcat                                                                                                                                                           >  2:115781/1‑289 (MQ=255)                                                                                                                          cTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGtcattca                                                                                                                                                        >  1:423131/1‑288 (MQ=255)                                                                                                                            aCCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCg                                                                                                                                                     >  1:143776/1‑289 (MQ=255)                                                                                                                                    ggCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAgccagcca                                                                                                                                                                                                                                                                     >  1:273285/1‑169 (MQ=255)                                                                                                                                    ggCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAgccagcca                                                                                                                                                                                                                                                                     <  2:273285/169‑1 (MQ=255)                                                                                                                                          aTGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGTTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCttt                                                                                                                                         <  2:303053/287‑1 (MQ=255)                                                                                                                                          aTGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCttt                                                                                                                                         <  2:143776/287‑1 (MQ=255)                                                                                                                                               aggCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTAAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGTGTGGCTGCGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCACTCATCGACTCATGCCTTTCACTGat                                                                                                                                  <  1:261790/287‑1 (MQ=255)                                                                                                                                                   ggCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATAtc                                                                                                                               >  1:469527/1‑288 (MQ=255)                                                                                                                                                   ggCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATAtc                                                                                                                               <  2:423131/288‑1 (MQ=255)                                                                                                                                                     cTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAag                                                                                                                                                                                                                                                                            >  2:378009/1‑145 (MQ=255)                                                                                                                                                     cTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAag                                                                                                                                                                                                                                                                            <  1:378009/145‑1 (MQ=255)                                                                                                                                                      ttGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATAtccctc                                                                                                                           <  2:314588/289‑1 (MQ=255)                                                                                                                                                                               tCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTa                                                                                                  >  1:302639/1‑289 (MQ=255)                                                                                                                                                                               tCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTa                                                                                                  >  2:385397/1‑289 (MQ=255)                                                                                                                                                                                   ctggctggTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCa                                                                                               <  1:390307/288‑1 (MQ=255)                                                                                                                                                                                            gcgcCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCa                                                                                                                                                                                                                                                       >  2:363689/1‑127 (MQ=255)                                                                                                                                                                                            gcgcCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCa                                                                                                                                                                                                                                                       <  1:363689/127‑1 (MQ=255)                                                                                                                                                                                             cgcCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAAc                                                                                                                                                                                       <  2:73624/190‑1 (MQ=255)                                                                                                                                                                                             cgcCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAAc                                                                                                                                                                                       >  1:73624/1‑190 (MQ=255)                                                                                                                                                                                                          tCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGTATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCg                                                                       <  2:505583/289‑1 (MQ=255)                                                                                                                                                                                                          tCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCg                                                                       <  2:302639/289‑1 (MQ=255)                                                                                                                                                                                                          tCA‑‑CTCTCAAACCCTGTTTGATGGCAGTGTAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATGGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCg                                                                       <  1:68448/289‑1 (MQ=255)                                                                                                                                                                                                                tttcAAACCCGGTTGTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAAGCACGACGGTGGGTGATATCGACGGTGGCCTGGGCAAGCCAGCCAAAGTACAGGTGGCATTGGATCTGGGTGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTGATTTCTAATTATCAGCGTTTTTGGCTAGCGGCGTAGCGATgc                                                                   <  1:115781/287‑1 (MQ=255)                                                                                                                                                                                                                tctcAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATg                                                                    >  1:435270/1‑288 (MQ=255)                                                                                                                                                                                                                     aaCCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTa                                                                                                             >  2:120715/1‑242 (MQ=255)                                                                                                                                                                                                                     aaCCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTa                                                                                                             <  1:120715/242‑1 (MQ=255)                                                                                                                                                                                                                                ggggCAGTGGAAACGCAGGGCGAATTTACCTCAGGGACGACGGTGGTTGATATCGACGGTTGCCAGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGGGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaaa                                                   <  1:458278/287‑1 (MQ=255)                                                                                                                                                                                                                                 tgGAAGTGGGAACTCAGGGCGAATTTACCTCAGGCACGCCGTTGGTTGATATCGACGGTTGCCAGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGaaa                                                   <  1:458274/288‑1 (MQ=255)                                                                                                                                                                                                                                                 gggCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTAc                                                                                                                                                                                                                                                               >  2:317673/1‑68 (MQ=255)                                                                                                                                                                                                                                                 gggCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTAc                                                                                                                                                                                                                                                               <  1:317673/68‑1 (MQ=255)                                                                                                                                                                                                                                                 gggCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTa                                  >  1:233893/1‑289 (MQ=255)                                                                                                                                                                                                                                                            ccTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCGTTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGaa                       <  2:412487/289‑1 (MQ=255)                                                                                                                                                                                                                                                            ccTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGaa                       >  1:376472/1‑289 (MQ=255)                                                                                                                                                                                                                                                                aGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATa                     >  2:190703/1‑287 (MQ=255)                                                                                                                                                                                                                                                                          gtggtTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATa                     >  1:286016/1‑277 (MQ=255)                                                                                                                                                                                                                                                                          gtggtTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATa                     <  2:286016/277‑1 (MQ=255)                                                                                                                                                                                                                                                                                  tatCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACa  >  1:503089/1‑288 (MQ=255)                                                                                                                                                                                                                                                                                             gCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGaa                       <  1:4523/256‑1 (MQ=255)                                                                                                                                                                                                                                                                                             gCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGaa                       >  2:4523/1‑256 (MQ=255)                                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                                  GTATTGAAATCGTCATGTGCGGTTTGGATGTCACCAATCAGGCAATATTAACTCCTGACTATCTCTCTACACTGCCGCAGTTAAACCGTACCGGGAAAATGCTTCACGCCCTGTTTAGCCACTACCGTAGCGGCAGTATGCAAAGCGGCTTGCGAATGCACGATCTCTGCGCCATCGCCTGGCTGGTGCGCCCGGACCTGTTCA‑‑CTCTCAAACCCTGTTTTGTGGCAGTGGAAACTCAGGGCGAATTTACCTCAGGCACGACGGTGGTTGATATCGACGGTTGCCTGGGCAAGCCAGCCAATGTACAGGTGGCATTGGATCTGGATGTGAAAGGCTTCCAGCAGTGGGTGGCTGAGGTGCTGGCTCTGGCGTCGTAACCTGTCACATGTTATTGGCATGCAGTCATTCATCGACTCATGCCTTTCACTGATATCCCTCCCTGTTTATCATTAATTTCTAATTATCAGCGTTTTTGGCTGGCGGCGTAGCGATGCGCTGGTTACTCTGAAAACGGTCTATGCAAATTAACAAAAGAGAATAGCTATGCATGATGCAAACA  >  NC_000913/27831‑28389

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑