Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913941,7130TG47.2% 98.2 / 39.2 72V252G (GTG→GGG) dmsAdimethyl sulfoxide reductase, anaerobic, subunit A
Reads supporting (aligned to +/- strand):  ref base T (4/34);  new base G (34/0);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.89e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941563‑941862
                                                                                                                                                      |                                                                                                                                                         
gCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGt                                                                                                                                                           <  2:698898/151‑1 (MQ=255)
gCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGt                                                                                                                                                           <  2:274712/151‑1 (MQ=255)
 cGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTg                                                                                                                                                          <  1:735268/151‑1 (MQ=255)
         cctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACt                                                                                                                                                  <  2:1304522/150‑1 (MQ=255)
          ctccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACt                                                                                                                                                  >  1:987606/1‑149 (MQ=255)
            ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATc                                                                                                                                              <  2:1144683/150‑1 (MQ=255)
               cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGa                                                                                                                                         <  1:285743/151‑1 (MQ=255)
               cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtttccatttt                                                                                                                                             >  1:284956/1‑141 (MQ=255)
               cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCtttatatttt                                                                                                                                             >  2:1102068/1‑141 (MQ=255)
               cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACCTTACATCTTgt                                                                                                                                         >  2:964111/1‑150 (MQ=255)
                  aaaTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACa                                                                                                                                      <  1:1468666/151‑1 (MQ=255)
                    aTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGAAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAgg                                                                                                                                    <  1:335534/151‑1 (MQ=255)
                       gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCa                                                                                                                                  <  2:696363/150‑1 (MQ=255)
                          gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGTATTACTATCTTTAAAAAGCAcccc                                                                                                                              >  1:785764/1‑151 (MQ=255)
                          gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTATTTTTAACAGGCACGcc                                                                                                                              >  1:540827/1‑151 (MQ=255)
                              gTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAACAGGCAcccc                                                                                                                              >  1:140477/1‑147 (MQ=255)
                                 tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAAt                                                                                                                        <  1:457130/150‑1 (MQ=255)
                                   aaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTa                                                                                                                     <  2:301514/151‑1 (MQ=255)
                                     cTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTTATAACTTTGAAAAGGACGCCAGAAAACtta                                                                                                                    >  1:1445314/1‑148 (MQ=255)
                                       ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAAt                                                                                                                   <  1:746497/149‑1 (MQ=255)
                                       ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTATATCTTTAAAAAGCAcccccaaaaactta                                                                                                                    >  1:1162449/1‑138 (MQ=255)
                                           cTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGATTACTATTTTTAAAAAGCCCGCCCAAAAATTAATGacc                                                                                                               >  2:604294/1‑149 (MQ=255)
                                            tACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTTCTATCTTTAAAAAGCACGCCCGAAAAATAATGccccc                                                                                                             >  1:1358354/1‑148 (MQ=255)
                                             aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCa                                                                                                            <  1:1148791/150‑1 (MQ=255)
                                             aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCa                                                                                                            <  2:646868/150‑1 (MQ=255)
                                             aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCAt                                                                                                           <  2:143578/151‑1 (MQ=255)
                                             aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCAt                                                                                                           <  2:3447/151‑1 (MQ=255)
                                             aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGGAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTATATTCTTAAAAAGGACGCCCGAAAtattaaaccc                                                                                                               >  2:622022/1‑138 (MQ=255)
                                              cggcggCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGCACCCCCGAAAACTAAATCCCGcct                                                                                                           >  1:1308228/1‑148 (MQ=255)
                                               ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGa                                                                                                         <  2:284376/151‑1 (MQ=255)
                                               ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTCTCTTGAACAAGCACACCAGAAATCTAATGcc                                                                                                                >  2:1144110/1‑144 (MQ=255)
                                                gcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATCTTTAACAAGCACGCCCGAAATCTAAAGCCCGCCTGat                                                                                                       >  2:1296951/1‑151 (MQ=255)
                                                gcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTTACAAGCACGCCCGAAAACTAATGCCCGCATGa                                                                                                         >  2:1457621/1‑150 (MQ=255)
                                                  ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATCTTTAAAAAGCACGCCCGAAAACTAATGCCCGCCTGa                                                                                                         >  2:642139/1‑148 (MQ=255)
                                                    cTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTACTATCTTTAACAGGCCCGCCCGAAAACTAATGCCCGc                                                                                                             >  1:600095/1‑142 (MQ=255)
                                                       ggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcatc                                                                                                <  2:199071/151‑1 (MQ=255)
                                                        gCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCg                                                                                               <  1:1224806/151‑1 (MQ=255)
                                                          aGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAACAGGCACGCCAGAAATCTAATGCCCCCCTGa                                                                                                         >  1:190775/1‑140 (MQ=255)
                                                              ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcg                                                                                          <  2:1415108/150‑1 (MQ=255)
                                                              ggCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCcg                                                                                          <  1:1364331/150‑1 (MQ=255)
                                                                  aCAGCCCGTCGGCTATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACGCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTAATATTTTGAACAAACGCGCCAGAAAACTAATGACCGCACGATCATCATCCATCGgcgc                                                                                       >  2:1429709/1‑149 (MQ=255)
                                                                    aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCtata                                                                                   <  2:750328/151‑1 (MQ=255)
                                                                       ccGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTTTCTTTAAAAAGCACGCCAGAAAACTAAAGTCCCCCCGATCATTATCGGTCCGCGCCATAcc                                                                                 >  2:839147/1‑150 (MQ=255)
                                                                        cGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCCGAAATCTAAAGCCCGCATGATCATCATCGATCCGCGCTATAcc                                                                                 >  1:9926/1‑149 (MQ=255)
                                                                         gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGac                                                                              <  1:876399/151‑1 (MQ=255)
                                                                         gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACCATCTTTAACAAGCACCCCAGAAAACTAATGGCCGCAAAATCAACATTGaatccc                                                                                         >  1:1213006/1‑135 (MQ=255)
                                                                           cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATTTTGGACAGGCACGCCCGAAAACTAAAGCCCGCATGATCATCATCGATCCGCGCCATACCGac                                                                              >  1:1243107/1‑149 (MQ=255)
                                                                               tatCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCCTGAACAGGCCCCCCAGAAAACTAAAAGCCCCACGATCATTATAAGTCCGCGCTATACCCGCAcc                                                                           >  1:91443/1‑148 (MQ=255)
                                                                                atCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTg                                                                       <  1:1054196/151‑1 (MQ=255)
                                                                                   gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCg                                                                    <  2:1399174/151‑1 (MQ=255)
                                                                                   gAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTTTCTTTAACAAGCACGCCCGAAATATAAAGCCCGCATGATCATCATCGCTCCGCCCTCTACCGACACCGGTGcc                                                                     >  2:768321/1‑150 (MQ=255)
                                                                                    aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTATTTTGAAAAGGCCCGCCAAAAATCTAATGCCCGCACGTTCATCATCCCTCCGCGCGATACCCACACCGGGGcccg                                                                   >  2:875951/1‑149 (MQ=255)
                                                                                     aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCggg                                                                  <  1:207278/151‑1 (MQ=255)
                                                                                        cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACGCCAGAAAATTAATGCCCGCATGATTATCATCGATCCGCGCTATACCGACAACGGGGGCggg                                                                  >  1:227301/1‑148 (MQ=255)
                                                                                              gCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGAt                                                         <  1:1125979/151‑1 (MQ=255)
                                                                                                tGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTAAAATCTTTAACAAGCAGGCCCGAAAAATCAATCCCCCGAGAGCTTCATCCATCCGCGCTATAAAGACCCCCGAGCCTGGGGCGGAggag                                                        >  2:1189348/1‑147 (MQ=255)
                                                                                                 ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGa                                                       <  1:657776/150‑1 (MQ=255)
                                                                                                 ggTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAg                                                      <  1:1166696/151‑1 (MQ=255)
                                                                                                          gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATTTTGAACAAGGACGCCCGAAAACTAAAGCCCCGACGATTATCAAAGATCCGCGCTATTACGACACCCGGGGCCGGCGGGGAGATGAGGGGATccc                                              >  1:411335/1‑150 (MQ=255)
                                                                                                             tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAGAGGCACGCCAGAAAACTAAAGTCCCCCTGGTCATCATCGATGCGCGCTATAACCATACCCGTGGCGGGGGCGAAGATGAGTGGAGCCCCAtt                                          >  2:1408930/1‑151 (MQ=255)
                                                                                                              ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAtt                                          <  2:242024/150‑1 (MQ=255)
                                                                                                                taataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAAGCACGCCAGAAATCTAAAAGCCCCATGATCATCATCGATCCGCGCTATACCGACACCGGGGGCGGGCGCGAAGAAGAGTGGGACCCTCtt                                          >  2:320240/1‑148 (MQ=255)
                                                                                                                 aataaCCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAAGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCCGTGGCGGGCGGGAAGATGAGTGGATCCCTATTCGt                                       >  1:729335/1‑150 (MQ=255)
                                                                                                                   taaCCCTGGCGAAACGCGAATGAGTGGCGGGCTGGTGGCCTTATATTTTGAACCAGCACGCCCCCAAATTAAATCCCCCCTGATTCTTCTTCCTTCGCGCCTTACCGACACCGGGTGCGGGCGCGAAGATGAGTGGAACCCTAGTCGt                                       >  1:543867/1‑148 (MQ=255)
                                                                                                                      cccTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGa                                 <  2:785685/151‑1 (MQ=255)
                                                                                                                       ccTGGCGAAACGCGAATGAGTGGCGGTGGTGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGaa                                <  2:361552/151‑1 (MQ=255)
                                                                                                                          ggCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAg                             <  1:175002/151‑1 (MQ=255)
                                                                                                                            cGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGAt                           <  1:722826/151‑1 (MQ=255)
                                                                                                                              aaaCGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAAGCCCCCCCGAAATCTAATGCCCGCAAGCAAAACATCGATCCGCCCCCTACCGCCCCCGGTGCCGGGCGCGAAGATGAGTGGGTCCCCTTTCGTCCCGGGACAGAt                           >  2:154308/1‑149 (MQ=255)
                                                                                                                                 cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGc                      <  1:1409038/151‑1 (MQ=255)
                                                                                                                                           tGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTaa             >  2:1234200/1‑150 (MQ=255)
                                                                                                                                                 tGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg      >  1:858482/1‑151 (MQ=255)
                                                                                                                                                  ggggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTg      <  2:1243008/150‑1 (MQ=255)
                                                                                                                                                    ggTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCg   >  2:1254667/1‑151 (MQ=255)
                                                                                                                                                      tGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGt  >  1:391717/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                         
GCGACTACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCTGGCGT  >  NC_000913/941563‑941862

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