Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,591,634 G→C 11.3% S116S (TCG→TCC) ‡ mccB → methylcrotonyl‑CoA carboxylase biotin‑containing subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,591,6340GC11.3% 958.8 / 88.1 309S116S (TCG→TCC) ‡mccBmethylcrotonyl‑CoA carboxylase biotin‑containing subunit beta
Reads supporting (aligned to +/- strand):  ref base G (127/146);  new base C (14/21);  total (141/168)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 313 total aligned reads displayed.

GGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCAT  >  NC_002947/4591485‑4591783
                                                                                                                                                     |                                                                                                                                                     
 gCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGt                                                                                                                                                      <  2:2588718/150‑1 (MQ=255)
  cTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGta                                                                                                                                                     <  1:2430625/150‑1 (MQ=255)
   tCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGtac                                                                                                                                                    >  2:6760525/1‑150 (MQ=255)
     gCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGtactc                                                                                                                                                  >  2:8868967/1‑149 (MQ=255)
      ccGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGtactac                                                                                                                                                 <  2:1730829/150‑1 (MQ=255)
      ccGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGtactac                                                                                                                                                 >  1:4338937/1‑150 (MQ=255)
      ccGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGtactac                                                                                                                                                 <  1:2208321/150‑1 (MQ=255)
      ccGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGtactac                                                                                                                                                 >  2:7227834/1‑150 (MQ=255)
          tCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGc                                                                                                                                             <  2:6682741/150‑1 (MQ=255)
          tCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGc                                                                                                                                             <  2:6134457/150‑1 (MQ=255)
            cTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTg                                                                                                                                           <  2:2236681/150‑1 (MQ=255)
             tCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGCATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGGTCGTACTACCCGCTGa                                                                                                                                          >  2:3359605/1‑150 (MQ=255)
             tCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGTCGGTTCGTACTACCCGCTGa                                                                                                                                          >  1:9115706/1‑150 (MQ=255)
             tCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGa                                                                                                                                          >  2:9107555/1‑150 (MQ=255)
             tCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGa                                                                                                                                          >  1:6643530/1‑150 (MQ=255)
             tCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGa                                                                                                                                          >  1:6643416/1‑150 (MQ=255)
             tCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGa                                                                                                                                          >  1:4538765/1‑150 (MQ=255)
               gagaTCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACg                                                                                                                                        >  2:2375663/1‑150 (MQ=255)
                   tCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTga                                                                                                                                    <  2:2528104/150‑1 (MQ=255)
                     ggCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTgaag                                                                                                                                  >  2:700571/1‑150 (MQ=255)
                      gCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTgaaga                                                                                                                                 <  2:6951231/150‑1 (MQ=255)
                       ccAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTgaagaa                                                                                                                                <  1:1084148/150‑1 (MQ=255)
                       ccAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTgaagaa                                                                                                                                <  2:3712923/150‑1 (MQ=255)
                       ccAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTgaagaa                                                                                                                                <  1:248632/150‑1 (MQ=255)
                         aGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGc                                                                                                                              <  2:4526612/150‑1 (MQ=255)
                            tGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCAcc                                                                                                                           >  2:834094/1‑150 (MQ=255)
                             ggCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCt                                                                                                                          <  2:1018679/150‑1 (MQ=255)
                                cgccCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCg                                                                                                                       >  1:159571/1‑150 (MQ=255)
                                 gccCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGt                                                                                                                      >  1:7446485/1‑150 (MQ=255)
                                   ccATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGc                                                                                                                    >  1:3629472/1‑150 (MQ=255)
                                      tGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGta                                                                                                                                                     >  1:9022096/1‑114 (MQ=255)
                                           tgtACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCa                                                                                                            <  2:3975740/150‑1 (MQ=255)
                                            gtACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCAt                                                                                                           <  2:4085517/150‑1 (MQ=255)
                                             tACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATc                                                                                                          >  1:8083783/1‑150 (MQ=255)
                                               cGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGc                                                                                                        >  2:3014792/1‑150 (MQ=255)
                                                ggCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGcc                                                                                                       >  1:6018601/1‑150 (MQ=255)
                                                  cGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCt                                                                                                     <  1:7064182/150‑1 (MQ=255)
                                                  cGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCt                                                                                                     >  1:4234563/1‑150 (MQ=255)
                                                   gAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTg                                                                                                    <  2:8443697/150‑1 (MQ=255)
                                                   gAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTg                                                                                                    <  2:2174301/150‑1 (MQ=255)
                                                     aGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCa                                                                                                  <  1:6702878/150‑1 (MQ=255)
                                                     aGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCa                                                                                                  <  1:4996391/150‑1 (MQ=255)
                                                      gACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAg                                                                                                 <  1:7860446/150‑1 (MQ=255)
                                                      gACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAg                                                                                                 >  1:6951436/1‑150 (MQ=255)
                                                            ccAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGc                                                                                           <  2:357990/150‑1 (MQ=255)
                                                             cAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGcc                                                                                                       >  2:3231742/1‑137 (MQ=255)
                                                             cAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGcc                                                                                          >  2:1172152/1‑150 (MQ=255)
                                                                     ggcgTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGc                                                                                  <  2:7772109/150‑1 (MQ=255)
                                                                      gcgTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTACa                                                                                 <  2:730811/150‑1 (MQ=255)
                                                                        gTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCatc                                                                               >  1:3732888/1‑150 (MQ=255)
                                                                         tGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCatct                                                                              >  1:2642931/1‑150 (MQ=255)
                                                                          gATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCatcta                                                                             >  2:7732520/1‑150 (MQ=255)
                                                                             caccGGCATCGGGCGCGTGTCTGGCCGTGACTGCATGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCatttatct                                                                          <  2:8002453/148‑1 (MQ=255)
                                                                              gCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtg                                                                         <  1:1329412/150‑1 (MQ=255)
                                                                               ccGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtgg                                                                        <  2:5887772/150‑1 (MQ=255)
                                                                               ccGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtgg                                                                        <  2:5008009/150‑1 (MQ=255)
                                                                                cGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggt                                                                       <  2:1520680/150‑1 (MQ=255)
                                                                                cGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggt                                                                       <  1:2680207/150‑1 (MQ=255)
                                                                                 ggCATCGGGCGCGTGTCTGGCCGTGACTGCATGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCtggtc                                                                      <  2:4762680/150‑2 (MQ=255)
                                                                                 ggCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggtg                                                                      <  2:7901094/150‑1 (MQ=255)
                                                                                 ggCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggtg                                                                      <  2:7045/150‑1 (MQ=255)
                                                                                 ggCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGACCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggtg                                                                      <  1:7297986/150‑1 (MQ=255)
                                                                                       ggCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTcc                                                                <  2:7377832/150‑1 (MQ=255)
                                                                                                  aGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCaa                                                     <  2:8151366/150‑1 (MQ=255)
                                                                                                   ggCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAAc                                                    <  1:2527217/150‑1 (MQ=255)
                                                                                                      gTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTg                                                 >  2:2858343/1‑150 (MQ=255)
                                                                                                       tGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGc                                                <  1:2257144/150‑1 (MQ=255)
                                                                                                       tGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGc                                                >  1:3066035/1‑150 (MQ=255)
                                                                                                         gAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGccc                                              <  2:1995387/150‑1 (MQ=255)
                                                                                                              cATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCa                                         <  1:6410176/150‑1 (MQ=255)
                                                                                                               aTGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAg                                        <  2:3367113/150‑1 (MQ=255)
                                                                                                               aTGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAg                                        >  2:4302922/1‑150 (MQ=255)
                                                                                                                tGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAgg                                       >  1:5862973/1‑150 (MQ=255)
                                                                                                                   tcgtctgcAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACg                                    >  1:8220108/8‑150 (MQ=255)
                                                                                                                   tCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACg                                    >  2:3644411/1‑150 (MQ=255)
                                                                                                                       ggCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAgt                                >  1:6511469/1‑150 (MQ=255)
                                                                                                                         ccAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAgtgt                              >  2:5365712/1‑150 (MQ=255)
                                                                                                                          cAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGtttt                             >  1:8466535/1‑150 (MQ=255)
                                                                                                                           aaCGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGttttt                            <  1:6286844/150‑2 (MQ=255)
                                                                                                                           aaCGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTc                            >  2:1251850/1‑150 (MQ=255)
                                                                                                                           aaCGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTc                            >  2:5452081/1‑150 (MQ=255)
                                                                                                                           aaCGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTc                            <  1:4009264/150‑1 (MQ=255)
                                                                                                                            acgacgCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTcc                           <  2:3272859/150‑1 (MQ=255)
                                                                                                                                   cACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACcg                    >  1:4072520/1‑150 (MQ=255)
                                                                                                                                   cACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACcg                    >  1:1901854/1‑150 (MQ=255)
                                                                                                                                   aaCGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGACCCTGCCCAATCAGGACGAGGTTTTTCCGGATcg                    >  2:2759988/2‑150 (MQ=255)
                                                                                                                                    aCGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACcgc                   <  1:1172152/150‑1 (MQ=255)
                                                                                                                                      ggTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGTTTTTCCGGATCGCGa                 <  2:4160407/150‑1 (MQ=255)
                                                                                                                                      ggTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGTTTTTCCGGATCGCGa                 <  1:2142356/150‑1 (MQ=255)
                                                                                                                                      ggTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGa                 <  2:7446485/150‑1 (MQ=255)
                                                                                                                                      ggTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGa                 >  2:2138678/1‑150 (MQ=255)
                                                                                                                                      ggTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGa                 <  1:8672011/150‑1 (MQ=255)
                                                                                                                                      ggTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGa                 <  1:5511809/150‑1 (MQ=255)
                                                                                                                                      ggTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGa                 >  2:9131423/1‑150 (MQ=255)
                                                                                                                                       gTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAg                <  1:7182018/150‑1 (MQ=255)
                                                                                                                                              gCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCg         <  1:414593/150‑1 (MQ=255)
                                                                                                                                               cGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCgg        <  2:763170/150‑1 (MQ=255)
                                                                                                                                               cGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCgg        <  1:2825430/150‑1 (MQ=255)
                                                                                                                                                ggTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCggg       >  1:3195192/1‑150 (MQ=255)
                                                                                                                                                     gTACTACCCGCTGACGGTGAAGAAGCGCCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCAt  >  2:5954252/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCTCGCCGTTCCTCGAGATCGGCCAGCTGGCCGCCCATGAGGTGTACGGCGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCAT  >  NC_002947/4591485‑4591783

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: