Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913941,7130TG50.0% 88.8 / 24.6 78V252G (GTG→GGG) dmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (14/25);  new base G (39/0);  total (53/25)
Fisher's exact test for biased strand distribution p-value = 1.77e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCG  >  NC_000913/941565‑941861
                                                                                                                                                    |                                                                                                                                                     
gACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGa                                                                                                                                                     <  1:39156/151‑1 (MQ=255)
    actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtt                                                                                                                                                  >  2:940795/1‑150 (MQ=255)
     ctcctcCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGAcct                                                                                                                                                  >  1:455983/1‑147 (MQ=255)
        ctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGGCTTAc                                                                                                                                                >  1:269587/1‑148 (MQ=255)
           cgcgCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGCCAGATGGAAAAAGCCCGCCGGATATCGAAAAAAGTAAGCTGGAAGGGCTGTTTGGAAATAACCCTGGCGAAACGCGAAGGGGGGGGGGGGGGGGGGATTCCTATCtt                                                                                                                                          >  1:12773/1‑151 (MQ=255)
            gcgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATCttt                                                                                                                                         >  2:186005/1‑150 (MQ=255)
             cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCtt                                                                                                                                          >  1:65068/1‑149 (MQ=255)
               cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAAc                                                                                                                                      <  1:440440/151‑1 (MQ=255)
                aaaTTGCGGAAGGTTTGAACTATGCCTAAGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACa                                                                                                                                     <  2:367149/151‑1 (MQ=255)
                 aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGaa                                                                                                                                       >  1:306782/1‑148 (MQ=255)
                 aaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCtttatatttttaa                                                                                                                                       >  2:60633/1‑137 (MQ=255)
                       ggAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGGTGTGGAAAACAGTAAGCCGGTAGGGCTGTTTGGTCATAACCCTGGCGAACGCCGCATGGGTGGCGGCGGGGtcggctacaatttt                                                                                                                                          >  2:380937/1‑126 (MQ=255)
                       ggAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGc                                                                                                                              >  1:72924/1‑151 (MQ=255)
                        gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGcc                                                                                                                             <  1:522348/151‑1 (MQ=255)
                        gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATCTTGAACAAGCACGcc                                                                                                                             >  2:271310/1‑151 (MQ=255)
                            gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGaa                                                                                                                         <  1:188614/151‑1 (MQ=255)
                             tttGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGGACGCCCGaa                                                                                                                         >  1:940773/1‑150 (MQ=255)
                                 aaCTATTCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTa                                                                                                                    <  2:318316/151‑1 (MQ=255)
                                        ccTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAAAAGGCACGCCCGAAATCTAAAGccc                                                                                                              >  1:105588/1‑150 (MQ=255)
                                         cTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGc                                                                                                            >  1:436629/1‑151 (MQ=255)
                                           aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCAt                                                                                                          <  1:415919/151‑1 (MQ=255)
                                             ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGTTGACTTACTATCTTGAACCGGCGCGCCAAAAATCTAATGCCCGCATGa                                                                                                        <  2:185866/151‑1 (MQ=255)
                                             ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGa                                                                                                        <  2:733580/151‑1 (MQ=255)
                                               cggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatc                                                                                                      <  1:457832/151‑1 (MQ=255)
                                                ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTTGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAAGCACGCCAGAAATCTAATGCCCGCATGatc                                                                                                      >  1:700698/1‑150 (MQ=255)
                                                    gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatca                                                                                                  >  1:576139/1‑150 (MQ=255)
                                                      gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCg                                                                                               <  1:266490/151‑1 (MQ=255)
                                                      gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGCACGCCAGAAATCTAATGCCCGCATGatcat                                                                                                    >  1:1061752/1‑146 (MQ=255)
                                                        aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGa                                                                                              >  1:239255/1‑150 (MQ=255)
                                                         gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTATATCTTTAAAAAGGACGCCCGAAAACTAATGCCCGCATGATTATCATCg                                                                                               >  2:212374/1‑148 (MQ=255)
                                                         gATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATc                                                                                            >  1:663888/1‑151 (MQ=255)
                                                                   gCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGCACGCCCGAAAACTAAAGCCCGCATGATCATCATCGATCCGCGCtat                                                                                   >  2:952274/1‑149 (MQ=255)
                                                                      cGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGATTACTATCTTGAAAAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATAcc                                                                                >  2:792319/1‑149 (MQ=255)
                                                                      cGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGa                                                                              <  2:260225/151‑1 (MQ=255)
                                                                      cGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAAAAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAc                                                                                 >  1:266883/1‑148 (MQ=255)
                                                                         cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGaca                                                                            <  2:345221/150‑1 (MQ=255)
                                                                          ggATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGCACCCCCGAAAACTAAAGCCCGCAAGATCATCATCGATCCGCGCCATACCGACAcc                                                                          >  2:1063595/1‑151 (MQ=255)
                                                                           gATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCg                                                                         <  2:346515/151‑1 (MQ=255)
                                                                            atatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCCATACCGACACCgg                                                                        >  2:828077/1‑151 (MQ=255)
                                                                             tatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAAGCACGCCAGAAATCTAAAGCCCGCATGATCATCATCGATCCCCGCTATACCGACACCgg                                                                        >  2:370270/1‑150 (MQ=255)
                                                                               tCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTg                                                                      >  1:186010/1‑150 (MQ=255)
                                                                                cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGGCCGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACAcccggg                                                                      >  2:723961/1‑145 (MQ=255)
                                                                                cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTGAAAAAGCACGCCAGAAAACTAATGCCCCCATGATCATTCTCGAACCGCGCCATACCGACACCGGGGcc                                                                    >  2:519306/1‑151 (MQ=255)
                                                                                 gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGAGGGGGGAATTACTATCTTTAACAAGCACACCAAAAATCTCATCCCCGCATTATTATCCTCGATCCGCGCTATACCGACACCCGTgac                                                                    >  2:961316/1‑148 (MQ=255)
                                                                                  aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCgg                                                                  <  1:1042347/151‑1 (MQ=255)
                                                                                            gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGAt                                                        <  1:471165/151‑1 (MQ=255)
                                                                                            gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCCCGCTATACCGACACCGGTGCCGGGCGCGa                                                            >  1:440662/1‑147 (MQ=255)
                                                                                               ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGCACGCCAGAAATCTAATGCCCCCATGATCATCATCGATCCGCGCTATACCCGCACCGGTGCCGGGGGCGAAGATg                                                       >  1:203898/1‑149 (MQ=255)
                                                                                                gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGCACGCCAGAAATCTAATGCCCCCATGATCATCATCGATCCGCGCTATACCGACACCCGTGCCGGGggcg                                                             >  1:463132/1‑142 (MQ=255)
                                                                                                     gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGGTc                                               >  2:36793/1‑151 (MQ=255)
                                                                                                        gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGGCTTACTATCTTTAACAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGGCGGGCGCGAAGATGAGGGGATccc                                             >  2:255955/1‑150 (MQ=255)
                                                                                                        gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATccc                                             >  1:266133/1‑150 (MQ=255)
                                                                                                          ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTa                                           <  1:783075/150‑1 (MQ=255)
                                                                                                          ttGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGGCGGGGGCGAAGATGAGTGGATccc                                             >  1:736057/1‑148 (MQ=255)
                                                                                                           tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAt                                          <  2:794/150‑1 (MQ=255)
                                                                                                            ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTc                                        <  1:112851/151‑1 (MQ=255)
                                                                                                                 taaCCCTGGCGAAACGCGAATGAGTGGCGGCGGGGGGGGTTATTATCTTTAACAGGCCGCCCAGAAATCTAATGCCCCCATGATCCTAAACGAACTCAGGTTTTACCACAACGGTGCCGGGCGCCGAGATGGGTggggccccaatcg                                       >  2:109186/1‑136 (MQ=255)
                                                                                                                   aCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCCTCATCGGTCCCCGCTATACCGACACCGGGGCCCGGGGCGGAGATGAGTGGATCCCCATTCGTCCggg                                 >  2:213775/1‑151 (MQ=255)
                                                                                                                     ccTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCggg                                 >  1:411907/1‑149 (MQ=255)
                                                                                                                      cTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATTTTTTACAGGCACGCCAGAAAACTAATGCCCGCCTGGTCATCATCGATCCGCGCCATAACCACACCGGGGCCGGGCGCGAAGATGAGGGGGTCCCCATTTCTCCgggga                               >  2:240198/1‑148 (MQ=255)
                                                                                                                        ggCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGGGCGAAGATGAGTGGATCCCTATTCGTCCGGGGAc                              >  1:724575/1‑149 (MQ=255)
                                                                                                                          cGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGAt                          <  2:48664/151‑1 (MQ=255)
                                                                                                                           gAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGAt                          >  1:939517/1‑150 (MQ=255)
                                                                                                                            aaaCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAAGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGGACAGATg                         >  2:580463/1‑150 (MQ=255)
                                                                                                                              aCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCg                      <  2:780431/151‑1 (MQ=255)
                                                                                                                               cgcgAATGAGTGGCGGTGGGGTGACCTACTATCTTGAACACGCCCGCCCGAAATCTAATGGCCCCGTGATCATTATCCATCCGCGCTATACCCACACCGGTGCCCGGCGCGAAGATGCGGGGATGCCCACTCTCCCggg                                 >  2:80220/1‑139 (MQ=255)
                                                                                                                               cgcgAATGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGGTCATCATCGATCCGCCCTATACCCACACCCGGGCCGGGGGCGGAGAAGAGGGGGTCCCCCTTCGTCCGGGGACAGATGcccc                     >  2:496706/1‑149 (MQ=255)
                                                                                                                                  gAATGAGTGGCGGTGGGGGGACTTAATATCTTGAACAGGCACGCCAGAAATCTAATGCCCCCATGATCATCATCGATCCCCGCTATACCGACACCCGGGCCGGGCGCGAAGATGAGGGGGTCCCCATTTGTCCGGGGACAGATGcc                       >  2:596669/1‑146 (MQ=255)
                                                                                                                                   aaTGAGTGGCGGTGGGGGGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGAAGcc                       >  1:573199/1‑145 (MQ=255)
                                                                                                                                       aGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCTCTGGtt              <  2:1030631/150‑1 (MQ=255)
                                                                                                                                          ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg          >  1:873681/1‑151 (MQ=255)
                                                                                                                                             ggTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTc       >  2:800360/1‑151 (MQ=255)
                                                                                                                                              gTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTc       >  1:1047664/1‑150 (MQ=255)
                                                                                                                                              gTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCt      >  1:558331/1‑151 (MQ=255)
                                                                                                                                               tGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGACGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg     >  1:766180/1‑151 (MQ=255)
                                                                                                                                                ggggTGTCTTTCTATCTTGAACAGGCACGCCACAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCt      <  1:940823/149‑1 (MQ=255)
                                                                                                                                                ggggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTgg    <  2:492390/151‑1 (MQ=255)
                                                                                                                                                ggggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTgg    >  1:927906/1‑151 (MQ=255)
                                                                                                                                                  ggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCg  <  1:662271/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCG  >  NC_000913/941565‑941861

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