Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	802,526	C→T	L214L (CTG→TTG)	ybhI →	putative tricarboxylate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	802,526	0	C	T	96.6%	175.7 / NA	59	L214L (CTG→TTG)	ybhI	putative tricarboxylate transporter
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (19/38);  total (21/38)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGCTCCGGCAACACCGTCTAACACCGCGCGTGCGGGCGGCATTGTGTTACCGATCATCAACAGCGTGGCGGTGGCTTTGGGGTCCGAACCGGAAAAAAGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACCT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCTGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGG  >  NC_000913/802299‑802762                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              gggCTCCGGCAACACCGTCTAACACCGCGAGTGCGGGCGGCATTGTGTTACCGATCATCAACAGCGTGGCGGTGGCTGTGGGGTCCGAACCGGAAAAAAGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCGCTT‑GCAAATTAGCTGGGGTGGATgg                                                                                                                                                                                                                       <  2:550804/250‑1 (MQ=255)  ggCTCCGGCAACACCGTCTAACACCGCGCGTGCGAGCGGCATTGTGTTACCGATCATCAACAGCGTGGCGGTGGCTTTGGGGTCCGAACCGGAAAAAAGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTGCCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATgg                                                                                                                                                                                                                       <  1:482865/250‑1 (MQ=255)      ccGGCAACACCGTCTAACACCGCGCGTGCGGGCGGCATTGTGTTACCGATCATCAACAGCGTGGCGGTGGCTTTGGGGTCCGAACCGGAAAAAAGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGTGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGgcg                                                                                                                                                                                                                    <  2:144019/249‑1 (MQ=255)      ccGGCAACACCGTCTAACACCGCGCGTGCGGGCGGCATTGTGTTACCGATCATCAACAGCGTGGCGGTGGCTTTGGGGTCCGAACCGGAAAAAAGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGgcg                                                                                                                                                                                                                    <  2:184358/249‑1 (MQ=255)            acacCGTCTACCACCGCGCGTGCGGGCGGCATTGCGTTACCGATCATCAACAGCGTGGCGGGGGCTTTGGGGTCCGAACCGGAAAAAAGCCCGCGTCGTATCGGACATTACCTGATGATGGCCATTTACATGGTCACCAAAACAACCAGCTATATGTTCTTTACCGCAATGGCGGGGGACATTTTTGCGCCCAGAATGAACGAAGA‑CAATCTGCCACT‑GCAAATTAGCTGTGGTGGATGGGCGCCGGCAGc                                                                                                                                                                                                            >  2:150853/1‑251 (MQ=255)                     aaCACCGCGCGTGCGGGCGGCATTGTGTTACCGATCATCAACAGCGTGGCGGTGGCTTTGGGGTCCGAACCGGAAAAAAGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGAtt                                                                                                                                                                                                      <  1:527168/248‑1 (MQ=255)                       cacCGCGCGTGCGGGCGGCATTGTGTTACCGATCATCAACAGCGTGGCGGTGGCTTTGGGGTCCGAACCGGAAAAAAGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCgg                                                                                                                                                                                                 <  2:36099/251‑1 (MQ=255)                        cccGCGCGTGCGGGCGGCATTGTGTTACCGATCAGCAACAGCGTGGCGGTGGCTTTGGGGTCCGAACCGGGAAAAAGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCACAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCggg                                                                                                                                                                                                <  2:521158/250‑1 (MQ=255)                             gcgTGCGGGCGGCATTGTGTTACCGATCATCAACAGCGTGGCGGTGGCTTTGGGGTCCGAACCGGAAAAAAGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGcatca                                                                                                                                                                                           <  2:622252/251‑1 (MQ=255)                                  cgggcggCATTGTGTTACCGATCATCAACAGCGTGGCGGTGGCTTTGGGGTCCGAACCCGAAAAAAATCCGCGTCGAGTCGGACATTACCTGATGATGGCCATTTACATGGTCACCAAAACCACCAGCTAAATGTTCTTTACAGCAAAGGCGGGGAACAATTTGGCGCTGGAAATGATCAAAGA‑CATTTTTGACTT‑TCAAAATAGGTGGGGGGGATTGGGGCTGGCAGGCGGAATGCCCGGcatccctatg                                                                                                                                                                                      >  1:400854/1‑245 (MQ=255)                                    ggcggcATTGTGTTACCGATCATCAACAGCGTGGCGGTGGCTTTGGGGTCCGAAACGGAAAAAAGTCCGCGTCGTGGCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGGAAAATCTTGCGCTCAAAATGATTCACGAACAT‑CTG‑ACCTTGGAAATTATCTGGGGGGGATGGGCGCTGGGCCCCGGATTGCCGGGGACCCTTAtgct                                                                                                                                                                                    >  2:262383/1‑251 (MQ=255)                                      cggcATTGTGTTACCGATCATCAACAGCGTGGCGGTGGCTTTGGGGTCCGAACCGGAAAAAAGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTAtgctgc                                                                                                                                                                                  <  1:631012/251‑1 (MQ=255)                                      cggcATTGTGTTACCGATCATCAACAGCGTGGCGGTGGCTTTGGGGTCCGAACCGGAAAAAAGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAAAGA‑CATTCTGCACTT‑TCAAATTTGCTGGGGGGGATGGGGGCTGGCAGCCGGATTGCCCGGCCTCATTAtgctgt                                                                                                                                                                                  >  1:298113/1‑250 (MQ=255)                                          aTTGTGTTACCGATCATCAACAGCGTGGCGGTGGCTTTGGGGTCCGAACCGGAAAAAAGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGt                                                                                                                                                                              <  2:238734/251‑1 (MQ=255)                                             gtntTACCGATCATCAACAGCGTGGCGGTGGCTTTGGGGTCCGAACCGGAAAAAAGGCCGCGTCGTGTCGTACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAAGTAGCTTGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATTCTGTTGGTCa                                                                                                                                                                            <  2:598435/250‑1 (MQ=255)                                               gtTACCGATCATCAACAGCGTGGCGGTGGCTTTGGGGTCCGAACCGGAAAAAAGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCAcc                                                                                                                                                                          <  2:583002/250‑1 (MQ=255)                                                   ccGATCATCAACAGCGTGGCGGTGGCTTTGGGGTCCGAACCGGAAAAAAGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATCATCAACTA‑CATTCTGCACTT‑GCAAATTAGCTGTGGTGGATGGGCGCTGGCAGCCTGATTGCCGGGCATCAGTATGCTGCTGGTCACCCCGCt                                                                                                                                                                     <  1:211186/251‑1 (MQ=255)                                                            aaCAGCGTGGCGGTGGCTTTGGGGTCCGAACCGGAAAAAAGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTTGGGCTGAAAATTATTAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGGGTGCCGGGGATAATTATTCTGCTGGGCACCCCGCTGGTGGttt                                                                                                                                                             >  1:351727/1‑250 (MQ=255)                                                            aaCAGCGTGGCGGTGGCTTTGGGGTCCGAACCGGAAAAAAGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTa                                                                                                                                                            <  1:158694/251‑1 (MQ=255)                                                            aaCAGCGTGGCGGTGGCTTTGGGGTCCGAACCGGAAAAAAGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTa                                                                                                                                                            <  2:582661/251‑1 (MQ=255)                                                                       ggTGGCTTTGGGGTCCGAACCGGAAAAACGTCCGCGACGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAATATTCTGGCCCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTCGCTGGGGTGGATGGGCGCGGGCAGCCGGGTTGCCGGGCAACCTTATGCTGCTGTTCACACTGCTGGTGATTTACACCATGTATc                                                                                                                                                 >  2:323082/1‑251 (MQ=255)                                                                        gTGGCTTTGGGGTCCGAACCGGAAAAAAGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCAGCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATcc                                                                                                                                                <  2:125831/251‑1 (MQ=255)                                                                                  ggTCCGAACCGGAAAAAAGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGGACATTCTTGCGCTTAAAATGATCAAAGA‑CATTCTGCACTT‑GGAAATTAGCTGGGGGGGATGGGCGCCGGGAACCGGATTGCCGGACATTATTATGCTTCTGGTCACCCCCCCGGTGATTTACCCCATGTCTCCACCCGAAAtt                                                                                                                                      >  2:491631/1‑251 (MQ=255)                                                                                         aCCGGAAAAAGGTCCGCGTCGTGTCGGACATTACCTGCTGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACCATCTTGCGCTGAAAATGCACAACAA‑CCTTCTGTCCTT‑GCAAATTAGCTGGGGGGGGTGGGCGAGCGCAGCCGGGTTGCCGGGCATCTTTATGCTGCTGGGCACCCCGCTGGGTATTTTCCCCACGTATCCCCCAAAAATTataacgg                                                                                                                               >  2:5039/1‑245 (MQ=255)                                                                                             gAAAAAAGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGGAAATTCTGGCGCTGAAAATGATCAACGA‑AATTCTTCACTT‑GAAAATTAACTGGGGGGGATGGGGGCCGGCAGCCCGATTGCAGGGCCTTATTTTTCTGCTGGGCACCCCCCTGGTGAATTTCACCATGTATCCACCAGAAATTAAGAAGGTgg                                                                                                                            >  2:518882/1‑250 (MQ=255)                                                                                              aaaaaaGTCCGCGTCGTGTCGGACATTTCCTGATGATTTCCCTTTACATGGTCACCAAAACCTCCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGT‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGAATTACACCATGTATCCACCAGAAATTAAGAAGGTGGAt                                                                                                                          <  2:332239/251‑1 (MQ=255)                                                                                                 aaaGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAAc                                                                                                                       <  2:306716/251‑1 (MQ=255)                                                                                                  aaGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGGGCTGAAAATGATTAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTTCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATACCCCAGAAATTAAGAAGGGGGATAAc                                                                                                                       >  2:365056/1‑250 (MQ=255)                                                                                                      ccGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACaaaa                                                                                                                   <  2:163773/250‑1 (MQ=255)                                                                                                          gtcgtGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGTATGGGCGCTGGCAGCCTGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAATCCATc                                                                                                              <  2:631003/251‑1 (MQ=255)                                                                                                          gtcgtGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATc                                                                                                              <  1:335421/251‑1 (MQ=255)                                                                                                                 cGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCgg                                                                                                                                                                                                         <  1:195048/153‑1 (MQ=255)                                                                                                                 cGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGGGGATGGGCGCTGGCAGCCgg                                                                                                                                                                                                         >  2:195048/1‑153 (MQ=255)                                                                                                                      aTTACCTGATGATTTCCATTTACATGGTCACCAAAACCACCAGATATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATCGCTAAAGCGGGc                                                                                                  <  1:339763/251‑1 (MQ=255)                                                                                                                         aCCTGATTATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATCGCTAAAGCGGGCCt                                                                                                <  2:414364/250‑1 (MQ=255)                                                                                                                           cTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACTATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATCg                                                                                                             <  1:570106/235‑1 (MQ=255)                                                                                                                           cTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGGCCCCCCGCTGGGGATTTACACTATGTATCCACCAGAAATTAAAAAGGTGGATAACAAAACCATCg                                                                                                             >  2:570106/1‑235 (MQ=255)                                                                                                                             gatgatGTCCATTTACATGGGCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTTAAAATGATCAACGA‑CATTGTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGGTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATCGCTAAAGCGGGCCTTGc                                                                                             <  2:186475/249‑1 (MQ=255)                                                                                                                                   gTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAg                                                                                     <  1:566584/251‑1 (MQ=255)                                                                                                                                     ccATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAgg                                                                                    <  2:471087/250‑1 (MQ=255)                                                                                                                                     ccATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGa                                                                                   <  2:369487/251‑1 (MQ=255)                                                                                                                                          tACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAAAGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGCTGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGGAATTAAGAAGGTGGATAACAAAAACATAGCTAAAGCGGGCCTTGCCGAACTAGGACCGAt                                                                              >  1:611996/1‑251 (MQ=255)                                                                                                                                                gTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGCACATTTTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGGTGGGGTGGGTGGGGGGTAGCAGCCGGATTTCCGGACCTCATTATGCGGCGGGTCCCCCCCCCAGGGACTTACACCATGTATCCACCAGCCATTAAGAAGGTGGCTAACAAAACCATCGCTCAAGGCGGCCTTTCCGacattggacccctgaaaac                                                                        >  2:407164/1‑233 (MQ=255)                                                                                                                                                   cccAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACTCCATGTATCCGCCAGAAATTAAGAAGGGGGAGAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCcg                                                                     <  1:87308/250‑1 (MQ=255)                                                                                                                                                               aGCTATATGTTCTTTACCGCAATGGCGGGGAACATTTCGGCGCTGAAAATGATCCACAA‑CATTCTACACTT‑TCAAATTAGCTGGGTGGGATTGGCGCTGGCAGCCCGATTGCCGGGCCTCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCTACCTGAAATTAAGAAGAGGGATAACAAAACCATCGCTAAAACGGCCCCTGCCCGAGTAGGACCGCTGAAAATCCGCGAAAAAAtgcg                                                         >  1:217611/1‑250 (MQ=255)                                                                                                                                                                  tataTGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGAGGATGGGCGCTGGCAGCCGGATTTCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGGTTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATCCCTAAAGCGGGGCTTGCCGAAATAGGACCGATGAAAAACCGCGAAAAAAtgcttcc                                                      >  1:256668/1‑248 (MQ=255)                                                                                                                                                                   atatGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATCGCTAAAGCGGGCCGTGCCGAACTAGGACCGAGGAAAATCGGCGAAAAAATgctgct                                                      <  2:400854/250‑1 (MQ=255)                                                                                                                                                                     atGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCTGCTCgg                                                   <  2:192100/251‑1 (MQ=255)                                                                                                                                                                              aCCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAAGGCTGCTCGGTGTCTTtg                                           <  1:530784/250‑1 (MQ=255)                                                                                                                                                                                   aaTGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAAt                                                                                                                                       <  1:69529/153‑1 (MQ=255)                                                                                                                                                                                       gCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCTGCTCGGTGTCTTTGTGCTGGCgct                                 <  1:300512/251‑1 (MQ=255)                                                                                                                                                                                             aaCATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCTGCTCGGTGTCTTTGTGCTGGCGCTgctggg                           <  1:377620/251‑1 (MQ=255)                                                                                                                                                                                               cATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATCGCt                                                                                                           >  1:505163/1‑169 (MQ=255)                                                                                                                                                                                               cATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATCGCt                                                                                                           <  2:505163/169‑1 (MQ=255)                                                                                                                                                                                                  tCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCTGCTCGGTGTCTTTGTGCTGGCGCTgctgggctgg                       >  2:125399/1‑250 (MQ=255)                                                                                                                                                                                                           gAAAATGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCTGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGt              >  1:383512/1‑250 (MQ=255)                                                                                                                                                                                                             aaaTGATCAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCTGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAg           <  2:162468/251‑1 (MQ=255)                                                                                                                                                                                                              aaTGATAAACGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGAATCATTATGATGCTGGTCACCCCGCTGGTGATTTACACAATGTATACACCAGAAATTAAGAAGGTGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCTGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGt          >  2:213926/1‑251 (MQ=255)                                                                                                                                                                                                                      aCGA‑CATTCTGCACTT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCTGCTCGGTGTCTTTGTGCTGGCACTGCTGGGCTGGATTTTCAGTAAGTCTCTgggg  >  2:477518/1‑251 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                              TGGCTCCGGCAACACCGTCTAACACCGCGCGTGCGGGCGGCATTGTGTTACCGATCATCAACAGCGTGGCGGTGGCTTTGGGGTCCGAACCGGAAAAAAGTCCGCGTCGTGTCGGACATTACCTGATGATGTCCATTTACATGGTCACCAAAACCACCAGCTATATGTTCTTTACCGCAATGGCGGGGAACATTCTGGCGCTGAAAATGATCAACGA‑CATTCTGCACCT‑GCAAATTAGCTGGGGTGGATGGGCGCTGGCAGCCGGATTGCCGGGCATCATTATGCTGCTGGTCACCCCGCTGGTGATTTACACCATGTATCCACCAGAAATTAAGAAGGTGGATAACAAAACCATCGCTAAAGCGGGCCTTGCCGAACTAGGACCGATGAAAATCCGCGAAAAAATGCTGCTCGGTGTCTTTGTGCTGGCGCTGCTGGGCTGGATTTTCAGTAAGTCTCTGGGG  >  NC_000913/802299‑802762

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑