Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,040,143	A→G	19.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	19.8%	104.1 / 2.3	47	V222V (GTA→GTG)	yedY	membrane‑anchored, periplasmic TMAO, DMSO reductase
Reads supporting (aligned to +/- strand):  ref base A (25/12);  new base G (2/8);  total (27/20)
Fisher's exact test for biased strand distribution p-value = 1.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TATGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GATG‑A‑CCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAG  >  NC_000913/2039906‑2040378                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               tctGCGCTGCGTGTAAGCGTGGTCGATTGTTGTGCCGTGGATTGGTTTTCCGCCGCACAAAGTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGGTCCGGAACAGATTCCAGTCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAATGATTGCGTCTCGACGAAGCAATGCAGCCGCTCACACT‑GATG‑A‑CCGTA‑g                                                                                                                                                                                                                                              <  2:693013/237‑1 (MQ=255)  aTGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GATG‑A‑CCGTA‑g                                                                                                                                                                                                                                              <  1:1267266/238‑1 (MQ=255)  aTGCGCTGCGGGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GATG‑A‑CCGTA‑g                                                                                                                                                                                                                                              <  2:1163559/238‑1 (MQ=255)             ggAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GATG‑A‑CCGTA‑GGTGTTTATGGc                                                                                                                                                                                                                                   <  1:329958/238‑1 (MQ=255)               aaGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GATG‑A‑CCGTA‑GGTGTTTATGGCaa                                                                                                                                                                                                                                 <  1:59380/238‑1 (MQ=255)                    tGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGGGGGGTGAAATATACTTTTTTTGAAGGGTGTCGTCTCGGCGAAGCAAAGCATTCCCTCCCACT‑GATG‑A‑GCGGA‑GGGGTGTTTGGGAAGGCGtt                                                                                                                                                                                                                           >  2:884065/1‑239 (MQ=255)                      gTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATTCTTATGTCGCAGGGTTTGGTCTCGGCGAAGCAAAGCATCCGCCCCCACT‑GATG‑A‑CCGCA‑GGGGGTTATGGGAAgggggta                                                                                                                                                                                                                          >  1:551336/1‑233 (MQ=255)                                 tgCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGTCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GATG‑A‑CCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAAt                                                                                                                                                                                                            >  1:918817/1‑239 (MQ=255)                                            gtttttCCGCTGCAAAATTTTCTGGCGCTTCGCGACCCTACCAGAAAGGCGAAGTATGTCTCTGTCAAAACAATTAATGCCCCGGAACAGACCCCGGGCCAGCGGGGCCGTTTAAGCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATGGCGTCTCGACGAAGAAATGCACCCGCTCACACT‑GATG‑A‑CCGTG‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGt                                                                                                                                                                                                 <  2:551336/237‑1 (MQ=255)                                               ttcccgcgccacaAATTGCTGGCGCTTGCGGAACCTCCCGGCAATGCGAAGTAGTTCGCTTTCGAAAAAATTTATGCACCGGAACAGATGCCAGGCCGGCAGGACCGCTTTATCGGGGGTGGGCTAAAATTTCTTTTTGTCGAAGGTTTGCGTCTCGACGAAGCAATGCACCCGCTCACACT‑GATG‑A‑CCGTG‑GGTGTTTAGGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCg                                                                                                                                                                                              <  2:610713/230‑1 (MQ=255)                                                  ccGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGAAGGACCGCTTTATCGGCGGTGGGGGGCAAAATACTTATTTCGAAGGGGTGTGGCTCCACCGAGCAACGGCCTCGCTCAACCT‑CATGCA‑CGGTG‑GGGGTTT‑TGGCAAGGGGTTACCGCCACCAAAAATGgcgccggggagggc                                                                                                                                                                                             >  2:23939/1‑229 (MQ=255)                                                       gCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGGAATATTCTTATGTCGGAGGATTGGGTCTCGGCGAAGCAATGCATCCGCCCACACT‑GATG‑A‑CCGGA‑GGTGTTTATTGTCAGGGGGTAACGCCCC‑‑AAAATTGGGCCGCGGTGGGGCTGCTtg                                                                                                                                                                                      >  1:589058/1‑239 (MQ=255)                                                        cacaAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGGAAAATCCTTTTTTTGCAGGAGTTCGTCTCGGCCCAGCAAAGCAACCGCTCCCACC‑TATG‑ATCCTGG‑GGGGTTTTTGGAAGGGGGTTACC‑CCCC‑‑AAAAAAGCGCGCCCGGGGGACTGATtgt                                                                                                                                                                                     >  2:876796/1‑239 (MQ=255)                                                          aaaatggtgggCGCTTGCAGACCCTCCCAGCAATGCGAAGTGTTCCGTTTTAAAAAAATTTTATGCCCCGGAAAAGATGCCAGGCCAGCAGGACCGTTTTATGGGCGGTGGTCTGAAATCTCCTTATGTCGAAGGTTGGCTTCTCGACGAAGAAATGCCTCCGCTCACACT‑GATG‑A‑CGGGA‑GTTGTTTGGGGCAGGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGc                                                                                                                                                                                   <  2:640036/230‑1 (MQ=255)                                                                  tGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGAAATATCCTTTTTTCGAAGGATTGTGTCTCGCCGAAGCAAAGCATCCGCCCACAAC‑TATG‑A‑CCGTAGGGGGTTTTTGGAAAGG‑GTTATCCCCAC‑‑AAAATGGCGCGCCGGTGCGGCTGGTTGTGTCGTGGgaa                                                                                                                                                                           >  1:276482/1‑239 (MQ=255)                                                                     cgttggCCAAACCTCCCGGCAATGCGAAGTATGTCGCTTCCGAAAAAATTTATGCACCGGAACAGATGCCGGGCCAGCAGGACCGCTTTATCGGGGGTGGGCTGAAATATCCTAATGTCGAAGGATGGCGTCTCGACGAGCCATTGCACCCGCTCACACT‑GATG‑C‑CCGTG‑GGTGTTTAGGGCAGGGCTTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAata                                                                                                                                                                         <  1:630057/233‑1 (MQ=255)                                                                     cgcTTGCCGAACCTACCAGCAATGGGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAAGGCATCCGCTCACACT‑GATG‑A‑CCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAata                                                                                                                                                                         <  1:884065/238‑1 (MQ=255)                                                                                          aTGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGGAAAAATCTTTTGTCGGAGGGTTGTGTCTCGCCGGAGCAAATAATCCGCTCACCCT‑ATTT‑AGCGGTA‑GG‑GTGTTTTGCGAGGAGTTTCCGCCCC‑‑CAAAAAGCCCGCCGGTGCGCCTCGTTGTGGCGGGGCGATATCGCTTTTCAGAGATTaaaact                                                                                                                                                   >  2:1311350/1‑236 (MQ=255)                                                                                           tGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGGTGTAAAATCCTTTTTTCGAAGGATTTGGGCTCGGCGGAAGAATGCGTACCCTCACACT‑CCTG‑A‑ACGTA‑GGGGTTTATGGGCAAGAGGTACCGCCCC‑‑ACAAAGGCGCGCCGGGGTGACGCATTTTTGCGTGGGAATATGGCCTTTAAGAGAGTAAATAGa                                                                                                                                                  >  2:391781/1‑239 (MQ=255)                                                                                                aGTATGCCGTTTTAGAAACATTTTATCCCCCGAAACAGATGCCGGGCCAGCAGGCCCGTTTTACGGGGGGGGGGATAAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCACACT‑GATG‑A‑CCGTG‑GGTGTTTATGGAAAGGCTTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGt                                                                                                                                              <  1:1167738/238‑1 (MQ=255)                                                                                                                             ccGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTTTTTCGAAGGATTTCGTCTCGACGAAGCAATGCATCCGCCCACACT‑GATG‑A‑CCGTA‑GGGGTTTATGGGAAGGCGGTACCCCCCC‑‑AAAAAGGCGCGCCGGGGCGACTGATTTTGCCGTGGGAATATTGCTTTTAAGGGATTAAAACGATAGTCAGTATTAAACTGACCCGCGAGCGTCCGc                                                                                                                 >  1:1143458/1‑238 (MQ=255)                                                                                                                                    agagGCCAGGCCAGCAGGCCCGTTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCACCCGCTCACACT‑GATG‑A‑CCGTA‑GGTGTTTAGGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAaccacc                                                                                                         <  1:1347964/239‑1 (MQ=255)                                                                                                                                       gtgcaggacagAAGGACCGATTTCGCGGGGGGGGGCTGAAATATCTTTATGTGGAAGGATGGCGATGCGAGAAAGAAATGCACCCGCTCA‑ACTGGATG‑A‑CTGTG‑GTTGTTTATGGCAGGGCTTAACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTg                                                                                                       <  2:1149170/235‑1 (MQ=255)                                                                                                                                               aagcggGACCGTTTTATCGGGGGGGGGTTAAAATATCCTTATGTGGAAGGATTGCGTCTCGCAGAAGAAATGCACCCGCTCACACG‑GATG‑A‑CCGTG‑GGTGTTTATGGCAGGGCTTTCCCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGc                                                                                              <  2:516168/238‑1 (MQ=255)                                                                                                                                                                cggcggGGGGCTGAAATATCCTTATGGCGAAGGATTGCGTCTCGACGAAGCAATGCCTCCGCTCACCCG‑GATG‑A‑CGGTA‑GGTGTTTAGGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATAc                                                                              <  2:1143458/238‑1 (MQ=255)                                                                                                                                                                           tGAAATATCTTTATGCCGAAGGATTGCGTCTCGACGAAGCAAGCCACCCGCTCACACT‑GATG‑C‑CCGTA‑GGTGTTTATGGAAGGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGc                                                                   <  1:1350485/238‑1 (MQ=255)                                                                                                                                                                           tGAAATACCTTTATGGGGAAGGTGTGCGTCGCGACGAAGCAATGACTCCGCTCACACT‑GATG‑C‑CGGTG‑GGTGTTTAGGGCAAGGCTTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGc                                                                   <  2:589058/238‑1 (MQ=255)                                                                                                                                                                                 atCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GATG‑A‑CCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGtt                                                            <  1:421143/239‑1 (MQ=255)                                                                                                                                                                                       aTGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GATG‑A‑CCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTc                                                                                                                                             >  1:386309/1‑152 (MQ=255)                                                                                                                                                                                       aTGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GATG‑A‑CCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTg                                                                                                       >  1:165025/1‑190 (MQ=255)                                                                                                                                                                                       aTGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GATG‑A‑CCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCg                                                      >  1:466371/1‑239 (MQ=255)                                                                                                                                                                                       aTGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GATG‑A‑CCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCg                                                      >  1:414718/1‑239 (MQ=255)                                                                                                                                                                                       aTGTCGAAGGATTGCGTCTCGAAGAAGCAATGCATCCGCTCACACT‑GATG‑A‑CCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCg                                                      >  2:468618/1‑239 (MQ=255)                                                                                                                                                                                            gAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GATG‑A‑CCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTCAGCTGACCCGCGAGCGTCCGCCACCCACCTGGAATCTGGCAGCGCCTGACGAATACGGATTTTACGCCAATGTTACTACGTATGt                                                 >  1:1090895/1‑239 (MQ=255)                                                                                                                                                                                             aaGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GATG‑A‑CCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGtt                                                >  2:1138231/1‑239 (MQ=255)                                                                                                                                                                                              aGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GATG‑A‑CCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTg                                               >  2:236384/1‑239 (MQ=255)                                                                                                                                                                                                         tcaccgaACCAATGCATCCGCTCACACT‑GATG‑A‑CCGTG‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCgcgc                                    <  1:658134/235‑1 (MQ=255)                                                                                                                                                                                                                 gCAATGCATCCGCTCACACT‑GATG‑A‑CCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTcaca                            >  1:1205098/1‑239 (MQ=255)                                                                                                                                                                                                                     tGCATCCGCTCACACT‑GATG‑A‑CCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCAGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCt                        >  2:478725/1‑239 (MQ=255)                                                                                                                                                                                                                          ccGCTCACACT‑GATG‑A‑CCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGACTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATAAGGTTTTTACGCCAACGTTAATCCGTATGTTGATCAACCGAGCTGGTCACAGGCGAcctc                   >  2:867746/1‑237 (MQ=255)                                                                                                                                                                                                                               cacaCT‑GATG‑A‑CCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGAt              >  1:320423/1‑239 (MQ=255)                                                                                                                                                                                                                                  acT‑GATG‑A‑CCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAg                                                                                                                                    >  1:1232066/1‑118 (MQ=255)                                                                                                                                                                                                                                    cggatg‑a‑CCGTA‑GGTGTTTAGGGAAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTAt          <  1:391781/237‑1 (MQ=255)                                                                                                                                                                                                                                       atg‑a‑CCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTAACGAACGATTTATTgg       >  2:548441/1‑239 (MQ=255)                                                                                                                                                                                                                                             ccGTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCa   >  2:590627/1‑239 (MQ=255)                                                                                                                                                                                                                                               gTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAg  >  1:519752/1‑238 (MQ=255)                                                                                                                                                                                                                                               gTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAg  >  1:116021/1‑238 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                               TATGCGCTGCGTGGAAGCGTGGTCGATGGTGGTGCCGTGGATTGGTTTTCCGCTGCACAAATTGCTGGCGCTTGCCGAACCTACCAGCAATGCGAAGTATGTCGCTTTCGAAACAATTTATGCACCGGAACAGATGCCAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACT‑GATG‑A‑CCGTA‑GGTGTTTATGGCAAGGCGTTACCGCCAC‑‑AAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAG  >  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 < 27 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑