Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,414,699 T→C intergenic (+27/‑48) ackA → / → pta acetate kinase A and propionate kinase 2/phosphate acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,414,6990TC90.0% 82.4 / ‑0.3 30intergenic (+27/‑48)ackA/ptaacetate kinase A and propionate kinase 2/phosphate acetyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (10/17);  total (11/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCC  >  NC_000913/2414562‑2414807
                                                                                                                                         |                                                                                                            
cGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTg                                                                                                             <  1:848800/139‑1 (MQ=255)
    tCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGgcgc                                                                                                            >  1:262239/1‑136 (MQ=255)
    tCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGgcgc                                                                                                            <  2:262239/136‑1 (MQ=255)
         aaaTCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGt                                                                                                    <  1:84853/139‑1 (MQ=255)
          aaTCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTa                                                                                                   <  1:249007/139‑1 (MQ=255)
                    catcaACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCaaa                                                                                         <  1:199419/139‑1 (MQ=255)
                     atcaACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTgtgaaat                                                                                                       <  2:863199/124‑7 (MQ=255)
                     atcaACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTgtgaaat                                                                                                       >  1:863199/1‑118 (MQ=255)
                      tcaACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATc                                                                                       >  2:759202/1‑139 (MQ=255)
                          caaAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATcggcg                                                                                   <  1:759202/139‑1 (MQ=255)
                                  gTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGaa                                                                           >  1:587434/1‑139 (MQ=255)
                                     cccGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAaga                                                                        >  2:597136/1‑139 (MQ=255)
                                     cccGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAaga                                                                        >  1:679250/1‑139 (MQ=255)
                                     cccGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAaga                                                                        >  2:88025/1‑139 (MQ=255)
                                         tCCTGCGGTGGTTATCCCAACCAACGATGAGCTTGTTATCGCGCAAGACGCTAGCCGCCTGCCTGCCTGAGTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGAt                                                                    <  1:813209/139‑1 (MQ=255)
                                         tCCTGCGGTGGTTATCCCAACCAACGAAGAACTAGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGAt                                                                    >  1:981355/1‑139 (MQ=255)
                                                 tggtTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCgtg                                                            <  1:979523/139‑1 (MQ=255)
                                                 tggtTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCgtg                                                            <  2:679250/139‑1 (MQ=255)
                                                   gtTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTc                                                          <  1:83657/139‑1 (MQ=255)
                                                   gtTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTc                                                          <  2:806121/139‑1 (MQ=255)
                                                      aTCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTAACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCg                                                       <  1:538940/139‑1 (MQ=255)
                                                              caacGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCg                                                       >  2:398779/1‑131 (MQ=255)
                                                              caacGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCg                                                       <  1:398779/131‑1 (MQ=255)
                                                                acGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATg                                             <  2:28760/139‑1 (MQ=255)
                                                                 cGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTc                                                          <  1:696862/125‑1 (MQ=255)
                                                                 cGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTc                                                          >  2:696862/1‑125 (MQ=255)
                                                                              aTCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCgg                               >  2:667560/1‑139 (MQ=255)
                                                                                      aGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGCTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCATGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCg                       <  1:355924/139‑1 (MQ=255)
                                                                                                           cTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGcc  <  1:98950/139‑1 (MQ=255)
                                                                                                           cTGATTTCACACCGCCAGCTCAGCTGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGcc  <  1:88025/139‑1 (MQ=255)
                                                                                                                                  ctGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGta                                                     <  1:98031/65‑1 (MQ=255)
                                                                                                                                  ctGGCGGCGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGta                                                     >  2:98031/1‑65 (MQ=255)
                                                                                                                                         |                                                                                                            
CGTTTCGGCAAATCTGGTTTCATCAACAAAGAAGGTACCCGTCCTGCGGTGGTTATCCCAACCAACGAAGAACTGGTTATCGCGCAAGACGCGAGCCGCCTGACTGCCTGATTTCACACCGCCAGCTCAGCTGGCGGTGCTGTTTTGTAACCCGCCAAATCGGCGGTAACGAAAGAGGATAAACCGTGTCCCGTATTATTATGCTGATCCCTACCGGAACCAGCGTCGGTCTGACCAGCGTCAGCC  >  NC_000913/2414562‑2414807

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