Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,591,684 T→G 14.0% L133R (CTG→CGG) ‡ 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,6840TG14.0% 809.4 / 110.0 292L133R (CTG→CGG) ‡mccBmethylcrotonyl‑CoA carboxylase biotin‑containing subunit beta
Reads supporting (aligned to +/- strand):  ref base T (109/142);  new base G (17/24);  total (126/166)
Fisher's exact test for biased strand distribution p-value = 8.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

Only 100 of 302 total aligned reads displayed.

CGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCG  >  NC_002947/4591535‑4591833
                                                                                                                                                     |                                                                                                                                                     
cGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCt                                                                                                                                                       >  1:4234563/1‑150 (MQ=255)
 gAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTg                                                                                                                                                      >  2:5426324/1‑150 (MQ=255)
  aaGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGc                                                                                                                                                     <  1:2375663/150‑1 (MQ=255)
   aGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCa                                                                                                                                                    <  1:7278583/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)
      cGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGaa                                                                                                                                                 <  1:2208271/150‑1 (MQ=255)
          ccAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGc                                                                                                                                             <  2:357990/150‑1 (MQ=255)
           cAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGcc                                                                                                                                            >  2:7718395/1‑150 (MQ=255)
                   ggcgTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGc                                                                                                                                    <  2:6951436/150‑1 (MQ=255)
                   ggcgTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGc                                                                                                                                    <  2:7772109/150‑1 (MQ=255)
                       tGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCatct                                                                                                                                >  1:2643102/1‑150 (MQ=255)
                        gATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCatcta                                                                                                                               >  1:6671419/1‑150 (MQ=255)
                           caccGGCATCGGGCGCGTGTCTGGCCGTGACTGCATGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCatttatct                                                                                                                            <  2:8002453/148‑1 (MQ=255)
                             ccGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtgg                                                                                                                          <  2:5887772/150‑1 (MQ=255)
                             ccGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtgg                                                                                                                          >  2:8654645/1‑150 (MQ=255)
                              cGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggt                                                                                                                         >  2:5080760/1‑150 (MQ=255)
                               ggCATCGGGCGCGTGTCTGGCCGTGACTGCATGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCtggtc                                                                                                                        <  2:217321/150‑2 (MQ=255)
                               ggCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggtg                                                                                                                        <  2:8360699/150‑1 (MQ=255)
                               ggCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggtg                                                                                                                        <  2:7045/150‑1 (MQ=255)
                                gCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggtgg                                                                                                                       >  2:2684912/1‑150 (MQ=255)
                                gCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCtggtgg                                                                                                                       >  2:2823575/1‑150 (MQ=255)
                                  aTCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGAc                                                                                                                     <  1:3932767/150‑1 (MQ=255)
                                       ccGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAAAGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCcgg                                                                                                                <  1:7366484/150‑1 (MQ=255)
                                           gTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCggcggc                                                                                                            <  1:3698405/150‑1 (MQ=255)
                                              gAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGcc                                                                                                         >  2:4687886/1‑150 (MQ=255)
                                               aaGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCa                                                                                                        >  1:865834/1‑150 (MQ=255)
                                                aGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCaa                                                                                                       <  2:8151471/150‑1 (MQ=255)
                                                   cGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACtt                                                                                                    >  2:1101705/1‑150 (MQ=255)
                                                    gTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTg                                                                                                   >  2:3346191/1‑150 (MQ=255)
                                                    cgtgactGCATGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTg                                                                                                   <  1:2783597/144‑1 (MQ=255)
                                                     gtgactGCATGATCGTCTGCAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGc                                                                                                  <  2:3269681/145‑1 (MQ=255)
                                                           gCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGcc                                                                                            >  2:7185182/1‑150 (MQ=255)
                                                            cATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCa                                                                                           <  1:5712008/150‑1 (MQ=255)
                                                             aTGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAg                                                                                          >  2:4302922/1‑150 (MQ=255)
                                                             aTGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAg                                                                                          >  2:1197905/1‑150 (MQ=255)
                                                             aTGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAg                                                                                          >  1:1781496/1‑150 (MQ=255)
                                                                   gtctgcAACGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGag                                                                                    >  2:168842/6‑149 (MQ=255)
                                                                   gTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGaa                                                                                    <  2:4848843/150‑1 (MQ=255)
                                                                     ggCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAgt                                                                                  >  1:6511469/1‑150 (MQ=255)
                                                                     ggCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAgt                                                                                  >  2:2317054/1‑150 (MQ=255)
                                                                         aaCGATGCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGttttt                                                                              <  1:6286844/150‑2 (MQ=255)
                                                                         aaCGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTc                                                                              >  2:7335031/1‑150 (MQ=255)
                                                                         aaCGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTc                                                                              <  1:4009264/150‑1 (MQ=255)
                                                                          acgatgCAACGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGtttttc                                                                             <  1:2521530/150‑3 (MQ=255)
                                                                          acgacgCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTcc                                                                             <  1:1116892/150‑1 (MQ=255)
                                                                          acgacgCCACGGTAAAAGGCGGTTCATACTCCCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTcc                                                                             <  1:4798055/150‑1 (MQ=255)
                                                                            gatgcaaCGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGTTTTtccg                                                                           >  1:3120258/7‑149 (MQ=255)
                                                                            gacgCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTcccc                                                                           >  1:6380672/1‑150 (MQ=255)
                                                                              tgcaaCGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGACTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGTTTTTCCGGa                                                                         <  2:724661/146‑1 (MQ=255)
                                                                               gcaaCGGTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGTTTTTCCGGat                                                                        >  2:2045248/4‑149 (MQ=255)
                                                                                ccACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGAcc                                                                       >  1:4503735/1‑150 (MQ=255)
                                                                                  aCGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACcgc                                                                     <  1:1172152/150‑1 (MQ=255)
                                                                                    ggTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGTTTTTCCGGATCGCGa                                                                   <  2:4160407/150‑1 (MQ=255)
                                                                                    ggTGAAGGGCGGCACCTATTATCCGGTCACGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGTTTTTCCGGATCGCGa                                                                   <  2:3958417/150‑1 (MQ=255)
                                                                                    ggTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGa                                                                   <  2:917657/150‑1 (MQ=255)
                                                                                    ggTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGa                                                                   <  2:7446485/150‑1 (MQ=255)
                                                                                    ggTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGa                                                                   >  2:7947204/1‑150 (MQ=255)
                                                                                    ggTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGa                                                                   <  1:2116838/150‑1 (MQ=255)
                                                                                     gTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAg                                                                  <  1:7182018/150‑1 (MQ=255)
                                                                                            gCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCg                                                           <  2:4519851/150‑1 (MQ=255)
                                                                                            gCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCg                                                           <  1:414593/150‑1 (MQ=255)
                                                                                              ggTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCggg                                                         <  1:1420390/150‑1 (MQ=255)
                                                                                               gTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATGTCGGgc                                                        >  1:9022698/1‑150 (MQ=255)
                                                                                                ttCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGgcg                                                       >  1:8084133/1‑150 (MQ=255)
                                                                                                 tCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGTCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCAGTTCGGgcgc                                                      >  2:8091418/1‑150 (MQ=255)
                                                                                                   gTACTACCCGCTGACGGTGAAGAAGCGCCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCAt                                                    >  2:5954252/1‑150 (MQ=255)
                                                                                                    tactacCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCAtt                                                   >  2:8998656/1‑150 (MQ=255)
                                                                                                    tactacCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCAtt                                                   <  1:6760525/150‑1 (MQ=255)
                                                                                                    nactacCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCAtt                                                   >  1:2253915/2‑150 (MQ=255)
                                                                                                     actacCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCAttt                                                  >  2:6584318/1‑150 (MQ=255)
                                                                                                       tacCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTttc                                                <  1:5426324/150‑1 (MQ=255)
                                                                                                       tacCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTttc                                                <  2:6054046/150‑1 (MQ=255)
                                                                                                          ccGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTttcttc                                             <  1:834094/150‑1 (MQ=255)
                                                                                                            gCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCaa                                           <  1:7042577/150‑1 (MQ=255)
                                                                                                             gtcaCGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGTTTTTCCGGATCGCGATCATTTCGGGCGCATCTTCTACaat                                          >  2:3069851/4‑149 (MQ=255)
                                                                                                             cTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAAc                                          <  2:405859/150‑1 (MQ=255)
                                                                                                               gACGGTGAAGAAGCACCTGCGTGCACAACCCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCa                                        <  2:6742610/150‑1 (MQ=255)
                                                                                                               caCGGTGAAGAAGCACCTGCGTGCGCAGGAAATCGCGGGCGAGAACCGTCTGCCCTGCATTTATCTGGTCGATTCCGGCGGCGCGAACCTGCCCAATCAGGACGAGGTTTTTCCGGATCGCGATCATTTCGGGCGCATCTTCTACAATCa                                        <  1:7037208/149‑1 (MQ=255)
                                                                                                                      aagaagCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAAc                                 <  1:1008395/150‑1 (MQ=255)
                                                                                                                      aagaagCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAAc                                 <  2:3648066/150‑1 (MQ=255)
                                                                                                                       agaagCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACa                                >  2:3216817/1‑150 (MQ=255)
                                                                                                                        gaagCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACAt                               >  1:1284904/1‑150 (MQ=255)
                                                                                                                         aagCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATg                              >  2:8257749/1‑150 (MQ=255)
                                                                                                                           gCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAg                            >  2:1494759/1‑150 (MQ=255)
                                                                                                                               cTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAgcgcg                        >  2:7898916/1‑150 (MQ=255)
                                                                                                                               cTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAgcgcg                        >  1:2585599/1‑150 (MQ=255)
                                                                                                                                 gCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCa                      >  1:3038743/1‑150 (MQ=255)
                                                                                                                                    tGCACACACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAgg                   <  1:7104143/150‑1 (MQ=255)
                                                                                                                                         aaaCCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATcc              >  1:5263060/1‑150 (MQ=255)
                                                                                                                                          caCCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATccc             <  1:3359605/149‑1 (MQ=255)
                                                                                                                                             cATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGATCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCa          <  2:6643466/150‑1 (MQ=255)
                                                                                                                                              aTCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAg         >  1:8845464/1‑150 (MQ=255)
                                                                                                                                               tCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGa        <  1:6584318/150‑1 (MQ=255)
                                                                                                                                                 gCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATc      <  2:2213100/150‑1 (MQ=255)
                                                                                                                                                 gCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTAGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATc      >  2:1306976/1‑150 (MQ=255)
                                                                                                                                                     tGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCg  >  2:1986880/1‑150 (MQ=255)
                                                                                                                                                     tGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCg  >  2:8816369/1‑150 (MQ=255)
                                                                                                                                                     tGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCg  >  1:5985423/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGAAGACGTGCCAGCGGCGGGCGTGATTGCCGGCATCGGCCGCGTGGAAGGCGTGGAATGCATGATCGTGGCCAACGACGCCACGGTAAAAGGCGGTTCGTACTACCCGCTGACGGTGAAGAAGCACCTGCGTGCACAAACCATCGCCCTGCAGAACCGCCTGCCGTGCATCTATCTGGTGGACTCCGGCGGCGCCAACTTGCCCCGCCAGGACGAAGTGTTCCCCGACCGCGAGCATTTCGGGCGCATTTTCTTCAACCAGGCCAACATGAGCGCGCAAGGCATCCCGCAGATCGCCG  >  NC_002947/4591535‑4591833

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: