Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	941,713	0	T	G	41.9%	101.8 / 18.3	62	V252G (GTG→GGG)	dmsA	dimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (5/31);  new base G (26/0);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.62e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941569‑941862                                                                                                                                                 |                                                                                                                                                       actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACt                                                                                                                                                    <  2:354437/149‑1 (MQ=255)  ctcctcCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtt                                                                                                                                                   >  1:358038/1‑149 (MQ=255)     ctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTAt                                                                                                                                              >  1:1030020/1‑151 (MQ=255)       ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATc                                                                                                                                             <  1:332152/150‑1 (MQ=255)        cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCtt                                                                                                                                           <  1:363273/151‑1 (MQ=255)         gcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCtt                                                                                                                                           >  2:678246/1‑150 (MQ=255)           gcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCtt                                                                                                                                           >  1:40834/1‑148 (MQ=255)            cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAc                                                                                                                                       <  1:1184422/151‑1 (MQ=255)              aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTATTATCTTTaa                                                                                                                                        >  2:222289/1‑148 (MQ=255)                   cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACg                                                                                                                                <  2:339269/151‑1 (MQ=255)                               aCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTa‑a                                                                                                                   <  1:1230764/151‑1 (MQ=255)                                 tataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATTTTGAAAAAGCACGCCCGAAATCTa‑a                                                                                                                   >  1:1091894/1‑149 (MQ=255)                                     ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTTCTATCTTTAACAAGCACCCCAGAAAACTA‑ATGccc                                                                                                              >  1:208005/1‑150 (MQ=255)                                         cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATg                                                                                                         <  1:142651/151‑1 (MQ=255)                                             ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTTGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACTCCAGAAATCTA‑ATGCCCGCATGat                                                                                                       <  2:430748/149‑1 (MQ=255)                                              gcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACCCCCGAAAACTA‑ATGCCCCCATGatca                                                                                                     >  1:932958/1‑150 (MQ=255)                                                 gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGatcatca                                                                                                  <  2:48621/150‑1 (MQ=255)                                                 gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGatcat                                                                                                    >  1:673716/1‑148 (MQ=255)                                                  ggCAGATGGCAACAGCCCGTCGGATATCGAAACCAGTAAGCCGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTCCTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGatcatcatc                                                                                                <  2:286173/151‑1 (MQ=255)                                                     aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGAt                                                                                             <  1:1231913/151‑1 (MQ=255)                                                         ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCcg                                                                                          <  2:1067657/150‑1 (MQ=255)                                                           cAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGTCTAATATCCTTAACAAACAAGCCACAAAACTA‑ATTCCCGCATGATCATCATCGGTCcg                                                                                          >  2:1117569/1‑148 (MQ=255)                                                            aaCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAAAAGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCgcgc                                                                                       >  1:995128/1‑150 (MQ=255)                                                             aCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCta                                                                                     <  2:563797/151‑1 (MQ=255)                                                             aCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAAGCACGCCCGAAAACTA‑ATGCCCGCATGATCATCatccatccc                                                                                          >  2:980624/1‑145 (MQ=255)                                                              cAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCtat                                                                                    <  1:921353/151‑1 (MQ=255)                                                                   cGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAAGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATAcc                                                                                 >  1:629563/1‑149 (MQ=255)                                                                      cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCcgc                                                                                         <  1:486978/138‑1 (MQ=255)                                                                      cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGGGTTACTTCCCTGTACAGGCACGCCAGAAAATTA‑ATGGCCCCATGGACATCATCGATcccgcaataggga                                                                                >  2:487181/1‑136 (MQ=255)                                                                         atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCgg                                                                         <  2:186589/151‑1 (MQ=255)                                                                                aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCggg                                                                  >  1:329827/1‑151 (MQ=255)                                                                                 aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAAAAAGGACGCCCGAAATCTT‑ATGCCCGCATGAGCATTATCGATCCGCGCCATACCGACCCCGGGGCCggg                                                                  >  1:800617/1‑150 (MQ=255)                                                                                    aGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTAATATCTTGAACAAGCCCGCCAGAAAATTA‑AAGCCCGCATGATTATCATCGAACCGCGCTTTACCCACACCCGTGCCGGgcg                                                                >  2:30600/1‑149 (MQ=255)                                                                                       aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGaa                                                            <  1:548546/150‑1 (MQ=255)                                                                                         gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGAt                                                         <  2:80704/151‑1 (MQ=255)                                                                                         gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGAt                                                         <  1:146856/151‑1 (MQ=255)                                                                                         gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTGAACAGGCACGCCCGAAAACTA‑ATGACCGCATGATCATCATCGATCCGCGCTATACCGACACCCGTGCCGGGCGCGGAg                                                           >  2:437027/1‑149 (MQ=255)                                                                                          cTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCCACACCCGTGCCGGGcgcgga                                                            >  1:655973/1‑145 (MQ=255)                                                                                            ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGa                                                       <  1:1082865/150‑1 (MQ=255)                                                                                            ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCTTGAACAGGCCCCCCCCAAAACTA‑AAGCCCCCCTGAACATCCTCGAACCGCGCCATACAGACACCCGTGGCGGGGGCGAAGa                                                          >  2:334964/1‑147 (MQ=255)                                                                                             gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGt                                                     <  1:578582/151‑1 (MQ=255)                                                                                               aGTGCTGCTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTTCTAACTTGAACAGGAACGCCCGAAATCTAGTTGCCCCCA‑GGTCATCATCGTCCCTCGCCATACCGACCCCGGCGCCGGGCAAGGAGATGAGttg                                                   >  2:925241/1‑149 (MQ=255)                                                                                                  gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGAt                                                 <  1:384384/150‑1 (MQ=255)                                                                                                     gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCt                                             <  1:1188251/151‑1 (MQ=255)                                                                                                            aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAGTCGt                                       <  2:467343/150‑1 (MQ=255)                                                                                                               aaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCg                                    <  1:588881/150‑1 (MQ=255)                                                                                                                 cccTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGa                                 <  1:956839/151‑1 (MQ=255)                                                                                                                  ccTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAAACTA‑ATGCCCGCATGATCATCATCGGTCCGCGCTATACCGCCCCCGGGGCCCGGCGCGGAGATGAGTGGGTCCCCATTTGTCCgg                                   >  2:127771/1‑148 (MQ=255)                                                                                                                      gCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGGACGCCAAAAAACTA‑AAGGCCGCATTATCATCATCGGTCCGCGCTATACCGAAACCGGTGCCGGGCGCGGAGATGAGTGGGTCCCCATTCGTCCggg                                  >  1:477254/1‑145 (MQ=255)                                                                                                                        gAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGGACAGAt                           >  2:834773/1‑150 (MQ=255)                                                                                                                          aaCGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGCACGCCAGAAAACTA‑ATGCCCGCATGGTCATTATTCATCCGCGCTATACCGACCCCGGGGCCGGGCGCGAAGATGAGTGGGTCCCCATTTGTCCCGGGACAGATg                          >  2:212744/1‑149 (MQ=255)                                                                                                                          aaCGCGAATGAGTGGCGGGGGGGGGACTTACTATCTTTAAAAGCGACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCCCGCTATAACGACACCGGTGCCGGGGGCGGAGATGAGTGGATCCCCATTCGTCCGGGAACAGaag                          >  2:543933/1‑147 (MQ=255)                                                                                                                           aCGCGAATGAGTGGCGGTGGGGGGAGTTACTTACCTGAAAAAGCACGCCAGAAAATTA‑AAGCCCGCAAGATCATCATCGCTCCGCGCTATACCGAAACCGGTGCCGGGGGGGGAGGTGAGTGGGTACCCATTCGTCCGGGGAAAGAttgc                        >  2:1257025/1‑147 (MQ=255)                                                                                                                               gAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCAGAAAACTA‑ATGCCCGCATGATAATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGGAGATGAGTGGATCCCTATTCGTCCGGGGAAAGATGACGCa                     >  2:772968/1‑149 (MQ=255)                                                                                                                                aaTGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTg                  <  2:1223550/151‑1 (MQ=255)                                                                                                                                    aGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTa              <  1:635382/151‑1 (MQ=255)                                                                                                                                    aGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTa              <  1:636492/151‑1 (MQ=255)                                                                                                                                       ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAAc            <  2:1297664/150‑1 (MQ=255)                                                                                                                                         cGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGt         >  1:918490/1‑151 (MQ=255)                                                                                                                                            tGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTTCGGGAACAGATGCCGCACTGGTTAACGGTCt       >  2:446021/1‑150 (MQ=255)                                                                                                                                            tGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg      >  2:658400/1‑151 (MQ=255)                                                                                                                                                gTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGt  <  2:629842/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                       ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTA‑ATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941569‑941862

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