Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,669,717 C→A 15.4% intergenic (+41/+34) 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,7170CA15.4% 98.1 / 8.3 39intergenic (+41/+34)prpD/acnB2‑methylcitrate dehydratase/aconitate hydratase B
Reads supporting (aligned to +/- strand):  ref base C (16/17);  new base A (4/2);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.70e-01

CCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTCGCG  >  NC_002947/2669580‑2669865
                                                                                                                                         |                                                                                                                                                    
ccGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCttt                                                                                                                                          <  1:724759/150‑1 (MQ=255)
        ggTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACaa                                                                                                                                  <  1:332384/150‑1 (MQ=255)
        ggTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACaa                                                                                                                                  <  1:503658/150‑1 (MQ=255)
         gtgtCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCGCAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAAc                                                                                                                                 <  1:702209/150‑1 (MQ=255)
                 atgTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACt                                                                                                                         <  1:672690/148‑1 (MQ=255)
                                        cagcagAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCggg                                                                                                  >  2:653765/1‑150 (MQ=255)
                                              aaGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTg                                                                                            <  1:118288/150‑1 (MQ=255)
                                              aaGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTg                                                                                            <  1:268623/150‑1 (MQ=255)
                                                    gAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTg                                                                                      >  1:549402/1‑150 (MQ=255)
                                                               tgtgCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTc                                                                           >  1:104674/1‑150 (MQ=255)
                                                                 tgCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTcgc                                                                         <  2:518738/150‑1 (MQ=255)
                                                                   cacaGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACg                                                                                                                <  2:446351/109‑1 (MQ=255)
                                                                   cacaGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACg                                                                                                                >  1:446351/1‑109 (MQ=255)
                                                                   cacaGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGc                                                                                                               >  2:69687/1‑110 (MQ=255)
                                                                   cacaGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGc                                                                                                               <  1:69687/110‑1 (MQ=255)
                                                                        gTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTc                                                                  >  2:491088/1‑150 (MQ=255)
                                                                                aCTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCt                                                          >  1:185963/1‑150 (MQ=255)
                                                                                                aaGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGt                                                       >  1:369150/1‑137 (MQ=255)
                                                                                                aaGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGt                                                       <  2:369150/137‑1 (MQ=255)
                                                                                                   cggccggcAATGTACCAAACAACAAAGCCCCGTCATCGCTGCCGGGGTTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGTGATCTTTGCTTTGGCTGCCGCTTCGCGGATCTCGGCAATCTGGTCCAAGCTCAGGTAGCg                                       <  2:248337/150‑1 (MQ=255)
                                                                                                     gccggcACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGAGCAATGCAGTCAACCTGTAACTCAGCGTACGGCTAAACCACCGGGATCTTGGCGTTTGTTACCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGt                                     >  1:511878/1‑150 (MQ=255)
                                                                                                      ccagcACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCgg                                      <  1:213790/145‑1 (MQ=255)
                                                                                                      ccagcACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCgg                                      >  2:213790/4‑148 (MQ=255)
                                                                                                                 gttgtaaacaaAGCCCCGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAg                         >  1:188043/7‑150 (MQ=255)
                                                                                                                          aaaGCCCCGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAAc                                                                    >  1:57154/1‑98 (MQ=255)
                                                                                                                          aaaGCCCCGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAAc                                                                    <  2:57154/98‑1 (MQ=255)
                                                                                                                             gCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTcgcg                                                                        <  1:211411/91‑1 (MQ=255)
                                                                                                                             gCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTcgcg                                                                        >  2:211411/1‑91 (MQ=255)
                                                                                                                             gCCCCGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGa             >  1:728439/1‑150 (MQ=255)
                                                                                                                             gCCCCGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGa             <  2:188043/150‑1 (MQ=255)
                                                                                                                              ccccGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGAt            >  1:299140/1‑150 (MQ=255)
                                                                                                                              ccccGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAACTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGAt            >  1:350265/1‑150 (MQ=255)
                                                                                                                                ccGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGt          <  1:491088/150‑1 (MQ=255)
                                                                                                                                ccGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGt          >  1:330483/1‑150 (MQ=255)
                                                                                                                                 cGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGt                <  1:374885/143‑1 (MQ=255)
                                                                                                                                 cGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGGCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGt                >  2:374885/1‑143 (MQ=255)
                                                                                                                                     gccGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGt                                                                                        >  2:577863/3‑67 (MQ=255)
                                                                                                                                     gccGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGt                                                                                        <  1:577863/65‑1 (MQ=255)
                                                                                                                                        gCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTcgcg  >  1:183462/1‑150 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
CCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTCGCG  >  NC_002947/2669580‑2669865

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