Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,040,143	A→G	21.8%	V222V (GTA→GTG)	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,143	0	A	G	21.8%	104.3 / 5.9	53	V222V (GTA→GTG)	yedY	membrane‑anchored, periplasmic TMAO, DMSO reductase
Reads supporting (aligned to +/- strand):  ref base A (27/13);  new base G (2/10);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 2.70e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

TATGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GAT‑G‑A‑CC‑GTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAG  >  NC_000913/2039906‑2040378                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 tctGCGCTGCGTGTAAGCGTGGTCGATTGTTGTGCCGTGGATTGGTTTTCCGCCGCACAAAGTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGGTCCGGAACAGATTCCAGTCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAATGATTGCGTCTCGACGAAGCAATGCAGCCGCTCACACT‑GAT‑G‑A‑CC‑GTA‑g                                                                                                                                                                                                                                                <  2:693013/237‑1 (MQ=255)  aTGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GAT‑G‑A‑CC‑GTA‑g                                                                                                                                                                                                                                                <  1:1267266/238‑1 (MQ=255)  aTGCGCTGCGGGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GAT‑G‑A‑CC‑GTA‑g                                                                                                                                                                                                                                                <  2:1163559/238‑1 (MQ=255)             ggAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GAT‑G‑A‑CC‑GTA‑GGTGTTTATGGc                                                                                                                                                                                                                                     <  1:329958/238‑1 (MQ=255)               aaGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GAT‑G‑A‑CC‑GTA‑GGTGTTTATGGCaa                                                                                                                                                                                                                                   <  1:59380/238‑1 (MQ=255)                    tGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGTGAAATATACTTTTTTTGAAGGGTGTCGTCTCGGCGAAGCAAAGCATTCCCTCCCACT‑GAT‑G‑A‑GC‑GGA‑GGGGTGTTTGGGAAGGCGtt                                                                                                                                                                                                                             >  2:884065/1‑239 (MQ=255)                      gTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATTCTTATGTCGCAGGGTTTGGTCTCGGCGAAGCAAAGCATCCGCCCCCACT‑GAT‑G‑A‑CC‑GCA‑GGGGGTTATGGGAAgggggta                                                                                                                                                                                                                            >  1:551336/1‑233 (MQ=255)                                 tgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGTCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GAT‑G‑A‑CC‑GTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAAt                                                                                                                                                                                                              >  1:918817/1‑239 (MQ=255)                                            gtttttCCGCTGCAAAATTTTCTGGCGCTTCGCGACCCTACCAGAAAGGCGAAGTATGTCTCTGTCAAAACAATTAATGCCCCGGAACAGACCCCGGGCCAGCGGGGCCGTTTAAGCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATGGCGTCTCGACGAAGAAATGCACCCGCTCACACT‑GAT‑G‑A‑CC‑GTG‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGt                                                                                                                                                                                                   <  2:551336/237‑1 (MQ=255)                                               ttcccgcgccacaAATTGCTGGCGCTTGCGGAACCTCCCGGCAATGCGAAGTAGTTCGCTTTCGAAAAAATTTATGCACCGGAACAGATGCCAGGCCGGCAGGACCGCTTTATCGGGGGTGGGCTAAAATTTCTTTTTGTCGAAGGTTTGCGTCTCGACGAAGCAATGCACCCGCTCACACT‑GAT‑G‑A‑CC‑GTG‑GGTGTTTAGGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCg                                                                                                                                                                                               <  2:610713/230‑1 (MQ=255)                                                  ccGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGAAGGACCGCTTTATCGGCGGTGGGGGGCAAAATACTTATTTCGAAGGGGTGTGGCTCCACCGAGCAACGGCCTCGCTCAACCT‑CAT‑GCA‑CG‑GTG‑GGGGTTT‑TGGCAAGGGGTTACCGCCACCAAAAATGgcgccggggagggc                                                                                                                                                                                               >  2:23939/1‑229 (MQ=255)                                                       gCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATTCTTATGTCGGAGGATTGGGTCTCGGCGAAGCAATGCATCCGCCCACACT‑GAT‑G‑A‑CC‑GGA‑GGTGTTTATTGTCAGGGGGTAACGCCCC‑‑AAAATTGGGCCGCGGT‑GGGGCTGCTtg                                                                                                                                                                                       >  1:589058/1‑239 (MQ=255)                                                        cacaAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGGAAAATCCTTTTTTTGCAGGAGTTCGTCTCGGCCCAGCAAAGCAACCGCTCCCACC‑TAT‑G‑ATCC‑TGG‑GGGGTTTTTGGAAGGGGGTTACC‑CCCC‑‑AAAAAAGCGCGCCCGG‑GGGACTGATtgt                                                                                                                                                                                      >  2:876796/1‑239 (MQ=255)                                                          aaaatggtgggCGCTTGCAGACCCTCCCAGCAATGCGAAGTGTTCCGTTTTAAAAAAATTTTATGCCCCGGAAAAGATGCCAGGCCAGCAGGACCGTTTTATGGGCGGTGGTCTGAAATCTCCTTATGTCGAAGGTTGGCTTCTCGACGAAGAAATGCCTCCGCTCACACT‑GAT‑G‑A‑CG‑GGA‑GTTGTTTGGGGCAGGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGc                                                                                                                                                                                    <  2:640036/230‑1 (MQ=255)                                                                  tGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCCTTTTTTCGAAGGATTGTGTCTCGCCGAAGCAAAGCATCCGCCCACAAC‑TAT‑G‑A‑CC‑GTAGGGGGTTTTTGGAAAGG‑GTTATCCCCAC‑‑AAAATGGCGCGCCGGT‑GCGGCTGGTTGTGTCGTGG‑gaa                                                                                                                                                                           >  1:276482/1‑239 (MQ=255)                                                                     cgttggCCAAACCTCCCGGCAATGCGAAGTATGTCGCTTCCGAAAAAATTTATGCACCGGAACAGATGCCGGGCCAGCAGGACCGCTTTATCGGGGGTGGGCTGAAATATCCTAATGTCGAAGGATGGCGTCTCGACGAGCCATTGCACCCGCTCACACT‑GAT‑G‑C‑CC‑GTG‑GGTGTTTAGGGCAGGGCTTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAata                                                                                                                                                                         <  1:630057/233‑1 (MQ=255)                                                                     cgcTTGCCGAACCTACCAGCAATGGGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAAGGCATCCGCTCACACT‑GAT‑G‑A‑CC‑GTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAata                                                                                                                                                                         <  1:884065/238‑1 (MQ=255)                                                                                          aTGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGGAAAAATCTTTTGTCGGAGGGTTGTGTCTCGCCGGAGCAAATAATCCGCTCACCCT‑ATT‑T‑AGCG‑GTA‑GG‑GTGTTTTGCGAGGAGTTTCCGCCCC‑‑CAAAAAGCCCGCCGGT‑GCGCCTCGTTGTGGCGGGG‑CGATATCGCTTTTCAGAGATTaaaact                                                                                                                                                   >  2:1311350/1‑236 (MQ=255)                                                                                           tGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGTAAAATCCTTTTTTCGAAGGATTTGGGCTCGGCGGAAGAATGCGTACCCTCACACT‑CCT‑G‑A‑AC‑GTA‑GGGGTTTATGGGCAAGAGGTACCGCCCC‑‑ACAAAGGCGCGCCGGG‑GTGACGCATTTTTGCGTGG‑GAATATGGCCTTTAAGAGAGTAAATAGa                                                                                                                                                  >  2:391781/1‑239 (MQ=255)                                                                                              gaattttgtcttttgcGAAAAAATTTGTGCCCGGAAACAGGTGCCGGGCCGGTAGGACCGCTTTCGGGGCGGTGGGCTGAAATATCTTTATGTCGAAGGTGTGCGTCGCGACGAAGCAAGGCTCCCGCTCAAACG‑GGT‑G‑G‑CG‑GTG‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAg                                                                                                                                               <  2:304541/224‑1 (MQ=255)                                                                                                aGTATGCCGTTTTAGAAACATTTTATCCCCCGAAACAGATGCCGGGCCAGCAGGCCCGTTTTACGGGGGGGGGGATAAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCACACT‑GAT‑G‑A‑CC‑GTG‑GGTGTTTATGGAAAGGCTTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                                              <  1:1167738/238‑1 (MQ=255)                                                                                                     gtctctttgaaaatattgtATTCACGGGAACAGTTGCGAGGCCAGCAGGACGGATTGACGGCCGGGGGGAGTAAATATCCTAATGTCAAAGGATGGTTTCCCGACGAAGGAATGCCCTCTCGCACGCT‑GGT‑G‑A‑CT‑GTA‑GTTTTTTAGGGTAGGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTa                                                                                                                                         <  2:918817/220‑1 (MQ=255)                                                                                                     gTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGCTGGATAGAAAAAACTTATGGCGCAGGAGTGTGTCGCGACGGAGAAAAGACACCGCGCACACT‑GATGG‑A‑CA‑GT‑‑GGGGTTTATGGGAAGGAGGTGCGGCCAC‑‑AAAAAAGAGGGCCGGTGGCGA‑GGATTGTGCCGGGGAAAAAATGGGTTTAAGGGGGTaaaacagcgatggctat                                                                                                                                         >  2:357845/1‑226 (MQ=255)                                                                                                                             ccGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTTTTTCGAAGGATTTCGTCTCGACGAAGCAATGCATCCGCCCACACT‑GAT‑G‑A‑CC‑GTA‑GGGGTTTATGGGAAGGCGGTACCCCCCC‑‑AAAAAGGCGCGCCGGG‑GCGACTGATTTTGCCGTGG‑GAATATTGCTTTTAAGGGATTAAAACGATAGTCAGTATTAAACTGACCCGCGAGCGTCCGc                                                                                                                 >  1:1143458/1‑238 (MQ=255)                                                                                                                              gcaaagcgttgccgcgccagCGGGACCGTTTTACGGGGGGGGGGCGGAAATCCCTTTAGGGAGGAGGATGGCGTCTGCAGGAACCAATGCACCCGCAC‑CAAT‑TAC‑G‑G‑CCGGTG‑GGTGTTTTGGGCAAGGCGCTGCCGGCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑ACATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCa                                                                                                               <  2:595380/224‑1 (MQ=255)                                                                                                                              cGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGGAAAATTCTTATGTCGCAGGATTGTGTCTCGGCGAAACAATGCATCCGCCCACACT‑GAT‑G‑A‑CC‑GCA‑GGGGGTTATGGCAAGGGGTTACCCCCCC‑‑AAAAAAGCGCGCCGGG‑GCGGCGGGTTTTGGCGTGG‑GGAAAAGGCTTTTAAAGGGGTAAAAAAATAATAAGTATTAAGcaggacccccagcgtgccccc                                                                                                               >  2:622794/1‑219 (MQ=255)                                                                                                                                    agagGCCAGGCCAGCAGGCCCGTTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCACACT‑GAT‑G‑A‑CC‑GTA‑GGTGTTTAGGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAaccacc                                                                                                         <  1:1347964/239‑1 (MQ=255)                                                                                                                                      aTGCCAGGCCAGCAGGACCGCTTTATCGGCGGCGGGGTTGAATCAACTTATTTTGCAGGAGTGTGTCTCCCCGAAACAATGCATCCCCCCCCACT‑GCT‑G‑T‑CC‑GTA‑GGTGGTTGTGGCAAGGTGGTACCCCCCC‑‑AAAAACGCGAGGCCGG‑GGGGCTGATTCTGCCGTGG‑AGATCTTGGGTTTAAGGGTGTAAAACGAGAGGGAGGATTGAGTTGGCCCCCGAGctgccccccaccccccg                                                                                                       >  2:1350485/1‑223 (MQ=255)                                                                                                                                       gtgcaggacagAAGGACCGATTTCGCGGGGGGGGGCTGAAATATCTTTATGTGGAAGGATGGCGATGCGAGAAAGAAATGCACCCGCTCA‑ACTGGAT‑G‑A‑CT‑GTG‑GTTGTTTATGGCAGGGCTTAACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTg                                                                                                       <  2:1149170/235‑1 (MQ=255)                                                                                                                                               aagcggGACCGTTTTATCGGGGGGGGGTTAAAATATCCTTATGTGGAAGGATTGCGTCTCGCAGAAGAAATGCACCCGCTCACACG‑GAT‑G‑A‑CC‑GTG‑GGTGTTTATGGCAGGGCTTTCCCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGc                                                                                              <  2:516168/238‑1 (MQ=255)                                                                                                                                                                cggcggGGGGCTGAAATATCCTTATGGCGAAGGATTGCGTCTCGACGAAGCAATGCCTCCGCTCACCCG‑GAT‑G‑A‑CG‑GTA‑GGTGTTTAGGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATAc                                                                              <  2:1143458/238‑1 (MQ=255)                                                                                                                                                                           tGAAATATCTTTATGCCGAAGGATTGCGTCTCGACGAAGCAAGCCACCCGCTCACACT‑GAT‑G‑C‑CC‑GTA‑GGTGTTTATGGAAGGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGc                                                                   <  1:1350485/238‑1 (MQ=255)                                                                                                                                                                           tGAAATACCTTTATGGGGAAGGTGTGCGTCGCGACGAAGCAATGACTCCGCTCACACT‑GAT‑G‑C‑CG‑GTG‑GGTGTTTAGGGCAAGGCTTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGc                                                                   <  2:589058/238‑1 (MQ=255)                                                                                                                                                                                 atCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GAT‑G‑A‑CC‑GTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGtt                                                            <  1:421143/239‑1 (MQ=255)                                                                                                                                                                                       aTGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GAT‑G‑A‑CC‑GTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTc                                                                                                                                             >  1:386309/1‑152 (MQ=255)                                                                                                                                                                                       aTGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GAT‑G‑A‑CC‑GTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTg                                                                                                       >  1:165025/1‑190 (MQ=255)                                                                                                                                                                                       aTGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GAT‑G‑A‑CC‑GTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCg                                                      >  1:466371/1‑239 (MQ=255)                                                                                                                                                                                       aTGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GAT‑G‑A‑CC‑GTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCg                                                      >  1:414718/1‑239 (MQ=255)                                                                                                                                                                                       aTGTCGAAGGATTGCGTCTCGAAGAAGCAATGCATCCGCTCACACT‑GAT‑G‑A‑CC‑GTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCg                                                      >  2:468618/1‑239 (MQ=255)                                                                                                                                                                                            gAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GAT‑G‑A‑CC‑GTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTCAGCTGACCCGCGAGCGTCCGCCACCCACCTGGAATCTGGCAGCGCCTGACGAATACGGATTTTACGCCAATGTTACTACGTATGt                                                 >  1:1090895/1‑239 (MQ=255)                                                                                                                                                                                             aaGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GAT‑G‑A‑CC‑GTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGtt                                                >  2:1138231/1‑239 (MQ=255)                                                                                                                                                                                              aGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GAT‑G‑A‑CC‑GTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTg                                               >  2:236384/1‑239 (MQ=255)                                                                                                                                                                                                         tcaccgaACCAATGCATCCGCTCACACT‑GAT‑G‑A‑CC‑GTG‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCgcgc                                    <  1:658134/235‑1 (MQ=255)                                                                                                                                                                                                                 gCAATGCATCCGCTCACACT‑GAT‑G‑A‑CC‑GTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTcaca                            >  1:1205098/1‑239 (MQ=255)                                                                                                                                                                                                                     tGCATCCGCTCACACT‑GAT‑G‑A‑CC‑GTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCAGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCt                        >  2:478725/1‑239 (MQ=255)                                                                                                                                                                                                                          ccGCTCACACT‑GAT‑G‑A‑CC‑GTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGACTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATAAGGTTTTTACGCCAACGTTAATCCGTATGTTGATCAACCGAGCTGGTCACAGGCGAcctc                   >  2:867746/1‑237 (MQ=255)                                                                                                                                                                                                                               cacaCT‑GAT‑G‑A‑CC‑GTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGAt              >  1:320423/1‑239 (MQ=255)                                                                                                                                                                                                                                  acT‑GAT‑G‑A‑CC‑GTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAg                                                                                                                                    >  1:1232066/1‑118 (MQ=255)                                                                                                                                                                                                                                    cggat‑g‑a‑CC‑GTA‑GGTGTTTAGGGAAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTAt          <  1:391781/237‑1 (MQ=255)                                                                                                                                                                                                                                       at‑g‑a‑CC‑GTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTAACGAACGATTTATTgg       >  2:548441/1‑239 (MQ=255)                                                                                                                                                                                                                                              cc‑GTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCa   >  2:590627/1‑239 (MQ=255)                                                                                                                                                                                                                                                 gTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAg  >  1:116021/1‑238 (MQ=255)                                                                                                                                                                                                                                                 gTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAg  >  1:519752/1‑238 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 TATGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GAT‑G‑A‑CC‑GTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGT‑GCGACTGATTGTGCCGTGG‑AAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAG  >  NC_000913/2039906‑2040378

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: ‑