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

ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGC  >  NC_000913/941569‑941860                                                                                                                                                  |                                                                                                                                                    actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGTCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTa                                                                                                                                               <  2:38470/151‑1 (MQ=255)    cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGGAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAACGAGTGGC‑GGGGGGGGGACTTa                                                                                                                                               >  2:475081/1‑148 (MQ=255)                  gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGc                                                                                                                                <  1:494845/149‑1 (MQ=255)                  gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCa                                                                                                                               <  1:354664/150‑1 (MQ=255)                   cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGGTTTCTTTCTTTAACa                                                                                                                                   >  2:448775/1‑145 (MQ=255)                     gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGGGGGGGGGGTTACTATCTTTAAAAGGCACGcc                                                                                                                           >  2:208271/1‑151 (MQ=255)                      aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCa                                                                                                                          <  2:589249/151‑1 (MQ=255)                      aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGGGGGGGGGCGTACTATCTTTAACaaggacccc                                                                                                                           >  1:694956/1‑142 (MQ=255)                              aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGAGTTA‑TATCTTGAACAAGGACGCCAGaaa                                                                                                                      >  1:953368/1‑147 (MQ=255)                               aCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTaa                                                                                                                 <  1:688850/151‑1 (MQ=255)                                cTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGGGGGGGGACTTACCATCTTTAACAGGCACGCCCGAAAACTaa                                                                                                                 >  2:998791/1‑150 (MQ=255)                                        aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCAt                                                                                                        <  1:4363/151‑1 (MQ=255)                                        aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGAGGGGGGGCATACTAATTTTAACAAGCCCGCCAAAAAAATAATTccccc                                                                                                          >  1:688942/1‑147 (MQ=255)                                         cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCAt                                                                                                        <  2:574690/150‑1 (MQ=255)                                          ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGGGGGGGGGCTTACTATTTTGAAAAGGCACGCCAGAAAACTAATGCCCGc                                                                                                          >  1:576228/1‑147 (MQ=255)                                            cggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatc                                                                                                    <  2:281724/151‑1 (MQ=255)                                            cggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatc                                                                                                    <  2:711809/151‑1 (MQ=255)                                                tGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTAATATCTTTAAAAAGCACGCCCGAAATATACAGTCCGCATGagcctc                                                                                                 >  2:526660/1‑145 (MQ=255)                                                    cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGa                                                                                            <  2:458528/151‑1 (MQ=255)                                                    cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGa                                                                                            <  1:666030/151‑1 (MQ=255)                                                           cAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCgcgc                                                                                     <  1:176858/151‑1 (MQ=255)                                                             aCAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCta                                                                                   <  2:456503/151‑1 (MQ=255)                                                                   cGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGAGTTACTATCTTGAAAAGGCACGCCCGAAATCTAATGCCCCCATGATAATCATCGATCCGCGCTATACCg                                                                              >  2:983993/1‑150 (MQ=255)                                                                    gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTACTATCTTGAACAAGCACGCCAGAAATTTAATGGCCGCATGATCATCATCGATCCGCGCTATACCGac                                                                            >  2:413485/1‑151 (MQ=255)                                                                      cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTACTATCTTGAAAAGGCACGCCAGAAATCTAATGCCCGCAAGATCATCATCGATCCGCGCCATACCGa                                                                             >  1:920968/1‑148 (MQ=255)                                                                        gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTACTATCTTGAAAAGGCACGCCAGAAATCTAAAGCCCCCATGATCATTATCGATCCGCGCCATACCGACAcc                                                                         >  1:756482/1‑150 (MQ=255)                                                                           atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGAATTTACATTTTTAAAAGGCAAGCCAGAAATCTAAAGCCCGCCTGATCATCATTGATCCGCGCTTTACCGACACCgg                                                                       >  1:654754/1‑149 (MQ=255)                                                                              gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTACTATCTTTAACCAGCACGCCCGAAATCTAATGCCCGCACGATCCTCATCGATCCGCGCCATAACCACACCGGTGcc                                                                   >  1:1003138/1‑150 (MQ=255)                                                                               aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTATTATCTTTAACAGGCACGCCAGCAAAATAATAGCCGCATGATCATTATCGATCCGCGCTATAACCACACCGGGGCCg                                                                  >  2:100300/1‑150 (MQ=255)                                                                                aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTTCTATCTTGAACAGGCACCCCAGAAATCTTATGCCCCCATGATCATCATCCCTCCCCGCTATAACCACACCCGTGCCggg                                                                >  1:856350/1‑151 (MQ=255)                                                                                  aCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgcg                                                              <  2:542499/151‑1 (MQ=255)                                                                                     gTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTTCTATCTTTAACAGGCAAGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCCCGCCATACCGACACCGGGGCCGGGggcg                                                            >  2:532017/1‑150 (MQ=255)                                                                                       aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGaa                                                          <  1:111209/150‑1 (MQ=255)                                                                                           tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTTCTAATTTTAAAAAGCCCGCCAGAAAACTAATGACCGCGTGATCCTCAACGATCCGCGCCGTTCCGAAACCGGTGCCGGGggcgga                                                          >  2:970744/1‑144 (MQ=255)                                                                                            ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTACTTTCTTTAACAAGGACCCCAGAAAACTAATGCCCGCATGATCATCATCGCTCCCCGCCATAACGACACCGGGGCCGGGCGCGGAGATg                                                      >  1:673105/1‑149 (MQ=255)                                                                                             gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGt                                                   <  2:966695/151‑1 (MQ=255)                                                                                              tAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGCTTTATATCTTGAACAGGGCCCCCAGAAAAATAATGCCCGCCTGATCATAATACATCCCCGCTATACCGACACCGGGTCCGGGCGCGAAGagta                                                     >  2:332715/1‑145 (MQ=255)                                                                                               aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTgg                                                 <  1:815618/151‑1 (MQ=255)                                                                                               aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTg                                                  <  2:19458/150‑1 (MQ=255)                                                                                                  gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGAt                                               <  1:1080035/150‑1 (MQ=255)                                                                                                  gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATc                                              <  1:609504/151‑1 (MQ=255)                                                                                                    tGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTTCTATCTTTAAAAAGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCCCGCTATACCCACCCCCGGGCCCGGGGCGAAGATGAGGGGATccc                                            >  1:525689/1‑151 (MQ=255)                                                                                                        tGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAt                                         <  1:385041/150‑1 (MQ=255)                                                                                                         ggTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAtt                                        <  2:1092507/150‑1 (MQ=255)                                                                                                         ggTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTc                                       <  2:549830/151‑1 (MQ=255)                                                                                                         ggTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTAATATTTTTAAAAAGCACGCCAGAAATCTAATGCCCGCCTGAACATCATCGATCCGCGCTATAACGACACCGGGGCCGGGGGCGAAGAAGAGTGGGTCCCTAtt                                        >  1:712080/1‑150 (MQ=255)                                                                                                           taataaCCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTACTATCTTGAACAGGCACGCCCGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGGACCCTAtt                                        >  1:1057907/1‑148 (MQ=255)                                                                                                                 cccTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGa                               <  1:250814/151‑1 (MQ=255)                                                                                                                  ccTGGCGAAACGCGAATGAGTGGC‑GGTGGGGGGGTTTTTTTTCTTTAACAGGCACGCCCGAAAACTAATGCCCGCCTGGTAATCCTCGCTCCGCGCTATACCGACACCGGGGTCGGGCGCGAAGATGAGTGGGACCCCATTCGTCCggg                                >  2:204124/1‑149 (MQ=255)                                                                                                                     ggCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACa                            <  1:1036419/150‑1 (MQ=255)                                                                                                                     ggCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAg                           <  2:211188/151‑1 (MQ=255)                                                                                                                       cGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGAt                         <  1:368811/151‑1 (MQ=255)                                                                                                                       cGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGAt                         <  1:656154/151‑1 (MQ=255)                                                                                                                       cGAAACGCGAATGAGTGGC‑GGTGGGGGGACTTAATATCTTTAAAAGGCCCGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCCCGCTATACCGACACCCGGGCCGGGCGCGAAGAAAAGGGGGACCCTATTCGTCCgggga                              >  2:868081/1‑144 (MQ=255)                                                                                                                        gcaaCGCGAATGAGTGGC‑GGTGGGGGGCCTTACTATCTTTAACAAGCACGCCAGCAATCTAATTCCCGCCTGATCCTCCTCCATCCGCGCCATACCCACCCCCGCGCCGGGGGCCCAGCAGAGGGGACCCCCCTTCGCCCGGGacccgacg                        >  2:328771/3‑144 (MQ=255)                                                                                                                        gAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATg                        <  1:723027/151‑1 (MQ=255)                                                                                                                         aaaCGCGAATGAGTGGC‑GGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGGCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACa                            >  1:897776/1‑146 (MQ=255)                                                                                                                            cgcgAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGc                    <  1:844519/151‑1 (MQ=255)                                                                                                                             gcgAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCa                   <  1:116362/151‑1 (MQ=255)                                                                                                                              cgAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAGCAGATGCCGCAc                  <  2:301993/151‑1 (MQ=255)                                                                                                                               gAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCAGCATCGATCCGCGCTATACCGACGCCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACt                 <  2:675024/151‑1 (MQ=255)                                                                                                                                 aTGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTgg               <  2:538445/151‑1 (MQ=255)                                                                                                                                   gAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGt              <  1:362298/150‑1 (MQ=255)                                                                                                                                       ggC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg         >  2:870971/1‑151 (MQ=255)                                                                                                                                               gggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGc  <  1:800211/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGC‑GGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGC  >  NC_000913/941569‑941860

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