Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,715,258 A→G 11.7% A21A (GCT→GCC) ‡ yfiE ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,715,2580AG11.7% 142.7 / 8.6 51A21A (GCT→GCC) ‡yfiEputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (29/16);  new base G (2/4);  total (31/20)
Fisher's exact test for biased strand distribution p-value = 1.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAA  >  NC_000913/2715114‑2715406
                                                                                                                                                |                                                                                                                                                     
aTAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTc                                                                                                                                                  >  1:267258/1‑149 (MQ=255)
       tgCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAGCTTCGCaaaaa                                                                                                                                           <  2:162129/149‑1 (MQ=255)
       tgCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGACACCGTCGATTGTGTACAGCACAATTTTTGCGAGCTTCTCaaaaa                                                                                                                                           <  4:238296/149‑1 (MQ=255)
                aaCTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCctc                                                                                                                                  >  4:103739/1‑149 (MQ=255)
                   ttttttCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTc                                                                                                                               >  3:625552/1‑149 (MQ=255)
                   ttttttCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTc                                                                                                                               >  3:796310/1‑149 (MQ=255)
                           ttCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACgg                                                                                                                       >  4:164213/1‑149 (MQ=255)
                                 ggTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTaa                                                                                                                 <  1:318592/149‑1 (MQ=255)
                                       gCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGaa                                                                                                                                       >  2:717346/1‑121 (MQ=255)
                                       gCACATCCGTCGGACAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGaa                                                                                                                                       <  1:717346/121‑1 (MQ=255)
                                            tCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATaaaa                                                                                                      <  3:480758/149‑1 (MQ=255)
                                              cGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACg                                                                                                    <  2:267258/149‑1 (MQ=255)
                                                             tcAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCCGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAAc                                                                                     >  2:73447/1‑149 (MQ=255)
                                                               gaaTAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGTACAATTTTTGCGAAGCTCGCaaaaa                                                                                                                                           >  4:693501/2‑93 (MQ=255)
                                                               aaaTAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCaaaaa                                                                                                                                           <  3:693501/93‑1 (MQ=255)
                                                                      tGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCt                                                                            >  1:147253/1‑149 (MQ=255)
                                                                       ggACTGATAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGCCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTa                                                                           <  3:688973/149‑1 (MQ=255)
                                                                           tGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGa                                                                       >  3:488423/1‑149 (MQ=255)
                                                                              aaaCTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGaaaa                                                                    <  4:796310/149‑1 (MQ=255)
                                                                              aaaCTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGaaaa                                                                    <  3:423079/149‑1 (MQ=255)
                                                                                 cTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAatat                                                                 <  3:164213/149‑1 (MQ=255)
                                                                                 cTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAatat                                                                 <  2:444256/149‑1 (MQ=255)
                                                                                        tCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCt                                                          >  4:303645/1‑149 (MQ=255)
                                                                                         cAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTa                                                         <  3:303645/149‑1 (MQ=255)
                                                                                                    atatGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATaaaaa                                              <  2:147253/149‑1 (MQ=255)
                                                                                                     tatGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAAc                                             <  3:103739/149‑1 (MQ=255)
                                                                                                       tGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAAcat                                           >  3:33380/1‑149 (MQ=255)
                                                                                                            aGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAAcattcatt                                      >  1:657097/1‑149 (MQ=255)
                                                                                                              tCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCAtttt                                    >  2:595902/1‑149 (MQ=255)
                                                                                                                 ccGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCAttttttt                                 >  3:678375/1‑149 (MQ=255)
                                                                                                                   gTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTaa                               >  1:120226/1‑149 (MQ=255)
                                                                                                                   gTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTaa                               >  2:499794/1‑149 (MQ=255)
                                                                                                                    tCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACACTCATTTTTTTaac                              <  4:692972/149‑2 (MQ=255)
                                                                                                                       aTTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGt                           >  2:545985/1‑149 (MQ=255)
                                                                                                                       aTTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGt                           >  2:103471/1‑149 (MQ=255)
                                                                                                                        ttGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAAc                                                                                     <  4:84761/90‑1 (MQ=255)
                                                                                                                        ttGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAAc                                                                                     >  3:84761/1‑90 (MQ=255)
                                                                                                                            gtACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCgt                      >  3:233752/1‑149 (MQ=255)
                                                                                                                             tACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCgtg                     >  1:656404/1‑149 (MQ=255)
                                                                                                                                aGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCg                  >  4:188893/1‑149 (MQ=255)
                                                                                                                                 gCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCgg                 >  2:649485/1‑149 (MQ=255)
                                                                                                                                   acaATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGt               >  4:20196/1‑149 (MQ=255)
                                                                                                                                      ctttttGCGAGATTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACt            >  4:355075/2‑149 (MQ=255)
                                                                                                                                      ctttttGCGAGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACt            >  4:131757/2‑149 (MQ=255)
                                                                                                                                       tttttGCGAGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATaaaaa                                              <  3:795898/114‑1 (MQ=255)
                                                                                                                                       tttttGCGAGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATaaaaa                                              >  4:795898/1‑114 (MQ=255)
                                                                                                                                       tttttGCGAGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCggg                <  3:536352/144‑1 (MQ=255)
                                                                                                                                       tttttGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTg           >  1:815298/1‑149 (MQ=255)
                                                                                                                                       tttttGCGAGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCGAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCggg                >  4:536352/1‑144 (MQ=255)
                                                                                                                                            gCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTAc      <  1:73447/149‑1 (MQ=255)
                                                                                                                                             cGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTAcc     >  1:537543/1‑149 (MQ=255)
                                                                                                                                                aGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCaaa  <  4:488423/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
ATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGCTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAA  >  NC_000913/2715114‑2715406

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