Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,040,157	A→G	18.4%	K227R (AAG→AGG)	yedY →	membrane‑anchored, periplasmic TMAO, DMSO reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,040,157	0	A	G	18.4%	130.8 / 4.0	68	K227R (AAG→AGG)	yedY	membrane‑anchored, periplasmic TMAO, DMSO reductase
Reads supporting (aligned to +/- strand):  ref base A (29/25);  new base G (2/12);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.41e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGC‑GGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGACCGTA‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGAT  >  NC_000913/2039920‑2040392                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               aaGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGC‑GGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGACCGTA‑GG‑TGTTTATGGCaa                                                                                                                                                                                                                                                <  1:1157343/238‑1 (MQ=255)        gTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGC‑GGAGGGCGGTAAAAATCTTATGTCTTAGGATTCTGGTTCGACGGAGC‑AAGGAAGCCGCCCACACCTTTGAACGTA‑GG‑TGGTTATTGCTGGGCGGTa                                                                                                                                                                                                                                         >  2:461622/1‑238 (MQ=255)           gATGGTGGTGCCGTGGGTTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGC‑GGTGGGGTGAAATATCCTTTTGTCGAAGGAGTGCGGCTCGCCGAAGC‑AAAGCATCCGCTCACACTGATGACCGTT‑GG‑GGTTTTTTGCAAAGCggttccc                                                                                                                                                                                                                                      >  1:553398/1‑232 (MQ=255)                       gTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGC‑GGTGGGGTGGAAAATCCTTATTTTGGAGGAGTGCGTCTCGGCGGAGG‑AAAGCATCCGCTCAAACTGATTACCGCA‑GG‑TGTTTTTGGCAAAGCGTTACCGCCACaaaaaagcg                                                                                                                                                                                                                         >  1:1352871/1‑234 (MQ=255)                          gATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGC‑GGTGGGGTGGAATATTCTTATGTCGAAGGAATGCGGCTCGACGCCAC‑AAAGCCTCCGCTCCCAATGATGAACCTA‑GG‑TGTTTTTGGGCAGGCGTTACCGCCACAAAAAGgcggg                                                                                                                                                                                                                       >  1:936689/1‑236 (MQ=255)                             tGGTTTTCCGCTGCACAATTTGCTGGCGCTTGCCGAACCTCCCAGAAAGGCAAAGTATGTCGTTTTCGAAACAATTTATGCCCCGGAACAGATGCCAGGCAAGGAGGACCGTTTTATCGGG‑GGGGGGCTAAAATATCCTTATGTCGAAGGATTGCGTCCCGACGAAGC‑AATGCATCCGCTCACACTGATGCCCGTG‑GG‑TGTTTATGGCAAGGCTTTACCGCCACAAAATGGCGC‑GCCgg                                                                                                                                                                                                                  <  2:553398/239‑1 (MQ=255)                                  ttCCGCTGCACAAATGGCTGGCGCTGCCCGAACCAACCGGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCGGCAGGCCCGCTTTATCGGC‑GGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGCCCGTA‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGa                                                                                                                                                                                                             <  1:1007656/239‑1 (MQ=255)                                         gCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGCGGAATAATCTTTAT‑TTGAAGGATTGTGTCTCGTCGAAGCAAATGAATCCG‑TCACCCCGTTGACCGGC‑GGACGTTTTTGGC‑AGGCGGTACCGCCCCacaaagagcggcccgctgggccgggtg                                                                                                                                                                                                       >  2:1444622/1‑216 (MQ=255)                                           acaAATTGCTGGCGCTTGCCGAACCTACCAGAAAGGCGAAGTATGTCGTTTGCGAAACAATTTATGCACCGGAACAGATCCCAGGCCAGCAGGACCGCTTTATCGGC‑GGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGACCGTA‑GG‑TGTTTATGGAAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATtgtg                                                                                                                                                                                                    <  1:1429148/239‑1 (MQ=255)                                           aaaatttGCTGGCGCTGCCCAACCCTCCCAGCAATGCGAAGTATGTCGTTTTCAAAAAAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGTTTTAGGGGG‑GGTGGGCTAAAATATCCTTGTGTCAAAGGATTGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGACCGTA‑GG‑TGTTTAGGGCAGGGGGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATtgtg                                                                                                                                                                                                    <  2:983481/234‑1 (MQ=255)                                           aaaaaTTGCTGCCCCTTGCCGAACCACCCAGCAAGGCGAAGGAGGTCGCTTTCGAAACAATTTATGCACCGGAAAAGATGCCAGGCCAGCAGGACCGCTTTACGGGC‑GGGGGGCTGAAATTTCCTTATGTCGAAGGATGGCGTCTCGACGAAGA‑AATGCATCCGCTCACACTGATGACCGTG‑GG‑TGTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATtgtg                                                                                                                                                                                                    <  2:11378/237‑1 (MQ=255)                                                    tGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGC‑GGTGGGGTGGAATAACCTTTTGTCGAAGGGGTGCGTCTCCGCGAAAC‑AAAGCAACCGCTCCCCCTGCTGACCGGAGGG‑GGTTTTTTGCAGGGGGTAACC‑CCCCAAAAAAGCGC‑CCCGGGGCGGCGGATtgtggccggggga                                                                                                                                                                                           >  2:616006/1‑230 (MQ=255)                                                               gAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACGGGAACAGATCCCAGGCCAGCAGGACCGCTTTATCGGC‑GGTGGGCTGAAATATCCTTATGTCGAGGGTTGGCGTCTCGACGAAGA‑AATGCACCCGCTCACACTGATGACGGTA‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTa‑a                                                                                                                                                                              <  2:570277/239‑1 (MQ=255)                                                               gAACCTACCAGAAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGC‑GGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGC‑AATGCACCCGCTCACACTGATGACCGTA‑GG‑TGTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTa‑a                                                                                                                                                                              <  1:616006/239‑1 (MQ=255)                                                               aaaCCTACAAGCATGGAGAAGTATGTCGCTTTCAAAAAAATTAATGCACGGGAACAGATGCAGGGCCAGCAGACCCTCTTTACCGGG‑GGTGGGCTAAAATATCCTTATGTCAAAGGATTGCGTCTCGACGAAGA‑AATGCATCCGCTCACACGGAGGACCGTG‑GG‑TGTTTATGGCAGGGGGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTa‑a                                                                                                                                                                              <  2:549915/238‑1 (MQ=255)                                                                       cAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATCCCCCGGAAAAGATGCCAGGCCAGCAGGACCGTTTTATCGGC‑GGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGACCGTG‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTa                                                                                                                                                                      <  1:944292/239‑1 (MQ=255)                                                                             tGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGC‑GGTGGGGTGAAAAATCCTTATTTCGAAGGATTTCGGCTCGCCGAAGC‑AAAGCATCCGCGCACACTGATGACCCCG‑GG‑GGGTTTTGGGAAAGGGGTAACCCCCCAAAAAAGGGG‑GCCGGGGGGACTGAGTTTGCCGGGGGAAAAATGGC‑TTA‑AAGGGGGTAAATAGa                                                                                                                                                                >  2:676972/1‑239 (MQ=255)                                                                             tGCGAAGTATGTCGCTTTCAAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGTTTTATCGGC‑GGTGGGCTAAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGCCCGTG‑GG‑TGTTTATGGAAAGGCTTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGa                                                                                                                                                                <  1:1529007/239‑1 (MQ=255)                                                                               cGAAGTATGTCGCTTTCGAAAAAATTTATGCCCCGAAACAGATGCCAGGCCAGCAGGCCCGTTTTATCGGC‑GGTGGGCTAAAATATCCTTAGGTCGAAGGATTGCGTCTCGACGAAGC‑AATGCACCCGCTCACACTGATGCCCGTA‑GG‑TGTTTAGGGAAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGa                                                                                                                                                                <  1:213262/237‑1 (MQ=255)                                                                                gAAGTATGCCGGTTCCGAAATAATTTATCCACCGAAACAGATGCCAGGGCAGCAGGACCGGTTTATGGGC‑GGGGGCCTGAAATTTCTTTATGCGGAAGGATGGCGTCTGGACGAAGA‑AATGCACCCGCTCACACTGAGGACCGTG‑GG‑TGTTTATGGAAGGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATa                                                                                                                                                              <  2:644045/238‑1 (MQ=255)                                                                                 aaGTGTGTGGCTTTAAAACAAATTGTTGCCCCGGAAAAGAGGCCAGGCCAGCAGGCCCCTTTTAGCGGC‑GGGGGGTTAAAAAATCCTTGTGTCGAGGGTTTGGGTCTCGAGGAAGA‑AATGCACCCGCTCACTCTGATGACCGTG‑GG‑TTTTTAGGGCAGGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGt                                                                                                                                                            <  2:1359363/239‑1 (MQ=255)                                                                                  agtattccgtgttgaaaaaaaTGTATGCCCCGGAAAAGATGCCAGGCAAGCAGGCCCGTTTTAGCGGC‑GGTGGGTGAAAATATCTTTGTGTCGAGGGATTGCGTCTCGACGAAGC‑AATGCACCCGCTCACACTGATGACCGTG‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGt                                                                                                                                                            <  2:936640/223‑1 (MQ=255)                                                                                  aGTATGTCGTTTCCAAAACATTTTTGGCCCCGAAAAAGATGCCGGGCCAGCAGGACCGTTTTATCGGC‑GGTGGGTTAAAATATCCTTATGTCGAAGGAGTGCGTCTCGACGAAGC‑AATGCACCCGCTCACACTGATGACCGTG‑GG‑TGTTTATGGAAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGt                                                                                                                                                            <  2:686675/238‑1 (MQ=255)                                                                                    tATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGC‑GGGGGGGTACAATATACATATTGATGAGGATTGAGTCTCGGCCGAAC‑AAAGCATCCGCACACACTCAAGACCGCA‑GG‑GGGTTATGGGAAGGCGGTACAGCCAAACAAAGGCGCGGCAGGTGGGGTTG‑TTGTGGC‑GTGGGAATATGGGTTTACAAgggtgagaaacaatgtcat                                                                                                                                                          >  2:37772/1‑223 (MQ=255)                                                                                      tcgcccttcggAAACAATTTAGGCACGGGAACAGACGCAAGGCCCGAAGCTCCGTTTTATCGGG‑GGTGGGCTGAAATACCCTTATGTCGAAGGATTGCGTCTCGACGAAGC‑AGTGCTCCCGCCCACACTGATGACGGTG‑GG‑TGTTTAGGGCAGGGTGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTa                                                                                                                                                       <  2:759338/229‑1 (MQ=255)                                                                                            tttCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGC‑GGTGGGGTGAAATATACTTTAGTCGAAGGAGTGCGTCTCGACGAAGC‑AATGCATTCGCCCACACTGATGAGCGTA‑GG‑TGGTTATTGCAAAGGGGTACAGCCCCACAAAGGGGC‑GCCGGGGCGGCTGGGTGTGTC‑GGGGGAATATGGCTTTT‑AAGGGGTTAAAAAGATAGACAGAATTaa                                                                                                                                                   >  2:448727/1‑237 (MQ=255)                                                                                            tttCAAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCGGGACCGTTTTACCGGC‑GGTGGGTTGAAATATCCTTATGTCGAAGGATGGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGACCGTA‑GG‑TGTTTATGGCAAGGCTTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGc                                                                                                                                                 <  1:1224109/239‑1 (MQ=255)                                                                                              tCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGC‑GGTGGGCTGCAAAATTCTTATTTCGAAAGGTTTCCTTTCGACGGAGC‑AATGCCTCCGCCCACACTGATGACCGTT‑GG‑GGGTTTTTGCAAGGCGTTACCACCACCAAAACGCGC‑GCCGGGGGGGATGATTGTGCC‑GTGGGAATATGGGTTTA‑AAAGGAGTAAATAGCTAATCAGTGTTaaattt                                                                                                                                               >  2:764446/1‑235 (MQ=255)                                                                                                       tttATGCCCCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGC‑GGGGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGACCGTA‑GG‑TGTTTATGGAAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGt                                                                                                                                                            <  1:543736/217‑1 (MQ=255)                                                                                                                cGGAACAGATGCCAGGACAGCAGGTCCGCTTTATCGGC‑GGTGGGCTGAAATACCCTTATGTCGAAGGATTGCGTCTCGACGAAGC‑AATGCACCCGCTCACCCTGATGACCGTA‑GG‑TGTTTATGGCAAGGCTTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGt                                                                                                                                                            <  1:739785/208‑1 (MQ=255)                                                                                                                 ggAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGC‑GGTGGGGTGGAAAATCCTTTTGTCGAAGGATTGCGTCTCGACGAAGC‑AATGCATCCGCTCCCACTGATGGCCGTA‑GG‑TGTTTTTGGCAAGGCGGTACCCCCCCAAAATGGGGC‑GCCGCGGGGACTGGTTGTTGC‑GGGGCAAAATTGCTTTT‑ATAGGGGTAAAACGATAGTCAGTATTAAGcctgtctctttttcacatctt                                                                                                                             >  2:543736/1‑218 (MQ=255)                                                                                                                   aaCAGATGCCAGGCCAGCAGGACCGTTTTATCGGC‑GGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGCCCGTA‑GG‑TGTTTATGGCAGGGCTTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAac                                                                                                                           <  1:1263253/238‑1 (MQ=255)                                                                                                                         gctgggcccagcgggcgctcTTTAGCGGC‑GGTGGGTCAACATCTCCTTATGGAGAAGGATGGCGTCTGGAGAAAGC‑AATGCATCCGCCCACACGGATGACGGGA‑GG‑TGTTTAGGGCAGGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTg                                                                                                                     <  2:156162/219‑1 (MQ=255)                                                                                                                                caggcagGACCGCTTTTCCGGG‑GGTGGGCTGAAATATCTTTATGCCGAAGGTTTGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGACCGTA‑GG‑TGTTTATGGCAAGGCTTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTgg                                                                                                             <  1:936154/236‑1 (MQ=255)                                                                                                                                    cagGACCGCTTTATCGGC‑GGTGGGGTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGA‑AATGCATCCCCTCCCCCTGATGACCGCA‑GG‑GGTTTTTTGCAAGGCGGTACCGCCCCCAAAAGGCGC‑GCCCGGGCGGCCGCTTGTGCC‑GTGGGAATATTGCTTTT‑AAGGGGTTAAATCGATAGGTAGTATTAAACTGACCCGCGCGCGTCCCCCCACCACCCGGAATCttggtgc                                                                                                         >  2:480151/1‑233 (MQ=255)                                                                                                                                          cgttttACCGGG‑GGGGGTCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACAAAGC‑AATGCACCCGCTCACACTGATGCCCGTG‑GG‑TGTTTATGAAAAGCGGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGa                                                                                                   <  1:1076157/236‑1 (MQ=255)                                                                                                                                                aagtggtgtaaaagagcaaaatCCTTATGGCGAAGGATGTCGTCTCCACGAACC‑AATGCCCCCCCCCCCCCTGAGCACCGTG‑GG‑TTTTTATGGAAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGaa                                                                                               <  2:857365/218‑1 (MQ=255)                                                                                                                                                         tGGGCTAAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGACCGTA‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGtttt                                                                                      <  1:545712/238‑1 (MQ=255)                                                                                                                                                                    atCCTTATGTGGAAGGATTGCTTCCCGACGAAGA‑AATGCATCCGCTCACACTGATGACGGTA‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGtt                                                                          <  1:461622/239‑1 (MQ=255)                                                                                                                                                                    atCCTTATGTCGAAGGATGGGGCCTCGAAGAGGC‑AATGCACCCGCTCACCCGGATGGCGGTG‑GG‑TGTTTTGGGCAGGGCTTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGtt                                                                          <  2:936689/239‑1 (MQ=255)                                                                                                                                                                    acccTTATTGCGAAGGATTGCGTCGCGACGAAGC‑AATGCATCCGCTCACACTGATGACCGTA‑GG‑TGTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGtt                                                                          <  1:480151/237‑1 (MQ=255)                                                                                                                                                                     tCCTTATGTCGAAGGATTGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGACCGTA‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTa                                                                         >  1:1060553/1‑239 (MQ=255)                                                                                                                                                                     tCCTTATGTCGAAGGATTGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGACCGTA‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTa                                                                         >  2:37631/1‑239 (MQ=255)                                                                                                                                                                     tCCTTATGTCGAAGGATTGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGACCGTA‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTa                                                                         >  2:1358750/1‑239 (MQ=255)                                                                                                                                                                       ctttggtcgagggttggggctccgacgaAGC‑AATGCACCCGCTCACGCTGATGACCGAG‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAAt                                                                       <  2:175662/217‑1 (MQ=255)                                                                                                                                                                       cTTATGTCGAGGGATTGCGTCCCGACGAAGC‑AATGCACCCGCTCACACTGATGACCGGA‑GG‑TTTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAAt                                                                       <  2:239051/239‑1 (MQ=255)                                                                                                                                                                       cTTATGTCGAAGGATTGCGTCCCGACGAAGC‑AATGCATCCGCTCACACTGATGACCGTA‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAAt                                                                       <  1:764446/239‑1 (MQ=255)                                                                                                                                                                          aTGTCGAAGGATTGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGACCGTA‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCg                                                                    >  2:820549/1‑239 (MQ=255)                                                                                                                                                                             tCGAAGGATTGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGACCGTA‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTAt                                                                 >  1:543419/1‑239 (MQ=255)                                                                                                                                                                               gAAGGCTTGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGACCGTA‑GG‑AGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGAGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCGGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGTAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGt                                                               >  2:572198/1‑239 (MQ=255)                                                                                                                                                                               gAAGGATTGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGACCGTA‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGttttt                                                                                     >  1:857365/1‑217 (MQ=255)                                                                                                                                                                                aaGGATTGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGACCGTA‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGtt                                                              >  1:99378/1‑239 (MQ=255)                                                                                                                                                                                      tGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGACCGTA‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCa                                                         >  2:1425473/1‑238 (MQ=255)                                                                                                                                                                                           ctggaggaAGC‑AATGCATCCGCTCACACGGATGACGGTG‑GG‑TGTTTAGGGCAGGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCcgc                                                    <  2:329873/232‑1 (MQ=255)                                                                                                                                                                                                       aaTGCATCCGCTCACAGTGATGACCGTA‑GG‑TGTTTATGGCAGGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTcaca                                          <  1:1351161/237‑1 (MQ=255)                                                                                                                                                                                                       aaTGCACCCGCCCACACTGATGCCGGTG‑GG‑TGTTTAGGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATc                                                                      <  2:785187/209‑1 (MQ=255)                                                                                                                                                                                                              ccGCTCACACTGATGACCGTA‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGa                                 >  1:635643/1‑239 (MQ=255)                                                                                                                                                                                                                 cTCACACTGATGACCGTA‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGa                                                            >  1:785187/1‑209 (MQ=255)                                                                                                                                                                                                                      acTGTTGACCGTA‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGATTGATTGTGTT‑GTGTAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGCCCGCCAACCACCTGGAATGTGTCAGCGCCTGACGAATACGGGTTTTACGCCAACGTTAATTCGTTTGCTGATCTACCGCGCTGGACACCGGCTTCAGTACGATTTa                         >  2:1178374/1‑239 (MQ=255)                                                                                                                                                                                                                                gTA‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAg                >  2:1056412/1‑238 (MQ=255)                                                                                                                                                                                                                                gTA‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAg                >  1:451451/1‑238 (MQ=255)                                                                                                                                                                                                                                    tgttttttAGGACAGGGCGTAACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAggc              <  2:101421/233‑1 (MQ=255)                                                                                                                                                                                                                                     g‑tgtTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAggcggc           >  2:1349740/1‑239 (MQ=255)                                                                                                                                                                                                                                     g‑tgtTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGAGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGAATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAggcggc           >  1:950138/1‑239 (MQ=255)                                                                                                                                                                                                                                        ttttATGGAAGGAGGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCAt         <  2:35808/238‑1 (MQ=255)                                                                                                                                                                                                                                           tATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGGTAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTAACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATcc       >  2:248011/1‑238 (MQ=255)                                                                                                                                                                                                                                             tGGAAAGGTGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCg    <  1:874395/239‑1 (MQ=255)                                                                                                                                                                                                                                              ggCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGa   >  2:1075482/1‑239 (MQ=255)                                                                                                                                                                                                                                               gCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGGTTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGAt  >  2:789088/1‑239 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                               AAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGC‑GGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGC‑AATGCATCCGCTCACACTGATGACCGTA‑GG‑TGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGC‑GCCGGTGCGACTGATTGTGCC‑GTGGAAATATGGCTTTA‑AAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGAT  >  NC_000913/2039920‑2040392

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