Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,965,610 G→A 7.6% intergenic (‑233/‑96) mnxG ← / → PP_3491 manganese‑oxidizing multicopper oxidase/enoly‑coenzyme A hydratase/isomerase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,965,6100GA 7.6% 1085.3 / 55.0 354intergenic (‑233/‑96)mnxG/PP_3491manganese‑oxidizing multicopper oxidase/enoly‑coenzyme A hydratase/isomerase family protein
Reads supporting (aligned to +/- strand):  ref base G (132/194);  new base A (12/15);  total (144/210)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 358 total aligned reads displayed.

TTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATGTACTGGCCGAGGTC  >  NC_002947/3965461‑3965759
                                                                                                                                                     |                                                                                                                                                     
ttGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCg                                                                                                                                                       >  1:2984215/1‑150 (MQ=255)
ttGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCg                                                                                                                                                       >  1:995533/1‑150 (MQ=255)
     cGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATtgt                                                                                                                                                  >  2:3334807/1‑150 (MQ=255)
      gATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATtgtg                                                                                                                                                 >  1:2590847/1‑150 (MQ=255)
        tCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCa                                                                                                                                               <  1:7425428/150‑1 (MQ=255)
        tCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATCGACAATTGTGCa                                                                                                                                               <  2:4012457/150‑1 (MQ=255)
        tCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCACGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCa                                                                                                                                               >  2:7000229/1‑150 (MQ=255)
         cAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCTGGGGACGGAATCCCGTGGGGGCAAATGATGCTCCACAAGCGCGCGCTGGCCTCACCCAAGGTAAAAACGCACGGTACACTTGCACAATTGTCGATTGTGCaa                                                                                                                                              <  2:4137230/150‑1 (MQ=255)
         cAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCaa                                                                                                                                              <  1:6446548/150‑1 (MQ=255)
         cAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCaa                                                                                                                                              <  2:5159107/150‑1 (MQ=255)
         cAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATCGACAATTGTGCaa                                                                                                                                              <  1:5107932/150‑1 (MQ=255)
          aaCCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAg                                                                                                                                             >  2:5018567/1‑150 (MQ=255)
          aaCCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAg                                                                                                                                             <  1:4603644/150‑1 (MQ=255)
           aCCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTCGCACAATCGACAATTGTGCAAgg                                                                                                                                            <  2:526403/150‑1 (MQ=255)
            ccGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGc                                                                                                                                           >  2:4346356/1‑150 (MQ=255)
             cGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCa                                                                                                                                          >  2:3768147/1‑150 (MQ=255)
               cAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCaaa                                                                                                                                        >  1:7289392/1‑150 (MQ=255)
                    gATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGc                                                                                                                                   <  2:1296539/150‑1 (MQ=255)
                    gATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATCCGc                                                                                                                                   <  2:1124039/150‑1 (MQ=255)
                         gggAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGcc                                                                                                                              >  1:3139430/1‑150 (MQ=255)
                         gggAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGcc                                                                                                                              <  2:443078/150‑1 (MQ=255)
                           gaggaTCCCCCAGCAGGGTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAg                                                                                                                            <  2:4683245/150‑1 (MQ=255)
                             ggaTCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGa                                                                                                                          <  2:6405548/150‑1 (MQ=255)
                              gaTCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGAt                                                                                                                         >  1:3841179/1‑150 (MQ=255)
                              gaTCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGAt                                                                                                                         >  2:6060411/1‑150 (MQ=255)
                               aTCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATc                                                                                                                        <  1:3334642/150‑1 (MQ=255)
                                 cccccAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCt                                                                                                                      >  1:6941661/1‑150 (MQ=255)
                                 cccccAGCAGGTACCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCt                                                                                                                      >  1:5632689/1‑150 (MQ=255)
                                   cccAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGc                                                                                                                    >  1:3515963/1‑150 (MQ=255)
                                   cccAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGc                                                                                                                    <  1:1033211/150‑1 (MQ=255)
                                        cagGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCa                                                                                                                             <  2:2488950/136‑1 (MQ=255)
                                            gtCCCCCGATTCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCa                                                                                                           >  1:1171359/2‑150 (MQ=255)
                                              cccccGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGc                                                                                                         <  2:5880030/150‑1 (MQ=255)
                                               ccccGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCg                                                                                                        >  1:6270030/1‑150 (MQ=255)
                                                 ccGAATCGGGTGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGaa                                                                                                      >  2:548252/1‑150 (MQ=255)
                                                 ccGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGaa                                                                                                      >  2:1241897/1‑150 (MQ=255)
                                                 ccGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGaa                                                                                                      <  2:3564208/150‑1 (MQ=255)
                                                 ccGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGaa                                                                                                      <  1:3187949/150‑1 (MQ=255)
                                                  cGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAAc                                                                                                     >  1:3957411/1‑150 (MQ=255)
                                                        gggggACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGa                                                                                               <  1:3541196/150‑1 (MQ=255)
                                                        gggggACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGa                                                                                               <  1:4397306/150‑1 (MQ=255)
                                                          gggACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCCAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCCGACt                                                                                             >  2:2329228/1‑150 (MQ=255)
                                                           ggACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACtt                                                                                            >  2:1512567/1‑150 (MQ=255)
                                                             aCGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCa                                                                                          <  2:1424855/150‑1 (MQ=255)
                                                              cGGAATCCCGTGGGGGCAAAGGCTGTTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCAc                                                                                         <  2:4891031/150‑1 (MQ=255)
                                                               ggAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATCGACAATTGTGCAAgg                                                                                                                                            <  1:6170745/98‑1 (MQ=255)
                                                                gAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACtt                                                                                       <  1:2401492/150‑1 (MQ=255)
                                                                 aaTCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTg                                                                                      >  1:6650293/1‑150 (MQ=255)
                                                                 aaTCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTg                                                                                      >  1:5383334/1‑150 (MQ=255)
                                                                   tCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGa                                                                                    <  2:3238012/150‑1 (MQ=255)
                                                                    cccGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGaa                                                                                   <  2:286607/150‑1 (MQ=255)
                                                                        tGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATc                                                                               >  2:2722946/1‑150 (MQ=255)
                                                                        tGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATc                                                                               <  1:5934683/150‑1 (MQ=255)
                                                                        tGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATc                                                                               <  1:2183724/150‑1 (MQ=255)
                                                                         gggggCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCa                                                                              <  1:3666011/150‑1 (MQ=255)
                                                                           gggCAAAGGCTGCTCCTCAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAAt                                                                            <  1:7348641/150‑1 (MQ=255)
                                                                           gggCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACACGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAAt                                                                            <  2:578116/150‑1 (MQ=255)
                                                                           gggCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAAt                                                                            <  1:4745382/150‑1 (MQ=255)
                                                                           gggCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAAt                                                                            <  2:6346963/150‑1 (MQ=255)
                                                                                    ctgctCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCa                                                                   <  1:4087126/150‑1 (MQ=255)
                                                                                     tgctCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAg                                                                  <  1:2152967/150‑1 (MQ=255)
                                                                                      gctCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAgg                                                                 <  1:4171638/150‑1 (MQ=255)
                                                                                       ctCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAggg                                                                >  2:2116535/1‑150 (MQ=255)
                                                                                         ccACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGgag                                                              >  2:6231018/1‑150 (MQ=255)
                                                                                           acaAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGaa                                                            <  2:2984215/150‑1 (MQ=255)
                                                                                           acaAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGaa                                                            >  2:6889318/1‑150 (MQ=255)
                                                                                            caAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAAc                                                           <  2:1177693/150‑1 (MQ=255)
                                                                                            caAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAAc                                                           <  2:5976239/150‑1 (MQ=255)
                                                                                              aGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGc                                                         >  1:6916747/1‑150 (MQ=255)
                                                                                                cgcgcgcTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCaa                                                       <  1:335247/150‑1 (MQ=255)
                                                                                                 gcgcgcTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAAt                                                      <  2:1054966/150‑1 (MQ=255)
                                                                                                     gcTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACt                                                  <  2:6501412/150‑1 (MQ=255)
                                                                                                         gCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTTCGATCCTCCAGCGAACTCCAGACTTCACTAGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGcac                                              >  1:5313450/1‑150 (MQ=255)
                                                                                                          ccTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAAc                                                                                                     >  2:4845152/1‑94 (MQ=255)
                                                                                                           cTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCacac                                            >  1:7187202/1‑150 (MQ=255)
                                                                                                            tCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACg                                           >  1:5577568/1‑150 (MQ=255)
                                                                                                            tCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACg                                           >  2:1807381/1‑150 (MQ=255)
                                                                                                             cacacaAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGc                                          >  1:1795174/1‑150 (MQ=255)
                                                                                                             cacacaAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGc                                          <  1:1292527/150‑1 (MQ=255)
                                                                                                             cacacaAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGc                                          <  2:1850256/150‑1 (MQ=255)
                                                                                                               cacaAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTc                                        <  1:1227389/150‑1 (MQ=255)
                                                                                                               cacaAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTc                                        >  1:4892105/1‑150 (MQ=255)
                                                                                                                  aaGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATa                                     >  2:1812829/1‑150 (MQ=255)
                                                                                                                  aaGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACGTCACTTGGAAAATCAATAACAAAGCAGGGAGACCGCAATGACTGCACACGCTCATa                                     >  2:7244889/1‑150 (MQ=255)
                                                                                                                     gTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATAccc                                  >  2:7410454/1‑150 (MQ=255)
                                                                                                                                cGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACa                       <  1:1592708/150‑1 (MQ=255)
                                                                                                                                 ggTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAg                      >  1:4938966/1‑150 (MQ=255)
                                                                                                                                     cccTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCa                  <  2:5304339/150‑1 (MQ=255)
                                                                                                                                     cccTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCa                  <  1:2728003/150‑1 (MQ=255)
                                                                                                                                     cccTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCa                  <  2:260856/150‑1 (MQ=255)
                                                                                                                                      ccTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATaa                                                                          >  2:1917614/1‑93 (MQ=255)
                                                                                                                                      ccTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATaa                                                                          <  1:1917614/93‑1 (MQ=255)
                                                                                                                                      ccTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCAt                 >  2:1596851/1‑150 (MQ=255)
                                                                                                                                      ccTTGCACAATCGACAATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCAt                 >  2:1197973/1‑150 (MQ=255)
                                                                                                                                      ccTTGCACAATCGACAATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCAt                 >  2:3119732/1‑150 (MQ=255)
                                                                                                                                       cTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATg                <  2:3163321/150‑1 (MQ=255)
                                                                                                                                       cTTGCACAATCGACAATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATg                >  1:5269799/1‑150 (MQ=255)
                                                                                                                                            acaATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATGTACTg           >  1:1662249/1‑150 (MQ=255)
                                                                                                                                             caATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATGTACTgg          <  1:515423/150‑1 (MQ=255)
                                                                                                                                                     gATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATGTACTGGCCGAGGTc  >  2:6627429/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTGTACGATCAACCGCAACGGATTGGGGAGGATCCCCCAGCAGGTTCCCCCGAATCGGGGGACGGAATCCCGTGGGGGCAAAGGCTGCTCCACAAGCGCGCGCTGGCCTCACACAAGGTAAAAACGCACGGTACCCTTGCACAATTGTCGATTGTGCAAGGCAAATTCGCATGCCAGGATCCTGCGATCCTCCAGCGAACTCCAGACTTCACTTGGAAAATCAATAACAAAGCAGGGAGAACGCAATGACTGCACACGCTCATACCCCGGCAGCAACAGACCATGTACTGGCCGAGGTC  >  NC_002947/3965461‑3965759

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