Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,669,714 T→G 8.0% intergenic (+38/+37) prpD → / ← acnB 2‑methylcitrate dehydratase/aconitate hydratase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,669,7140TG 8.0% 275.6 / 14.1 99intergenic (+38/+37)prpD/acnB2‑methylcitrate dehydratase/aconitate hydratase B
Reads supporting (aligned to +/- strand):  ref base T (49/42);  new base G (5/3);  total (54/45)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

GCGACCCGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTC  >  NC_002947/2669566‑2669862
                                                                                                                                                    |                                                                                                                                                    
gCGACCCGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGGACAGGCCCAAACAACAAAGCCCCGGCATc                                                                                                                                                     >  2:1614339/1‑150 (MQ=255)
  gACCCGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCg                                                                                                                                                    <  1:1234371/149‑1 (MQ=255)
  gACCCGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCg                                                                                                                                                    >  2:1234371/1‑149 (MQ=255)
  gACCCGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGc                                                                                                                                                   <  1:2092791/150‑1 (MQ=255)
  gACCCGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGCGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGc                                                                                                                                                   <  1:587306/150‑1 (MQ=255)
     ncGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGc                                                                                                                                                >  2:1197825/2‑150 (MQ=255)
       gCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCg                                                                                                                                              >  2:561597/1‑150 (MQ=255)
           tGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGc                                                                                                                                          <  2:1119152/150‑1 (MQ=255)
             cccGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCtt                                                                                                                                        >  2:615952/1‑150 (MQ=255)
              ccGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTTGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCttt                                                                                                                                       <  1:255994/150‑1 (MQ=255)
              ccGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCttt                                                                                                                                       <  1:740958/150‑1 (MQ=255)
              ccGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCttt                                                                                                                                       <  1:1614339/150‑1 (MQ=255)
              ccGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCAGCGATGCCGGGGCttt                                                                                                                                       <  2:1695622/150‑1 (MQ=255)
                ggCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGt                                                                                                                                     <  1:561597/150‑1 (MQ=255)
                   aGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGTCCGGCACAGGTGCAAACAACAGATCCCGGTCATCGCTGCCTGGGCTTTGTTTa                                                                                                                                  >  2:1288726/1‑150 (MQ=255)
                    gCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTAc                                                                                                                                 <  2:1643777/150‑1 (MQ=255)
                            agcagGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTa                                                                                                                                  <  1:564814/141‑1 (MQ=255)
                            agcagGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTa                                                                                                                                  >  2:564814/1‑141 (MQ=255)
                            agcagGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTa                                                                                                                         >  1:1702946/1‑150 (MQ=255)
                             gcagGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTaa                                                                                                                        >  1:432714/1‑150 (MQ=255)
                              cagGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAAc                                                                                                                       >  1:548891/1‑150 (MQ=255)
                               agGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACt                                                                                                                      <  1:615952/150‑1 (MQ=255)
                                 gTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCCGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCa                                                                                                                    >  1:1278321/1‑150 (MQ=255)
                                  ttttCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAg                                                                                                                   >  1:397560/1‑150 (MQ=255)
                                   tttCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGc                                                                                                                  >  2:1936285/1‑150 (MQ=255)
                                    ttCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCt                                                                                                                 >  1:2254996/1‑150 (MQ=255)
                                         aCTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGc                                                                                                            <  2:397560/150‑1 (MQ=255)
                                            gtgCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTg                                                                                                         >  2:1344418/1‑150 (MQ=255)
                                              gCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGaa                                                                                                       <  1:585319/150‑1 (MQ=255)
                                               cAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAac                                                                                                      >  1:806633/1‑150 (MQ=255)
                                                taGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGGCTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAacc                                                                                                     >  2:207502/2‑150 (MQ=255)
                                                aaGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTAGTGGCCTTGTCTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAacc                                                                                                     >  1:285531/1‑150 (MQ=255)
                                                 aGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAacca                                                                                                    >  1:1261085/1‑150 (MQ=255)
                                                      cagcagAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCggg                                                                                               >  1:1392730/1‑150 (MQ=255)
                                                      cagcagAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCggg                                                                                               >  2:626231/1‑150 (MQ=255)
                                                        gcagAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGAt                                                                                             <  2:548891/150‑1 (MQ=255)
                                                        gcagAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGAt                                                                                             >  2:436129/1‑150 (MQ=255)
                                                        gcagAAGCTGGAAGGCACGGCTGTGCACAGGTTCGAGGACTTGTTTGTGATATGAAGCGGCCGGCACAGCCCAAAACAACAATGTCCCGGCAGCGCGTCCGGGGCTTTGTTTACAACCTGTAGCTCAGCTTACGCCTGAACCACCGCGAt                                                                                             >  2:266228/1‑150 (MQ=255)
                                                           gAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTaa                                                                                                                        >  2:1318675/1‑120 (MQ=255)
                                                           gAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTaa                                                                                                                        <  1:1318675/120‑1 (MQ=255)
                                                           gAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCtt                                                                                          >  1:1027174/1‑150 (MQ=255)
                                                           gAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCtt                                                                                          <  1:436129/150‑1 (MQ=255)
                                                            aaGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTg                                                                                         <  1:410208/150‑1 (MQ=255)
                                                              gCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGc                                                                                       <  1:1936285/150‑1 (MQ=255)
                                                                tGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGt                                                                                     >  2:1292935/1‑150 (MQ=255)
                                                                       cACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGc                                                                              >  1:2148718/1‑150 (MQ=255)
                                                                                  acaGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGa                                                                   >  1:586042/1‑150 (MQ=255)
                                                                                       ttCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCg                                                              <  2:432714/150‑1 (MQ=255)
                                                                                       ttCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCg                                                              <  1:2244626/150‑1 (MQ=255)
                                                                                        tCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCgg                                                             <  1:1292935/150‑1 (MQ=255)
                                                                                          gTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTg                                                                                         >  2:1846322/1‑120 (MQ=255)
                                                                                          gTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTg                                                                                         <  1:1846322/120‑1 (MQ=255)
                                                                                           tGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGa                                                          >  1:1816794/1‑150 (MQ=255)
                                                                                            ggACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGAt                                                         >  2:1383226/1‑150 (MQ=255)
                                                                                             gACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTNANCCACCGNGATCTTGGCGTTNGCTGCCGCTTCGCGGAACTCGGCGNTc                                                        <  1:1197825/150‑1 (MQ=255)
                                                                                              aCTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCt                                                       >  2:358123/1‑150 (MQ=255)
                                                                                                       tgATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGc                                              >  1:2196980/1‑150 (MQ=255)
                                                                                                        gATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTc                                                   <  2:1431122/144‑1 (MQ=255)
                                                                                                        gATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTc                                                   >  1:1431122/1‑144 (MQ=255)
                                                                                                         atatGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCAGCGATGCCGGGGCtt                                                                                                                                        >  1:457152/1‑58 (MQ=21)
                                                                                                         atatGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCAGCGATGCCGGGGCtt                                                                                                                                        <  2:457152/58‑1 (MQ=21)
                                                                                                           atGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAg                                          <  2:2148718/150‑1 (MQ=255)
                                                                                                             gAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGt                                        >  2:265914/1‑150 (MQ=255)
                                                                                                               aGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGc                                                                                       <  2:1493028/101‑1 (MQ=255)
                                                                                                               aGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGc                                                                                       >  1:1493028/1‑101 (MQ=255)
                                                                                                                 cggccggcACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCg                                    <  2:2254996/150‑1 (MQ=255)
                                                                                                                   gccggcACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCttt                                                                                                                                       >  1:622276/1‑49 (MQ=255)
                                                                                                                   gccggcACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCttt                                                                                                                                       <  2:622276/49‑1 (MQ=255)
                                                                                                                   gccggcACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGt                                  <  2:1702946/150‑1 (MQ=255)
                                                                                                                       gcACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAAc                              <  2:1282873/150‑1 (MQ=255)
                                                                                                                       gcACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAAc                              <  1:1344418/150‑1 (MQ=255)
                                                                                                                              cccAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCa                       >  2:72655/1‑150 (MQ=255)
                                                                                                                              cccAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCa                       >  1:89983/1‑150 (MQ=255)
                                                                                                                                cAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGc                     >  1:1828329/1‑150 (MQ=255)
                                                                                                                                 aaaCAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGc                                                                                                            <  1:1258778/62‑1 (MQ=255)
                                                                                                                                 aaaCAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGc                                                                                                            >  2:1258778/1‑62 (MQ=255)
                                                                                                                                 aaaCAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAAGGTCGGCAGcc                    <  2:806633/150‑1 (MQ=255)
                                                                                                                                     aacaaAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGc                >  1:2011204/1‑150 (MQ=255)
                                                                                                                                     aacaaAGCCCCGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCgg                                                                    <  2:483650/98‑1 (MQ=255)
                                                                                                                                     aacaaAGCCCCGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCgg                                                                    >  1:483650/1‑98 (MQ=255)
                                                                                                                                       caaAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTg              >  2:1934201/1‑150 (MQ=255)
                                                                                                                                         aaGCCTGGGCATCGGTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTTGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTc            <  2:669893/150‑1 (MQ=255)
                                                                                                                                         aaGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGaa                                                                                                       >  2:1011836/1‑59 (MQ=255)
                                                                                                                                         aaGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGaa                                                                                                       <  1:1011836/59‑1 (MQ=255)
                                                                                                                                         aaGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTc            <  2:1392730/150‑1 (MQ=255)
                                                                                                                                         aaGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTc            <  1:626231/150‑1 (MQ=255)
                                                                                                                                           gCCCCGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGa          >  1:1700112/1‑150 (MQ=255)
                                                                                                                                            ccccGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGc                                                            <  2:2106909/99‑1 (MQ=255)
                                                                                                                                            ccccGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGc                                                            >  1:2106909/1‑99 (MQ=255)
                                                                                                                                            ccccGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGAt         <  2:89983/150‑1 (MQ=255)
                                                                                                                                            ccccGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGAt         >  1:19882/1‑150 (MQ=255)
                                                                                                                                             cccGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAg                                               <  2:58222/111‑1 (MQ=255)
                                                                                                                                             cccGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAg                                               >  1:58222/1‑111 (MQ=255)
                                                                                                                                             cccGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATg        >  2:518142/1‑150 (MQ=255)
                                                                                                                                              ccGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGt       >  1:500274/1‑150 (MQ=255)
                                                                                                                                                   aTCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGCTGTCTTTc  >  2:1710690/1‑150 (MQ=255)
                                                                                                                                                   aTCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTc  <  2:586042/150‑1 (MQ=255)
                                                                                                                                                   aTCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTc  <  2:1828329/150‑1 (MQ=255)
                                                                                                                                                   aTCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTc  >  1:1017394/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGACCCGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTC  >  NC_002947/2669566‑2669862

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