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

ATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAA  >  NC_000913/2715114‑2715406
                                                                                                                                                |                                                                                                                                                     
aTAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑Tc                                                                                                                                                  >  1:267258/1‑149 (MQ=255)
       tgCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAGCT‑TCGCaaaaa                                                                                                                                           <  2:162129/149‑1 (MQ=255)
       tgCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGACACCGTCGATTGTGTACAGCACAATTTTTGCGAGCT‑TCTCaaaaa                                                                                                                                           <  4:238296/149‑1 (MQ=255)
                aaCTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCctc                                                                                                                                  >  4:103739/1‑149 (MQ=255)
                   ttttttCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTc                                                                                                                               >  3:625552/1‑149 (MQ=255)
                   ttttttCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTc                                                                                                                               >  3:796310/1‑149 (MQ=255)
                           ttCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACgg                                                                                                                       >  4:164213/1‑149 (MQ=255)
                                 ggTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTaa                                                                                                                 <  1:318592/149‑1 (MQ=255)
                                       gCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGaa                                                                                                                                       >  2:717346/1‑121 (MQ=255)
                                       gCACATCCGTCGGACAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGaa                                                                                                                                       <  1:717346/121‑1 (MQ=255)
                                            tCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATaaaa                                                                                                      <  3:480758/149‑1 (MQ=255)
                                              cGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACg                                                                                                    <  2:267258/149‑1 (MQ=255)
                                                             tcAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCCGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAAc                                                                                     >  2:73447/1‑149 (MQ=255)
                                                               gaaTAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGTACAATTTTTGCGAAGC‑TCGCaaaaa                                                                                                                                           >  4:693501/2‑93 (MQ=255)
                                                               aaaTAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCaaaaa                                                                                                                                           <  3:693501/93‑1 (MQ=255)
                                                                      tGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCt                                                                            >  1:147253/1‑149 (MQ=255)
                                                                       ggACTGATAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGCCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTa                                                                           <  3:688973/149‑1 (MQ=255)
                                                                           tGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGa                                                                       >  3:488423/1‑149 (MQ=255)
                                                                              aaaCTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGaaaa                                                                    <  4:796310/149‑1 (MQ=255)
                                                                              aaaCTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGaaaa                                                                    <  3:423079/149‑1 (MQ=255)
                                                                                 cTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAatat                                                                 <  3:164213/149‑1 (MQ=255)
                                                                                 cTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAatat                                                                 <  2:444256/149‑1 (MQ=255)
                                                                                        tCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCt                                                          >  4:303645/1‑149 (MQ=255)
                                                                                         cAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTa                                                         <  3:303645/149‑1 (MQ=255)
                                                                                                    atatGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATaaaaa                                              <  2:147253/149‑1 (MQ=255)
                                                                                                     tatGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAAc                                             <  3:103739/149‑1 (MQ=255)
                                                                                                       tGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAAcat                                           >  3:33380/1‑149 (MQ=255)
                                                                                                            aGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAAcattcatt                                      >  1:657097/1‑149 (MQ=255)
                                                                                                              tCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCAtttt                                    >  2:595902/1‑149 (MQ=255)
                                                                                                                 ccGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCAttttttt                                 >  3:678375/1‑149 (MQ=255)
                                                                                                                   gTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTaa                               >  1:120226/1‑149 (MQ=255)
                                                                                                                   gTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTaa                               >  2:499794/1‑149 (MQ=255)
                                                                                                                    tCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACACTCATTTTTTTaac                              <  4:692972/149‑2 (MQ=255)
                                                                                                                       aTTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGt                           >  2:545985/1‑149 (MQ=255)
                                                                                                                       aTTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGt                           >  2:103471/1‑149 (MQ=255)
                                                                                                                        ttGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAAc                                                                                     <  4:84761/90‑1 (MQ=255)
                                                                                                                        ttGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAAc                                                                                     >  3:84761/1‑90 (MQ=255)
                                                                                                                            gtACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCgt                      >  3:233752/1‑149 (MQ=255)
                                                                                                                             tACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCgtg                     >  1:656404/1‑149 (MQ=255)
                                                                                                                                aGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCg                  >  4:188893/1‑149 (MQ=255)
                                                                                                                                 gCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCgg                 >  2:649485/1‑149 (MQ=255)
                                                                                                                                   acaATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGt               >  4:20196/1‑149 (MQ=255)
                                                                                                                                      ctttttGCGAGAT‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACt            >  4:355075/2‑149 (MQ=255)
                                                                                                                                      ctttttGCG‑AGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACt            >  4:131757/2‑149 (MQ=255)
                                                                                                                                       tttttGCGAGCT‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATaaaaa                                              <  3:795898/114‑1 (MQ=255)
                                                                                                                                       tttttGCGAGCT‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATaaaaa                                              >  4:795898/1‑114 (MQ=255)
                                                                                                                                       tttttGCGAGCT‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCggg                <  3:536352/144‑1 (MQ=255)
                                                                                                                                       tttttGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTg           >  1:815298/1‑149 (MQ=255)
                                                                                                                                       tttttGCG‑AGCTTCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCGAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCggg                >  4:536352/1‑144 (MQ=255)
                                                                                                                                            gCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTAc      <  1:73447/149‑1 (MQ=255)
                                                                                                                                             cGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTAcc     >  1:537543/1‑149 (MQ=255)
                                                                                                                                                aGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCaaa  <  4:488423/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
ATAAATGTGCGGCAGCAACTTTTTTCCTTCGCGGGTAAGGCACATCCGTCGGCCAATTTTCTCAAATAACTGGACTGAAAACTCCTGCTCAAGCTGCTGAATATGAAAAGTCACCGTCGATTGTGTACAGCACAATTTTTGCGAAGC‑TCGCAAAAAGGAACCCTCTTCCACCACGGTTTTAAGCGTAATAAAACGGCGCAGATCCATAACCCCAAACCTATCGAAAATATCGAATCTAGAATATAAAAACATTCATTTTTTTAAATGTTCCGTGTCGGGTACTGTCTACCAAA  >  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: ‑