Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,097,250 G→C 13.3% intergenic (+241/+51) PP_mr66 → / ← acsA‑I P6/acetyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,097,2500GC13.3% 77.9 / 4.2 30intergenic (+241/+51)PP_mr66/acsA‑IP6/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base G (14/12);  new base C (2/2);  total (16/14)
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.49e-01

TGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGG  >  NC_002947/5097103‑5097388
                                                                                                                                                    |                                                                                                                                          
tGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGa                                                                                                                                           <  2:77708/149‑1 (MQ=255)
       aaaCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGGGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTc                                                                                                                                    >  5:341837/1‑149 (MQ=255)
         aCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTcgc                                                                                                                                  <  3:170705/149‑1 (MQ=255)
         aCGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTcgc                                                                                                                                  <  2:80996/149‑1 (MQ=255)
            ggCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCggg                                                                                                                               <  3:137617/149‑1 (MQ=255)
            ggCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTAAGCTGTCCCCTGACCTCACGCATTCGCggg                                                                                                                               <  4:331544/149‑1 (MQ=21)
              ccGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAGAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTAAGCTGTCCCCTGACCTCACGCATTcgc                                                                                                                                  >  1:165733/1‑144 (MQ=18)
              ccGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAGAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTAAGCTGTCCCCTGACCTCACGCATTcgc                                                                                                                                  <  2:165733/144‑1 (MQ=18)
                 gAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCTGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGcttgc                                                                                                                          >  1:47910/1‑145 (MQ=255)
                         tCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCTGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCTTGCCCGcgaa                                                                                                                   <  2:47910/148‑4 (MQ=255)
                                                                          gatTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGcc                                                                                                                         <  3:151909/93‑1 (MQ=255)
                                                                          gatTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGcc                                                                                                                         >  4:151909/1‑93 (MQ=255)
                                                                                            cAAAAGCGGGTTTGCCCGGTAATGCGTGTGGTCAGGGGACAGCTTACCCTGAACTGACGCATTCGCGGGAAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAg                                                <  7:192691/148‑1 (MQ=255)
                                                                                             aaaaGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGt                                              >  2:93210/1‑149 (MQ=255)
                                                                                                       ttGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACcac                                    >  4:345992/1‑149 (MQ=255)
                                                                                                           ccGCGAATGCGTCAGTTCAGGCTAAGCTGTCCCCTGACCTCACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTGCATCGCGTTGTGGGTTTCGAACCGGTGCTGTCCCacacgc                                >  8:146024/1‑147 (MQ=17)
                                                                                                                     gTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGcc                       >  6:125412/1‑148 (MQ=255)
                                                                                                                      tGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGcc                       >  2:328502/1‑147 (MQ=255)
                                                                                                                      ggAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCaa                     <  6:341837/148‑1 (MQ=255)
                                                                                                                          gTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGAAAGCCCGCTCCTGCAGGGGGTATTGACGGGAAGGGTTTACGCCGAAGCCAGGTTCATAGCCTTGTGGGTTTCGATCAGCTGCAGAGCAACACCCGCGTCAGCCAAGGTc                 >  7:246711/1‑149 (MQ=255)
                                                                                                                             aGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGttt                                                                                        <  6:96620/75‑1 (MQ=255)
                                                                                                                             aGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGttt                                                                                        >  5:96620/1‑75 (MQ=255)
                                                                                                                             aGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAgg                   >  5:151974/1‑144 (MQ=255)
                                                                                                                             aGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAgg                   <  6:151974/144‑1 (MQ=255)
                                                                                                                               gggACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGCAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCtaga             >  1:98114/1‑148 (MQ=255)
                                                                                                                               gggACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCgaga             >  1:256913/1‑148 (MQ=255)
                                                                                                                                       ttAGCCTGAACTGACACATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGcc      >  3:49537/1‑147 (MQ=255)
                                                                                                                                           gCCTGAACTGACGCATTCGGGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAgg  <  1:93210/148‑1 (MQ=255)
                                                                                                                                             cTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCga               <  4:358518/133‑1 (MQ=255)
                                                                                                                                             cTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCga               >  3:358518/1‑133 (MQ=255)
                                                                                                                                                    |                                                                                                                                          
TGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCCCAGG  >  NC_002947/5097103‑5097388

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