Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	805,325	T→C	5.7%	V608A (GTG→GCG)	ybhJ →	aconitase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	805,325	0	T	C	5.7%	266.8 / 4.5	85	V608A (GTG→GCG)	ybhJ	aconitase family protein
Reads supporting (aligned to +/- strand):  ref base T (23/56);  new base C (2/3);  total (25/60)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGG  >  NC_000913/805191‑805463                                                                                                                                       |                                                                                                                                           ttATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACCGAGGTGTTTg                                                                                                                                       >  2:1649771/1‑140 (MQ=255)   aTCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAACGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTgcg                                                                                                                                     <  2:1222882/140‑1 (MQ=255)   aTCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTgcg                                                                                                                                     >  2:1269076/1‑140 (MQ=255)   aTCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTgcg                                                                                                                                     <  2:713087/140‑1 (MQ=255)         cAATTCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCGGCGATCCCGTTTATGTTAGCAGAAGTAAAGCGTCTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCAtt                                                                                                                               <  2:1361721/140‑1 (MQ=255)         cAAATCCGATTGGTCTGGCGGAGTTTACCCTTTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCAtt                                                                                                                               <  2:1861203/140‑1 (MQ=255)         cAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCGGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCAtt                                                                                                                               <  2:591474/140‑1 (MQ=255)         cAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCAtt                                                                                                                               <  2:552799/140‑1 (MQ=255)         cAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCAtt                                                                                                                               <  2:514242/140‑1 (MQ=255)         cAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCAtt                                                                                                                               <  2:472915/140‑1 (MQ=255)         cAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCAtt                                                                                                                               <  1:690417/140‑1 (MQ=255)         cAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCAtt                                                                                                                               <  2:1461330/140‑1 (MQ=255)         cAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCAtt                                                                                                                               <  1:1854266/140‑1 (MQ=255)         cAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCAtt                                                                                                                               <  2:1696023/140‑1 (MQ=255)         cAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCAtt                                                                                                                               <  2:1798182/140‑1 (MQ=255)         cAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCAtt                                                                                                                               <  2:1834650/140‑1 (MQ=255)         cAAATCAGATTGGTCTGGCGGAGTTTTCCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAATATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCAtt                                                                                                                               <  2:830760/140‑1 (MQ=255)         cAAAACCGATTGGTCTGGCGGAGTTAACCCTGACTCCCCGCGATCCCGGTAATGTTAGCACAAGTAAAGCGACTGCTGAGCTGCAAAATCAGCCTCTGTCGGGGAATGTCAGCGAGCTGACAGAGCTCATAGCGCGCAtt                                                                                                                               <  2:573345/140‑1 (MQ=255)                 aTTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGAt                                                                                                                       <  2:327627/140‑1 (MQ=255)                 aTTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGAt                                                                                                                       >  1:1562294/1‑140 (MQ=255)                 aTTGGTCTCGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGAt                                                                                                                       >  1:372866/1‑140 (MQ=255)                   tGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTg                                                                                                                     <  2:1634986/140‑1 (MQ=255)                   tGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTg                                                                                                                     <  2:1695525/140‑1 (MQ=255)                         ggCGGAGTTTACCATGTCTCGCCGCGTTCCCGGTTATGTTTGCAGTTGTAAAGCGTCTGATGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCTGATTGCTGGTc                                                                                                               <  2:1709418/140‑1 (MQ=255)                           cGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAg                                                                                                             <  2:637034/140‑1 (MQ=255)                             gattttACCGTGTATTGCCGCGATGCCGTTTATGTTTGCTGAAGTAAAGCGACTGCTGTTCTTGAATATCAGAGTCTGGCGGGGATTGTCAGCTAGATGACTTAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCTGGa                                                                                                           <  2:809532/137‑1 (MQ=255)                             gAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGa                                                                                                           <  1:116394/140‑1 (MQ=255)                             gAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGa                                                                                                           <  1:116389/140‑1 (MQ=255)                             gAGTTTACCCGGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGa                                                                                                           <  2:1846098/140‑1 (MQ=255)                              aGTTTACGCTGTCTCGCCGCGATCCCGCTTATGTTAGCAGAAGTAAAGCGTCGGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAg                                                                                                          <  2:1553943/140‑1 (MQ=255)                                tttACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCa                                                                                                        <  2:1025545/140‑1 (MQ=255)                                tttACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCa                                                                                                        >  1:1346700/1‑140 (MQ=255)                                tttACCCTGTCTCGCCGCGATCCCGGGTAGGTTAGCAGTCGTAAAGCGACTGCTGAGCGGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGTGCTGGCTGAGGTGTTTGCGCGCTTTAAGCAGATTGCTGGTCAGGAGCa                                                                                                        >  1:798881/1‑140 (MQ=255)                                  tACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCata                                                                                                      >  2:853042/1‑140 (MQ=255)                                       tGTCTCGCCGCGATCCGGGTTATGTTCGCAGACGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGAt                                                                                                 <  1:1660374/140‑1 (MQ=255)                                       tGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGAt                                                                                                 <  1:1763788/140‑1 (MQ=255)                                       tGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGAt                                                                                                 <  1:1763686/140‑1 (MQ=255)                                        gTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATc                                                                                                >  2:862163/1‑140 (MQ=255)                                        gTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATc                                                                                                >  2:1827119/1‑140 (MQ=255)                                        gTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGCGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAg                                                                                                          >  1:1556298/1‑130 (MQ=255)                                        gTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGCGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAg                                                                                                          >  1:1556331/1‑130 (MQ=255)                                        gTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGCGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAg                                                                                                          <  2:1556298/130‑1 (MQ=255)                                        gTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGCGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAg                                                                                                          <  2:1556331/130‑1 (MQ=255)                                                   aTCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCATGAGCATATTGATCCGCTGCAAACt                                                                                     >  1:1078915/1‑140 (MQ=255)                                                      ccGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGaa                                                                                  <  2:319889/140‑1 (MQ=255)                                                           tATGTTAGCAGTAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTgg                                                                             <  1:1027504/140‑1 (MQ=255)                                                            aTGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgc                                                                            <  1:561362/140‑1 (MQ=255)                                                              cctAGCCGAACTAAAGCCCATGCTGAGCAGGAACATCAGCGACTGGCGGCGAAGGTCATCGAGCTGACAGAGGAGTTAGCGCGCAATAAGCAGATTGCTGGTCAGGAGCCTATTGATCCGCTGCAAACTGAACTGGgcag                                                                          <  2:260709/138‑1 (MQ=255)                                                               ttagatgAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCTCTGCAAACTGAAATTGgcagc                                                                         <  2:605364/134‑1 (MQ=255)                                                               ttAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcagc                                                                         <  2:556743/140‑1 (MQ=255)                                                               ttAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcagc                                                                         <  2:1828827/140‑1 (MQ=255)                                                               ttAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcagc                                                                         <  2:1369315/140‑1 (MQ=255)                                                               ttAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcagc                                                                         <  1:65874/140‑1 (MQ=255)                                                               ttAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcagc                                                                         <  1:826431/140‑1 (MQ=255)                                                               ttAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcagc                                                                         <  2:486338/140‑1 (MQ=255)                                                               ttAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcagc                                                                         <  2:486339/140‑1 (MQ=255)                                                                tAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGgcagca                                                                        <  2:1634076/140‑1 (MQ=255)                                                                 aGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCAt                                                                       <  2:1749615/140‑1 (MQ=255)                                                                       ccaaacgCGACAGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTa                                                                 <  1:1843181/134‑1 (MQ=255)                                                                        gTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTAt                                                                <  1:1203870/140‑1 (MQ=255)                                                                            aGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCgg                                                            >  1:933354/1‑140 (MQ=255)                                                                                aCTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCAtt                                                                                                                               >  1:484924/1‑69 (MQ=255)                                                                                aCTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCAtt                                                                                                                               <  2:484924/69‑1 (MQ=255)                                                                                aCTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGaa                                                        >  1:25898/1‑140 (MQ=255)                                                                                aCTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGaa                                                        >  1:488662/1‑140 (MQ=255)                                                                                   gctgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAAcc                                                     >  1:1850707/1‑140 (MQ=255)                                                                                    ctgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCa                                                    <  2:955632/140‑1 (MQ=255)                                                                                    ctgAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCa                                                    <  1:399097/140‑1 (MQ=255)                                                                                      gAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAgg                                                  <  1:1563754/140‑1 (MQ=255)                                                                                         cTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCACGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGa                                               <  2:1078915/140‑1 (MQ=255)                                                                                              aaaTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGtt                                          >  2:219371/1‑140 (MQ=255)                                                                                               aaTCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTc                                         <  2:146230/140‑1 (MQ=255)                                                                                                         cggCGGGGAAAGGCAGCGAGCTGACAGAGGCGTTTGCGCGCATTAAGCAGAGTGCTGGACAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGaa                               <  2:802181/139‑1 (MQ=255)                                                                                                               ggAATGTCAGCGAGCTGACAGAGGTGTTGGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCTGGCGATGGTTCTGCGCGTGACCAggcg                         >  2:250512/1‑140 (MQ=255)                                                                                                                     tCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGACCAGGCGGCGAGc                   >  2:113212/1‑140 (MQ=255)                                                                                                                     tCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGc                   >  2:168085/1‑140 (MQ=255)                                                                                                                        gCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGc                <  2:1197920/140‑1 (MQ=255)                                                                                                                        gCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGc                <  2:1198114/140‑1 (MQ=255)                                                                                                                         cGAGCTGACAGAGGTCTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCAGCAAACTGAAATTGGCAGCATGGTCTTTGCGGTGAAACCAGGCGATGGTTCTGCTCGTGAACAGGCGGCGAGCTGcc               <  2:543694/140‑1 (MQ=255)                                                                                                                           aGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCTAGCTGCCAg             >  2:1840174/1‑140 (MQ=255)                                                                                                                           aGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCTAGCTGCCAg             >  2:1922963/1‑140 (MQ=255)                                                                                                                           aGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCTAGCTGCCAg             >  2:1922805/1‑140 (MQ=255)                                                                                                                           aGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCTAGCTGCCAg             >  2:1840161/1‑140 (MQ=255)                                                                                                                            gCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGc            <  1:298179/140‑1 (MQ=255)                                                                                                                                   gagGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGAt     <  1:952491/140‑1 (MQ=255)                                                                                                                                      gtgtTTGCGCGCATTTAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTgg  <  2:618364/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                                           TTATCGTTCAAATCCGATTGGTCTGGCGGAGTTTACCCTGTCTCGCCGCGATCCCGGTTATGTTAGCAGAAGTAAAGCGACTGCTGAGCTGGAAAATCAGCGTCTGGCGGGGAATGTCAGCGAGCTGACAGAGGTGTTTGCGCGCATTAAGCAGATTGCTGGTCAGGAGCATATTGATCCGCTGCAAACTGAAATTGGCAGCATGGTCTATGCGGTGAAACCAGGCGATGGTTCTGCGCGTGAACAGGCGGCGAGCTGCCAGCGTGTGATTGG  >  NC_000913/805191‑805463

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