Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,655,070 C→T 8.0% intergenic (+222/‑36) aceA → / → PP_4118 isocitrate lyase/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,655,0700CT 8.0% 137.2 / 3.0 50intergenic (+222/‑36)aceA/PP_4118isocitrate lyase/membrane protein
Reads supporting (aligned to +/- strand):  ref base C (23/23);  new base T (2/2);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

GCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCAAGCCACCTGATGATTCTCGCCCTGACCGCCGCTGCCCTCTATATATATGCACTGGCCAGCGACAACA  >  NC_002947/4654929‑4655217
                                                                                                                                              |                                                                                                                                                   
gctgGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGca                                                                                                                                                >  1:69779/1‑147 (MQ=255)
   gccaaaTTTTAACCCGCGGTGCACGGTTCCCG‑CGCAAACGTATTCGCCGGTAAACCCGCCCCTTCACTGTTCCGACTGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGc                                                                                                                                           <  1:280420/145‑1 (MQ=16)
         tttttACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGTCGGCACAGGCCACAg                                                                                                                                      >  8:13655/1‑148 (MQ=255)
                     gTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGAt                                                                                                                         >  5:266000/1‑149 (MQ=255)
                         aCCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAgg                                                                                                                     >  1:126640/1‑149 (MQ=255)
                              ttCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGAT‑CGGGATGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGTACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTAt                                                                                                                <  2:326938/149‑1 (MQ=255)
                                 cgccgcTACCGTATTCGCGGGTAAACCCGCACCTACACTGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGc                                                                                                              <  4:153747/148‑1 (MQ=255)
                                 cgccgcTACCGTATTCGCGGGTAAACCCGCACCTACACTGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCt                                                                                                             >  1:119512/1‑149 (MQ=255)
                                     gcTACCGTATTCGCGGGTAAACCCGCACCTACACTGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTtgtg                                                                                                         <  2:14490/149‑1 (MQ=255)
                                      cTACCGTATTCGCGGGTAAACCCGCACCTACACTGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTgg                                                                                                        <  7:82336/149‑1 (MQ=255)
                                       tACCGTATTCGCGGGTAAACCCGCACCTACACTGCT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCGGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGc                                                                                                       >  5:178682/1‑149 (MQ=255)
                                       tACCGTATTCGCGGGTAAAACCGCACCTACACTGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGTCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGc                                                                                                       >  6:13096/1‑149 (MQ=255)
                                         ccGTATTCGCGGGTAAACCCGCACCTACACTGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAg                                                                                                     <  2:119512/149‑1 (MQ=255)
                                          cGTATTCGCGGGTAAACCCGCACCTACACTGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTACACCCGCGAAAGGGCCGGAACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTCTGCCAga                                                                                                    >  8:56653/1‑148 (MQ=255)
                                                 gcgGGTAAACCCGCACCTACACTGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACACGCCACAGCAAAACCCCGGATCAGGGATATGCTTGTGGCAGCCcagcag                                                                                             >  2:115318/1‑149 (MQ=255)
                                                        aaCCCGCACCTACACTGAT‑CGGGAGGTGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGa                                                                                       <  7:56653/148‑1 (MQ=255)
                                                           ccGCGCCTACACTGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGc                                                                                   >  5:263465/1‑149 (MQ=255)
                                                                    cacTGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGTCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCt                                                                           <  5:13096/148‑1 (MQ=255)
                                                                    cacTGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATg                                                                                                               <  1:170371/112‑1 (MQ=255)
                                                                    cacTGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATg                                                                                                               >  2:170371/1‑112 (MQ=255)
                                                                     acTGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAg                                                                                         >  5:265157/1‑133 (MQ=255)
                                                                     acTGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAg                                                                                         <  6:265157/133‑1 (MQ=255)
                                                                     acTGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCg                                                                         >  1:174234/1‑149 (MQ=255)
                                                                       tGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTc                                                                       >  5:135615/1‑149 (MQ=255)
                                                                                aggagACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCAAGCCAc                                                               <  2:69779/149‑1 (MQ=255)
                                                                                  gagaCACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCaa                                                                    <  4:131196/142‑1 (MQ=255)
                                                                                  gagaCACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCaa                                                                    >  3:131196/1‑142 (MQ=255)
                                                                                   agaCACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGTCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCAAGCCACCtg                                                            <  7:13655/149‑1 (MQ=255)
                                                                                     acacacATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCGGCAGCAAg                                                                                         >  5:234812/1‑118 (MQ=255)
                                                                                     acacacATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCGGCAGCAAg                                                                                         <  6:234812/118‑1 (MQ=255)
                                                                                         acaTCGGCGATTGGGCGCAGGCCGCTGCCGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCAAGCCACCTGATGAtt                                                      <  6:218893/149‑1 (MQ=255)
                                                                                            tCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCAAGCCACCTGATGATtctc                                                   <  6:266000/149‑1 (MQ=255)
                                                                                              ggCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCAAGCCACCTGATGATTCTCGc                                                 >  3:88078/1‑149 (MQ=255)
                                                                                               gCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCAAGCCACCTGATGATTCTCGcc                                                <  2:126640/149‑1 (MQ=255)
                                                                                                cGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCAAGCCACCTGATGATTCTCGccc                                               <  6:135615/149‑1 (MQ=255)
                                                                                                             gCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCAt                                                                               <  1:236894/104‑1 (MQ=255)
                                                                                                             gCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCAt                                                                               >  2:236894/1‑104 (MQ=255)
                                                                                                             gCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCAAGCCACCTGATGATTCTCGCCCTGACCGCCGCTGc                                  >  2:285594/1‑149 (MQ=255)
                                                                                                                 cTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCAAGCCACCTGATGATtc                                                     <  4:15259/126‑1 (MQ=255)
                                                                                                                 cTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCAAGCCACCTGATGATtc                                                     >  3:15259/1‑126 (MQ=255)
                                                                                                                 cTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCAAGCCACCTGATGATTCTCGCCCTGACCGCCGCTGCCctc                              >  5:313712/1‑149 (MQ=255)
                                                                                                                     aGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCAAGCCACCTGATGATTCTCGCCCTGACCGCCGCTGCCCTCtata                          >  3:105712/1‑149 (MQ=255)
                                                                                                                      ggAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATTCCTCGTCCAAGCCACCTGAGGATTCTCGCCCTGACCGCCGCTGCCCTCtatat                         <  6:178682/149‑1 (MQ=255)
                                                                                                                          cGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGTATCAGGGCTATGCTTGTGGCAGCCCAGCCGCAAGGACCGCCCATGCCTCGTCCAAGCCACCTGATGATTCTCGCCCTGACCGCCGCTGCCCTCTatatatat                     <  5:98505/149‑1 (MQ=255)
                                                                                                                          cGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCAAGCCACCTGATGATTCTCGCCCTGACCGCCGCTGCCCTCTatatatat                     <  6:263465/149‑1 (MQ=255)
                                                                                                                          cGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCAAGCCACCTGATGATTCTCGCCCTGACCGCCGCTGCCCTCTatatatat                     <  2:174234/149‑1 (MQ=255)
                                                                                                                                  cccGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCAAGCCACCTGATGATTCTCGCCCTGACCGCCGCTGCCCTCTATATATATGCACTGGc             >  5:222962/1‑149 (MQ=255)
                                                                                                                                     gcgAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCAAGCCACCTGATGATTCTCGCCCTGACCGCCGCTGCCCTCTATATATATGCACTGGCCAg          >  6:231375/1‑149 (MQ=255)
                                                                                                                                        aaaGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCAAGCCACCTGATGATTCTCGCCCTGACCGCCGCTGCCCTCTATATATATGCACTGGCCAGCGa       <  5:231375/149‑1 (MQ=255)
                                                                                                                                             gCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCAAGCCACCTGATGATTCTCGCCCTGACCGCCGCTGCCCTCTATATATATGCACTGGCCAGCGacaaca  >  2:164159/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GCTGGCCAATTTTTACCAGCCGTTCACCGTTTCCGCCGCTACCGTATTCGCGGGTAAACCCGCACCTACACTGAT‑CGGGAGGAGACACACATCGGCGATTGGGCGCAGGCCGCTGCAGGAACGGGTTCACCCGCGAAAGGGCCGGCACAGGCCACAGCAAAACCCCGGATCAGGGCTATGCTTGTGGCAGCCCAGCAGCAAGGACCGCCCATGCCTCGTCCAAGCCACCTGATGATTCTCGCCCTGACCGCCGCTGCCCTCTATATATATGCACTGGCCAGCGACAACA  >  NC_002947/4654929‑4655217

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