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

ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCT  >  NC_000913/941569‑941857
                                                                                                                                                 |                                                                                                                                                  
actCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACt                                                                                                                                                <  1:952580/149‑1 (MQ=255)
      ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTAt                                                                                                                                          <  1:1527708/149‑1 (MQ=255)
      ccGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCt                                                                                                                                       <  2:394691/151‑1 (MQ=255)
         cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGATTTCTATCtt                                                                                                                                      >  2:1369816/1‑149 (MQ=255)
         cgcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTTCTATCtt                                                                                                                                      >  1:94424/1‑149 (MQ=255)
          gcAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCTTACTAtttt                                                                                                                                      >  2:354929/1‑148 (MQ=255)
           cAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTaa                                                                                                                                   >  2:1430227/1‑150 (MQ=255)
              aTTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTaaaa                                                                                                                                 >  2:602754/1‑147 (MQ=255)
                 gCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCa                                                                                                                             <  1:340888/150‑1 (MQ=255)
                   ggAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGc                                                                                                                          <  1:1197528/151‑1 (MQ=255)
                   ggAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTAATATCTTTaaaaag                                                                                                                               >  2:743040/1‑142 (MQ=255)
                    gAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGCGTTCTATCTTGAAAAAGGACGc                                                                                                                          >  1:1418827/1‑150 (MQ=255)
                     aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCa                                                                                                                        <  2:825777/151‑1 (MQ=255)
                     aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCa                                                                                                                        <  1:274460/151‑1 (MQ=255)
                     aaGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGACTTACTATCTTTAACAGGGACGcc                                                                                                                         >  2:259058/1‑150 (MQ=255)
                           ttaaCTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACCTTATTTCTTTAACCAGCACGCCCGaa                                                                                                                     >  2:214027/3‑148 (MQ=255)
                           tGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTTCTATCTTTAAAAAGCACGCCCGaaa                                                                                                                    >  2:1249515/1‑149 (MQ=255)
                               cTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTTCTATCTTTAAAAAGCACGCCCGAAAACTaa                                                                                                               >  1:486332/1‑150 (MQ=255)
                                 ataCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAAt                                                                                                              <  1:763570/149‑1 (MQ=255)
                                  taCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTTTTATCTTTTACAAGCACGCCACAAAACTAATgg                                                                                                            >  1:400882/1‑149 (MQ=255)
                                  taCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGGTTTCTATCTTTAAAAAGGACGCCCGAAATCTaaatgc                                                                                                           >  1:1267368/1‑147 (MQ=255)
                                   aCCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGGTTACTTTTTTTAACAACGGCGCCAGAAAACTTAAGcc                                                                                                           >  1:1585909/1‑150 (MQ=255)
                                       aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCa                                                                                                       <  2:951278/150‑1 (MQ=255)
                                       aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCa                                                                                                       <  2:1331352/150‑1 (MQ=255)
                                       aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCAt                                                                                                      <  1:1467394/151‑1 (MQ=255)
                                       aCGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCAt                                                                                                      <  2:889770/151‑1 (MQ=255)
                                         ggcggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATTTTGAACAGGCACGCCCGAAATCTAATGCCCGCCTg                                                                                                     >  2:906047/1‑150 (MQ=255)
                                           cggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTAGGGTGGCTTATTTTCTTGAAACGGCACGCCCGAAATCTAATGcaccca                                                                                                       >  1:1097823/1‑141 (MQ=255)
                                            ggcTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatca                                                                                                 <  1:437795/151‑1 (MQ=255)
                                                gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatcat                                                                                             <  2:1247868/151‑1 (MQ=255)
                                                gggCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGatcatca                                                                                              <  2:1106747/150‑1 (MQ=255)
                                                   cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTTTATCTTCAACAGGCCCGCCCCAAAACTAATGCCCGCCTGatcatcttc                                                                                            >  2:1403128/1‑149 (MQ=255)
                                                   cAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGTCGTAATATATTTGAACAGCCCGGCCGAAATTCAATGCCCGCTCCatcatcatc                                                                                            >  1:1402744/1‑149 (MQ=255)
                                                             cAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCtat                                                                                <  2:529521/151‑1 (MQ=255)
                                                              aGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCtata                                                                               <  1:235328/151‑1 (MQ=255)
                                                               gCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAc                                                                              <  1:259134/151‑1 (MQ=255)
                                                                   gTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGTCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGac                                                                          <  2:1510553/151‑1 (MQ=255)
                                                                     cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGaca                                                                         <  2:104310/150‑1 (MQ=255)
                                                                     cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATTTTGAACAGGGACGCCAGAAATCTAATGCCCGCCTGAATATCATCGATCCGCGCTATACCGacac                                                                       >  2:1082347/1‑151 (MQ=255)
                                                                     cGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAACAGGCACCCCCGAAAACTAATGCCCGCATGATCATCATCGATCCCCGCTATACCCacac                                                                       >  1:277963/1‑151 (MQ=255)
                                                                      ggATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACAcc                                                                      <  2:146865/151‑1 (MQ=255)
                                                                            cGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGc                                                                 <  1:1449724/150‑1 (MQ=255)
                                                                              aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCgg                                                              <  1:1645177/151‑1 (MQ=255)
                                                                              aaaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCCGAAATCTAATGCCCGCCTGATCATCATCGATCCGCGCTATACCCACACCCGGGCCgg                                                              >  2:1495131/1‑151 (MQ=255)
                                                                               aaaCAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGGTCCCCGCTATACCGCCACccgggccggg                                                             >  2:408904/1‑151 (MQ=255)
                                                                                  cAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGgcgc                                                          <  2:34555/151‑1 (MQ=255)
                                                                                    gTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGAATTACTATTTTTAACAAGCACGCCAGAAAACTAATGCCCCCAAGATCATCATCGATCCCCGCCATAACGGCACCGGGGGCGGgcgc                                                          >  2:441453/1‑149 (MQ=255)
                                                                                      aaGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGaa                                                       <  2:758498/150‑1 (MQ=255)
                                                                                            gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTTAAAAGGCACGCCCGAAATCTAATGGCCGCAAGATCATCATCGATCCGCGCTCTACCCACACCCGTGCCGGGCGCGAAGATGAg                                                 >  1:1372460/1‑150 (MQ=255)
                                                                                            gTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGGGACTTTCTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCCATACCGACACCGGTGCCGGGggcg                                                         >  2:640263/1‑142 (MQ=255)
                                                                                              aGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTg                                               <  2:1262806/150‑1 (MQ=255)
                                                                                                 gCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGAt                                            <  2:473461/150‑1 (MQ=255)
                                                                                                   tGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGGGGGGGGGATTACTATCTTTAAAAGGCACGCCAGAAATCTAATGCCCGCCTGATCATCATCGGTCCGGCCCATACCGGCACCCGGGGCCGGGGCGGAGATGAGGGGGTccc                                         >  2:954647/1‑151 (MQ=255)
                                                                                                    gTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCt                                        <  2:1090408/151‑1 (MQ=255)
                                                                                                       tGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTAt                                      <  2:1485039/150‑1 (MQ=255)
                                                                                                        ggTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTc                                    <  2:1374198/151‑1 (MQ=255)
                                                                                                         gTAATAACCCTGGCGAAACGCGAATGAGTGGCGGATGGGGGGATTAACAACTTTAACAGGGCAGCCACAAAAATAATGGGCGCAGGGGCCTCATCGAACCCCGTCTTAACGAACCCCGGGGCGGGGGCGAAGATCAGTgggacc                                          >  1:1077376/1‑140 (MQ=255)
                                                                                                          taataaCCCTGGCTAAACGCGAATGAGTGTCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGGCACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCg                                   <  1:397530/150‑1 (MQ=255)
                                                                                                              aaCCCTGGCGAAACGCGAATGAGTGGCGGAGGGGTGGCTTACTTTATTGAAAACGCCCGCCCGAAAATTAAAGCCCGCCTGAGAATCCACGGTGCGCGCCATATCCACCCCCGGGGCGGGGGCGAAGATGAGCGGAGGCCCATTCGTCCgg                              >  2:1673644/1‑151 (MQ=255)
                                                                                                                cccTGGCGAAACGCGAATGAGTGGCGGTTGGGGGTCCTACTTTCTTTAACAGGCAAGCCAGAAAACTAATGCCCGCCTGATCATCATCGATCCGCGCCATAACAACACCGGGGCCGGGCGCGGAGATGAGTGGGTCCCCATTCGTCCggg                             >  1:946070/1‑150 (MQ=255)
                                                                                                                 ccTGGCGAAACGCGAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGATCCCCCTTCGTCCggg                             >  2:1436675/1‑149 (MQ=255)
                                                                                                                 ccTGGCGAAACGCGAATGAGTGGCGGCGGGGGGATTTATTATCCTCAACAGGCACGCCAGAAATCTAATGCCCGCCCGATTATCATCGCACCGCGCTATACCCCCACCGCGGCCCGGCGCGAAGATGAGTGGGTCCCTATTCGTCCggg                             >  1:8387/1‑149 (MQ=255)
                                                                                                                  cTGGCGAAACGCGAATGAGTGGCGGTGGGGGTTATTACTCTCTTTAACAAGCACGCCAGCAAACTAATGGCCGCCACATCCTCATCCATTCCCGCCATACCGACCACCGCGCCGGGGGCCAAGATTGGGGGATCCCTATTCGTCCgggga                           >  2:1046225/1‑148 (MQ=255)
                                                                                                                     gCGAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGCACGCCAGAAAACTAATGCCCGCATGATCATCATCGATCCGCGCTATACCCACACCGGGGCCGGGGGCGAAGATGAGTGGATCCCTATTCGTCCGGGAAc                          >  2:619932/1‑148 (MQ=255)
                                                                                                                      cGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGAt                      <  1:1354010/151‑1 (MQ=255)
                                                                                                                       gAAACGCGAATGAGTGGCGGTGGGGGGACTTACTATCTTGAAAAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGGGCCGGGCGCGAAGATGAGTGGGTCCCTATTCGTCCGGGGACAGa                       >  2:239582/1‑149 (MQ=255)
                                                                                                                        aaaCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGc                    <  2:1578904/151‑1 (MQ=255)
                                                                                                                         aaCGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGcc                   <  1:826266/151‑1 (MQ=255)
                                                                                                                           cgcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGc                 <  1:1249470/151‑1 (MQ=255)
                                                                                                                            gcgAATGAGTGGCGGTGGGGTGGCCTACTATCTTGAAAAGGCACGCCCGAAAGCTCATGCCCGCGTGATCATAATCGATCCGCGCCATACCGACAACGGTGGCGGGGGCGGAGAGGAGGGGGTCCCCCTTCGTCCGGGGACAGATGCCg                  >  1:59497/1‑149 (MQ=255)
                                                                                                                            gcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCa                <  1:483950/151‑1 (MQ=255)
                                                                                                                            gcgAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCa                <  1:1006136/151‑1 (MQ=255)
                                                                                                                            gcgAATGAGTGGCGGTGGGGGGACTTACTATCTTGAACAAGGACGCCCGAAATATAATGCCCGCATGATCATCATCGATCCGCGCTATAACGGCACCGGTGCCGGGGGCGAAGATGAGTGGGTCCCCATTCGTCCGGGGACAGATGcc                   >  2:94703/1‑148 (MQ=255)
                                                                                                                             cgAATGAGTGGCGGTGGGGTGACTTACTATCTTTAACAGGCACGCCAGAAATCTAATGCCCGCATGATCCTCATCGATCCGCGCTATACCGACACCGGGGCCCGGCGCGGAGATGAGGGGGTCCCCCTTCGTCCCGGGACAGATGCCGc                 >  1:1427784/1‑149 (MQ=255)
                                                                                                                              gAATGAGTGGCGGTGGGGGGGCTTACTATCTTGAACAGGCACGCCCGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATAACCACACCGGTGCCGGGCGCGAAGATGAGTGGGTCCCTATTTGTCCGGGAACAGATTCCGCa                >  2:594166/1‑149 (MQ=255)
                                                                                                                                aTGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTgg            <  2:330669/151‑1 (MQ=255)
                                                                                                                                 tGAGTGGCGGTGGGGTGGATTTCTATCTTGAACAAGCACGCCCGAAAATTAATGCCCGCATGCTCATCCTCCATCCGCGCGATACCCCCACCCGTGCCGGGCGCGACGATGAGGGGGTCCCCATTCGTCCGGGAACAGATGCCGCCCTg             >  1:1064129/1‑149 (MQ=255)
                                                                                                                                      ggCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACg      >  1:999432/1‑151 (MQ=255)
                                                                                                                                          gTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCt  <  1:471608/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ACTCCTCCGCGCAAATTGCGGAAGGTTTGAACTATACCTACGGCGGCTGGGCAGATGGCAACAGCCCGTCGGATATCGAAAACAGTAAGCTGGTAGTGCTGTTTGGTAATAACCCTGGCGAAACGCGAATGAGTGGCGGTGGGGTGACTTACTATCTTGAACAGGCACGCCAGAAATCTAATGCCCGCATGATCATCATCGATCCGCGCTATACCGACACCGGTGCCGGGCGCGAAGATGAGTGGATCCCTATTCGTCCGGGAACAGATGCCGCACTGGTTAACGGTCT  >  NC_000913/941569‑941857

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: