Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913941,7130TG48.5% 89.2 / 21.8 66V252G (GTG→GGG) dmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (6/28);  new base G (32/0);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 4.05e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGG  >  NC_000913/941569‑941859
                                                                                                                                                |                                                                                                                                                     
actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACtt                                                                                                                                                  <  1:469493/150‑1 (MQ=255)
 ctcctcCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGAGGCGGCGGGGGGAAtt                                                                                                                                                  >  2:446042/1‑149 (MQ=255)
  tcctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtt                                                                                                                                                  >  1:326356/1‑148 (MQ=255)
   cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACt                                                                                                                                               <  2:419201/150‑1 (MQ=255)
      ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATc                                                                                                                                            <  1:972208/150‑1 (MQ=255)
      ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCt                                                                                                                                           <  2:400831/151‑1 (MQ=255)
          gcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGaa                                                                                                                                       <  1:288078/151‑1 (MQ=255)
              aTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTACTATCTTGAACa                                                                                                                                     >  2:442720/1‑149 (MQ=255)
                 gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCa                                                                                                                                 <  1:375345/150‑1 (MQ=255)
                 gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCAc                                                                                                                                <  2:160879/151‑1 (MQ=255)
                  cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACg                                                                                                                               <  1:510612/151‑1 (MQ=255)
                  cGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTaa                                                                                                                                       >  2:794534/1‑143 (MQ=255)
                    gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGGACGcc                                                                                                                             >  2:819017/1‑151 (MQ=255)
                     aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGGACGCCa                                                                                                                            <  1:724208/151‑1 (MQ=255)
                     aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTGAACAGGCACGcc                                                                                                                             >  2:12370/1‑150 (MQ=255)
                       ggTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑Ga                                                                                                                         <  1:347601/151‑1 (MQ=255)
                            gAACTATACCTACGGCGGCTGGGCAGACGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTTCTAAATTTTAGAGGAACGCCA‑GAAAta                                                                                                                     >  1:97401/1‑149 (MQ=255)
                                  taCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAGTTACTATATTTAAAAAGCACGCCA‑GAAAACTAATG‑cc                                                                                                             >  1:273534/1‑151 (MQ=255)
                                   aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGCACGCCA‑GAAAACTTAAG‑ccc                                                                                                            >  1:88037/1‑151 (MQ=255)
                                   aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATCTTTAACAGGCACGCCA‑GAAATCCaaagacc                                                                                                             >  1:610126/1‑146 (MQ=255)
                                        cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATg                                                                                                       <  1:86257/151‑1 (MQ=255)
                                        cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGTGTTTCTATCTTTAACAGGCACGCCA‑GAAATCTAATG‑CCCGCCTg                                                                                                       >  1:376609/1‑151 (MQ=255)
                                         ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCGCGGCC‑GAAATCTAATG‑CCCGc                                                                                                          >  1:113030/1‑147 (MQ=255)
                                             gcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTGAACAGGCCCGCCC‑GAAATCTAATG‑CCCCCATGatca                                                                                                   >  1:371319/1‑150 (MQ=255)
                                                gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGatcatca                                                                                                <  1:231870/150‑1 (MQ=255)
                                                gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAACAAGCACCCCC‑GAAAACTAAAG‑ACCGCCTGattatc                                                                                                 >  1:426904/1‑149 (MQ=255)
                                                  gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCg                                                                                             <  2:77899/151‑1 (MQ=255)
                                                   cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAAGCACCCCA‑GAAATCTAATG‑CCCGCATGatcatcatc                                                                                              >  2:124923/1‑149 (MQ=255)
                                                    aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGAt                                                                                           <  2:781030/151‑1 (MQ=255)
                                                    aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACCCCC‑GAAATCTAATG‑CCCGCATGatcctc                                                                                                 >  2:983932/1‑145 (MQ=255)
                                                    aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGa                                                                                            >  2:104493/1‑150 (MQ=255)
                                                       tGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATcc                                                                                         >  2:753490/1‑150 (MQ=255)
                                                        ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCcgc                                                                                       <  2:266530/151‑1 (MQ=255)
                                                               gCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATAc                                                                                <  1:358491/151‑1 (MQ=255)
                                                               gCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACCCCA‑AGAAACTAAAGCCCCCCATATTCA‑CATCGATCCGCGCtat                                                                                  >  2:636129/1‑149 (MQ=255)
                                                                cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATAcc                                                                               <  1:940035/151‑1 (MQ=255)
                                                                cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATAcc                                                                               <  2:81146/151‑1 (MQ=255)
                                                                cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATAcc                                                                               <  1:475955/151‑1 (MQ=255)
                                                                cccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGGACGCCA‑GAAAACTAATG‑GCCGCATGATCATCATCGATCCGCGCtataac                                                                               >  2:243659/1‑149 (MQ=255)
                                                                   gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCG‑a                                                                            >  2:749422/1‑150 (MQ=255)
                                                                    tCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGTGGGTTAATATCTTGAAAAAGGACGCCA‑GAAATATAATG‑CCCCCGAGATAATCATCGATCCGCGCcata                                                                                 >  1:172853/1‑145 (MQ=255)
                                                                     cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCG‑aca                                                                          <  2:911662/150‑1 (MQ=255)
                                                                       gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGTGGGGTTATATATTTAAAAGGGGCGACAGGAAATTTAA‑G‑GCCGCACGATGATCCTCAATGCGCGCTATACCC‑CCAcc                                                                        >  1:536581/1‑150 (MQ=255)
                                                                             gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAATAGGGACGCCA‑GGAATCTAATG‑CCCCCATGATCATCATCGATACGCGCTACACCG‑CCAACGGGGCCg                                                                 >  2:480940/1‑151 (MQ=255)
                                                                               aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCC‑ACACCGGTGCCggg                                                               >  2:296/1‑151 (MQ=255)
                                                                                aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGgc                                                              <  1:98286/151‑1 (MQ=255)
                                                                                aaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAAAAGGCACGCCC‑GAAATCTAATG‑ACCGCATGATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCgg                                                                >  1:726072/1‑149 (MQ=255)
                                                                                  cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAAGCACGCCA‑GAAATCTAATG‑CCCGCCTGATCATCATCGATTCGCGCTATACCC‑CCACCGGGGCCGGgcg                                                             >  1:414401/1‑150 (MQ=255)
                                                                                     tAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCG‑ACACCGGGGCCGGGcgcg                                                           >  2:417625/1‑149 (MQ=255)
                                                                                       aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTACTATCTTGAACAGGCACGCCC‑GAAAACTAATG‑CCCGCCTGATCATCCTCGATCCGCGCCATACCC‑ACACCGGGGCCCGgcgc                                                            >  1:497447/1‑146 (MQ=255)
                                                                                       aGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGcgcg                                                           >  2:647580/1‑147 (MQ=255)
                                                                                            gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCCCCTTATTATCATCCATCCGCGCTATACCC‑ACACCCGAGCCGGGCGCGgaga                                                       >  2:536347/1‑146 (MQ=255)
                                                                                                 gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGAt                                              <  2:282282/150‑1 (MQ=255)
                                                                                                    gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATccc                                           <  2:260601/150‑1 (MQ=255)
                                                                                                    gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCt                                          <  1:965622/151‑1 (MQ=255)
                                                                                                       tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGAAATATTTTCTTGAAAAGGCACGGCA‑CAAATCTTATG‑TCCGCAGGATGATCATCATTTCGCGCTGTAACC‑AAACCGGAGGCGGGGGCGCAGATGGGTGGATACCt                                          >  2:725608/1‑148 (MQ=255)
                                                                                                        ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTc                                      <  1:52634/151‑1 (MQ=255)
                                                                                                            ataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCC‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCCATACCG‑ACCCCCGGGCCGGGGGCGAAGATGAGTGGATCCCTATTTGTcc                                  >  2:445905/1‑151 (MQ=255)
                                                                                                               aCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCggg                               <  2:326436/151‑1 (MQ=255)
                                                                                                                cccTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGa                              <  1:12369/151‑1 (MQ=255)
                                                                                                                 ccTGGCGAAACGCGAATGAGTGGCGGTCGGGGGTCTTTCTATCTTTGACAAGCGCGCCA‑GAAAACTAAAT‑GCCGCCTGATCATCAACCGGCCGCGCTTTACCGCCCACCGG‑GCCGGGCGGGGAGATGAGTGGGTCCCCATTCGTCCggg                               >  1:367481/1‑149 (MQ=255)
                                                                                                                      cGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGCACGCCA‑GAAAACTAATG‑GCCGCATGATCATCATCGATCCGCGCTATACCG‑ACACCCGGGCCGGGCGCGAAGAAGAGGGGATCCCTATTCGTCCgggga                             >  1:35652/1‑144 (MQ=255)
                                                                                                                           cgcgAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGCCCGCCA‑GAAAAATAATG‑GCCGCATGATTATCATCGATTCGCGCTATACCC‑CCCCCGGGGCCGGGCGCGAAGATGGGGGGATCCCTCTTTGTCCGGGGACAGATGcc                     >  1:821892/1‑149 (MQ=255)
                                                                                                                               aaTGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTg               <  2:753401/151‑1 (MQ=255)
                                                                                                                                         ggTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTc     >  2:133799/1‑151 (MQ=255)
                                                                                                                                            ggggTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTgg  >  2:820455/1‑151 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCA‑GAAATCTAATG‑CCCGCATGATCATCATCGATCCGCGCTATACCG‑ACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGG  >  NC_000913/941569‑941859

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