Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 648,290 A→C 11.9% intergenic (+56/+35) acoR → / ← accC acetoin catabolism regulatory protein/acetyl‑CoA carboxylase biotin carboxylase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947648,2900AC11.9% 84.2 / 5.0 35intergenic (+56/+35)acoR/accCacetoin catabolism regulatory protein/acetyl‑CoA carboxylase biotin carboxylase subunit
Reads supporting (aligned to +/- strand):  ref base A (15/15);  new base C (2/2);  total (17/18)
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.92e-01

GGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCC  >  NC_002947/648141‑648435
                                                                                                                                                     |                                                                                                                                                 
ggCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGca                                                                                                                                                   >  2:612108/1‑150 (MQ=255)
 gCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGcac                                                                                                                                                  >  1:184241/1‑150 (MQ=255)
                gCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCCCTGTGCCGGGGttttt                                                                                                                                   <  2:473057/150‑4 (MQ=255)
                         aTTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATcac                                                                                                                          <  2:184241/150‑1 (MQ=255)
                          ttgcgcgcgcgACGCTTCATCGGCGGGTGGATGTGCTCGGCATCGACATGCCGCGGATGGAGGGCTAGGGTGGTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAACCCCCGGCACAGGGCCGGGGTTCTGGTACATcaca                                                                                                                         <  1:736443/147‑1 (MQ=255)
                                gggCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTacac                                                                                                                   <  2:458956/150‑1 (MQ=255)
                                   cGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTAAACGAt                                                                                                                <  2:610138/150‑1 (MQ=255)
                                                               ccGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTgg                                                                                                                                  <  2:663838/104‑1 (MQ=255)
                                                               ccGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTgg                                                                                                                                  >  1:663838/1‑104 (MQ=255)
                                                                  cATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGt                                                                                 >  2:503075/1‑150 (MQ=255)
                                                                   aTCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTa                                                                                >  2:556940/1‑150 (MQ=255)
                                                                            ccGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGtt                                                                       <  1:612108/150‑1 (MQ=255)
                                                                              gcgGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTgg                                                                                                       >  1:46856/1‑116 (MQ=255)
                                                                              gcgGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTgg                                                                                                       <  2:46856/116‑1 (MQ=255)
                                                                                      aGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCt                                                             >  2:511780/1‑150 (MQ=255)
                                                                                         ggTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGTCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTg                                                          <  1:758936/150‑1 (MQ=255)
                                                                                          gTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTg                                                                                                                                   >  2:19201/1‑76 (MQ=255)
                                                                                          gTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTg                                                                                                                                   <  1:19201/76‑1 (MQ=255)
                                                                                           tAGAGCCCTGGGGCCGCTACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCa                                                        <  1:121628/150‑1 (MQ=255)
                                                                                              agCCCTGGGGCCGCGACGCGGCACCAGATTCTGGCAGGCAATAAAAAACCCCGGCTCGGGGCCGGGGTTCTGGTACAGCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGaa                                                     <  1:8173/150‑1 (MQ=255)
                                                                                                   tGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCtt                                                <  1:503075/150‑1 (MQ=255)
                                                                                                                     ccAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTcgc                              >  2:382794/1‑150 (MQ=255)
                                                                                                                             tGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGc                                                   <  2:530459/121‑1 (MQ=255)
                                                                                                                             tGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGc                                                   >  1:530459/1‑121 (MQ=255)
                                                                                                                              ggcaggcaATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGtt                                                                       <  1:144531/100‑1 (MQ=255)
                                                                                                                              ggcaggcaATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGtt                                                                       >  2:144531/1‑100 (MQ=255)
                                                                                                                                 aggcaATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCCCCAGGTCGCGCTGCAGCgggg                  >  2:375905/1‑149 (MQ=255)
                                                                                                                                 aggcaATAAAAAACCCCGGCACAGGGCCAGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTg                                                                                                    >  1:586658/1‑68 (MQ=255)
                                                                                                                                 aggcaATAAAAAACCCCGGCACAGGGCCAGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTg                                                                                                    <  2:586658/68‑1 (MQ=255)
                                                                                                                                   gcaATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGAc                                                                    <  2:199399/98‑1 (MQ=255)
                                                                                                                                   gcaATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGAc                                                                    >  1:199399/1‑98 (MQ=255)
                                                                                                                                     aaTAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGtt              >  2:792294/1‑150 (MQ=255)
                                                                                                                                          agaaCCCCGGCCCTGTGCCGGGGTTTTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAg                                   <  2:440303/122‑1 (MQ=255)
                                                                                                                                          agaaCCCCGGCCCTGTGCCGGGGTTTTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAg                                   >  1:440303/3‑124 (MQ=255)
                                                                                                                                          aaaaCCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCt         >  1:228315/1‑150 (MQ=255)
                                                                                                                                             aCCCCGGCCCTGTGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGa      >  1:186582/1‑150 (MQ=255)
                                                                                                                                                 cGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGcc  >  1:683769/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GGCTGATGCGGCACGGGCTTTGGGCATTTCGCGGGCGACGCTTTATCGGCGAGTGCATGAGCACCGCATCGACATGCCGCGGATGAAGGGGTAGAGCCCTGGGGCCGCGACGCGGCCCCAGATTCTGGCAGGCAATAAAAAACCCCGGCACAGGGCCGGGGTTCTGGTACATCACACTACACGATCACTCCTGGTTGGCCAGTTTGTGCTCGAGGTAGTGAATGTTGACGCCGCCTTTGCAGAAGCCTTCATCACGCACCAGGTCGCGGTGCAGCGGGATGTTGGTCTTGATGCC  >  NC_002947/648141‑648435

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