Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,097,240 A→C 13.3% intergenic (+231/+61) 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,2400AC13.3% 75.0 / 3.4 29intergenic (+231/+61)PP_mr66/acsA‑IP6/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base A (14/11);  new base C (2/2);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCC  >  NC_002947/5097096‑5097384
                                                                                                                                                 |                                                                                                                                                
aCTCCGATGCACCTAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGc                                                                                                                                                <  1:311844/147‑1 (MQ=255)
  tCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTg                                                                                                                                             >  4:170705/1‑148 (MQ=255)
       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                                                                                                                              <  2:165733/144‑1 (MQ=18)
                     ccGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAGAGCGGGCTTGCCCGCGAATGCGTCAGTTCAGGCTAAGCTGTCCCCTGACCTCACGCATTcgc                                                                                                                              >  1:165733/1‑144 (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                                                                                    >  5:96620/1‑75 (MQ=255)
                                                                                                                                    aGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGttt                                                                                    <  6:96620/75‑1 (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)
                                                                                                                                                 |                                                                                                                                                
ACTCCGATGCACCGAAACGGGCCGGAATTGTTTCCTGGCGCCCCGAACCTGTGCGGGGCCGGTGACAGAACAGTAACCGATGATTGTCGGGTAACAATGCAAAAGCGGGCTTGCCCGCGAATGCGTGAGGTCAGGGGACAGCTTAGCCTGAACTGACGCATTCGCGGGCAAGCCCGCTCCTGCAGGGGGTATTGCCGGGAAGGGTTTACGCCGAAGCCAGGTTCATCGCCTTGTGGGTATCGATCAGGTGCTGTACCACACCCGGGTCAGCCAAGGTCGAGATATCGCC  >  NC_002947/5097096‑5097384

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: