Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,060,344 C→T 9.1% R164C (CGC→TGC)  iciA → chromosome initiation inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,060,3440CT 9.1% 119.5 / 2.3 44R164C (CGC→TGC) iciAchromosome initiation inhibitor
Reads supporting (aligned to +/- strand):  ref base C (22/18);  new base T (2/2);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.69e-01

AACAGAACGTGCTGACCGACCTGGTGGTTGAGGACCAGGAAGTGGGCCTGAAACGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAGCCGCCTGGCCCGTACTCCAGCGATCGTGTACGGCCCGGACGATTTCCTGCAGCATCGCTACCTGGCATCGCTGGGCA  >  NC_002947/1060204‑1060485
                                                                                                                                            |                                                                                                                                             
aacaGAACGTGCTGACCGACCTGGTGGTTGAGGACCAGGAAGTGGGCCTGAAACGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCgg                                                                                                                                       >  6:149097/1‑149 (MQ=255)
aacaGAACGTGCTGACCGACCTGGTGGTTGAGGACCAGGAAGTGGGCCTGAAACGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCgg                                                                                                                                       >  7:10916/1‑149 (MQ=255)
  caGAACGTGCTGACCGACCTGGTGGTTGAGGACCAGGAAGTGGGCCTGAAACGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCggg                                                                                                                                      >  4:307161/1‑148 (MQ=255)
        gTGCTGACCGACCTGGTGGTTGAGGACCAGGAAGTGGGCCTGAAACGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGc                                                                                                                               <  4:274160/149‑1 (MQ=255)
                       gtggtTGAGGACCAGGAAGTGGGCCTGAAACGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCAt                                                                                                                <  7:209600/149‑1 (MQ=255)
                        tggtTGAGGACCAGGAAGTGGGCCTGAAACGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATg                                                                                                               <  3:307161/149‑1 (MQ=255)
                          gtTGAGGACCAGGAAGTGGGCCTGAAACGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGgc                                                                                                             >  3:309524/1‑149 (MQ=255)
                          gtTGAGGACCAGGAAGTGGGCCTGAAACGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGACATTGGCCAGTCCCGTGTTCATGgc                                                                                                             >  7:275157/1‑149 (MQ=255)
                               ggACCAGGAAGTGGGCCTGAAACGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCgg                                                                                                         <  5:149097/148‑1 (MQ=255)
                                   cAGGAAGTGGGCCTGAAACGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACtt                                                                                                    >  7:160564/1‑149 (MQ=255)
                                    aGGAAGTGGGCATGAAACGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTc                                                                                                   <  6:315076/149‑1 (MQ=255)
                                            ggCCTGAAACGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATgc                                                                                                                                               <  7:75757/97‑1 (MQ=255)
                                            ggCCTGAAACGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATgc                                                                                                                                               >  8:75757/1‑97 (MQ=255)
                                            ggCCTGAAACGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGTGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAgg                                                                                           >  4:311602/1‑149 (MQ=255)
                                            ggCCTGAAACGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGTGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAgg                                                                                           >  8:153612/1‑149 (MQ=255)
                                            ggCCTGAAACGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTc                                                                                                                  <  3:313758/126‑1 (MQ=255)
                                            ggCCTGAAACGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTc                                                                                                                  >  4:313758/1‑126 (MQ=255)
                                              ccTGAAACGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGtt                                                                                         <  1:233644/149‑1 (MQ=255)
                                                  aaaCGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGTGGTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGt                                                                                     <  7:153612/149‑1 (MQ=255)
                                                  aaaCGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGTGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGt                                                                                     <  3:311602/149‑1 (MQ=255)
                                                  aaaCGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGt                                                                                     >  3:268478/1‑149 (MQ=255)
                                                   aaCGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTc                                                                                    <  8:160564/149‑1 (MQ=255)
                                                      gCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGcc                                                                                 >  7:152221/1‑149 (MQ=255)
                                                            gTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAgccgc                                                                           >  4:302292/1‑149 (MQ=255)
                                                               ccGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAGCCGCCTg                                                                        >  5:279538/1‑149 (MQ=255)
                                                                  gCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAGCCGCCTGGcc                                                                     >  2:30482/1‑149 (MQ=255)
                                                                    gAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAGCCGCCTGGCCCg                                                                   >  1:68924/1‑149 (MQ=255)
                                                                       gTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCa                                                                                                                             >  8:254795/1‑88 (MQ=255)
                                                                       gTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCa                                                                                                                             <  7:254795/88‑1 (MQ=255)
                                                                        tGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAGCCGCCTGGCCCGTACt                                                               >  2:25829/1‑149 (MQ=255)
                                                                         ggcggcCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAGCCGCCTGGCCCGTACTc                                                              >  7:107036/1‑149 (MQ=255)
                                                                                   ctgTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAGCCGCCTGGCCCGTACTCCAGCGATCgt                                                    <  8:152221/149‑1 (MQ=255)
                                                                                                      cGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAGCCGCCTGGCCCGTACACCAGCGATCGTGTACGGCCCGGACGATTTc                                 <  1:30482/149‑1 (MQ=255)
                                                                                                       ggTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAGCCGCCTGGCCCGTACTCCAGCGATCGTGTACGGCCCGGACGATTTcc                                <  8:10916/149‑1 (MQ=255)
                                                                                                       ggTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAGCCGCCTGGCCCGTACTCCAGCGATCGTGTACGGCCCGGACGATTTcc                                <  4:268478/149‑1 (MQ=255)
                                                                                                        gTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAGCCGCCTGGCCCGTACTCCAGCGATCGTGTACGGCCCGGACGATTTCCt                               >  5:253863/1‑149 (MQ=255)
                                                                                                         tAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAGCCGCCTGGCCCGTACTCCAGCGATCGTGTACGGCCCGGACGATTTCCTg                              <  6:253863/149‑1 (MQ=255)
                                                                                                                 gcacgcaGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAGCCGCCTGGCCCGTACTCCAGCGATCGTGTACGGCCCGGACGATTTCCTgcaattcgg                      >  6:95396/1‑143 (MQ=255)
                                                                                                                     gcaGTTTGCCACTGGGGGCCATGCGCTACCGGACATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAGCCGCCTGGCCCGTACTCCAGCGATCGTGTACGGCCCGGACGATTTCCTGCAGCATCGCTAc                  <  8:275157/149‑1 (MQ=255)
                                                                                                                              cACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAGCCGCCTGGCCCGTACTCCCGCGATCGTGTACGGCCCGGACGATGTCCTGCCGCATCGCTCCCTGGCATCg         >  4:103719/1‑149 (MQ=255)
                                                                                                                                 tGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAGCCGCCTGGCCCGTACTCCAGCGATCGTGTACGGCCCGGACGATTTCCTGCAGCATCGCTACCTGGCATCGCTg      >  1:281759/1‑149 (MQ=255)
                                                                                                                                 tGGGGGCCATGCGCTACCGGGCATGGGCCCGGCCCGTGGTCAGGGCGCGGCACTGCCCCCAAGGTGTTGTCGCCAGCCGCCTGGGCCGGACGGCTGCGCGCGTGTCCGGCCCGGCGGCTGGCCGGCCGCAGCGCGACCTGGCATCGcgg      >  4:249104/1‑147 (MQ=255)
                                                                                                                                   ggggCCATGCGCTACCGTGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAGCCGCCTGGCCCGTACTCCAGCGATCGTGTACGGCCCGGACGATTTGCTGCAGCATCGCTACCTGGCATCGCTgg     <  1:184920/148‑1 (MQ=255)
                                                                                                                                   ggggCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAGCCGCCTGGCCCGTACTCCAGCGATCGTGTACGGCCCGGACGATTTCCTGCAGCATCGCTACCTGGCATCGCTgg     <  8:107036/148‑1 (MQ=255)
                                                                                                                                    gggCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAGCCGCCTGGCCCGTACTCCAGCGATCGTGTACGGCCCGGACGATTTCCTGCAGCATCGCTACCTGGCATCGCTGGGc   <  1:25829/149‑1 (MQ=255)
                                                                                                                                        cATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAGCCGCCTGGCCCGTACTCCAGCGATCGTGTACGGCCCGGACGATTTCCTGCAGCATCGCTACCTGGCATCGCTGGGCa  >  8:258317/1‑146 (MQ=255)
                                                                                                                                            |                                                                                                                                             
AACAGAACGTGCTGACCGACCTGGTGGTTGAGGACCAGGAAGTGGGCCTGAAACGCATGCGTGCCGGCGAAGTGGCGGCCTGCCTGTGTGGGAGTGAGCGGCCGGTAGCCGGCGCACGCAGTTTGCCACTGGGGGCCATGCGCTACCGGGCATTGGCCAGTCCCGTGTTCATGGCGCGGCACTTCCCCCAAGGTTTTGTCGCCAGCCGCCTGGCCCGTACTCCAGCGATCGTGTACGGCCCGGACGATTTCCTGCAGCATCGCTACCTGGCATCGCTGGGCA  >  NC_002947/1060204‑1060485

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: ‑