Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,518,275 T→G 8.9% intergenic (‑44/‑192) alaW ← / → yfeC tRNA‑Ala/putative DNA‑binding transcriptional regulator YfeC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,518,2750TG 8.9% 126.4 / 2.2 48intergenic (‑44/‑192)alaW/yfeCtRNA‑Ala/putative DNA‑binding transcriptional regulator YfeC
Reads supporting (aligned to +/- strand):  ref base T (27/16);  new base G (2/3);  total (29/19)
Fisher's exact test for biased strand distribution p-value = 3.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAA  >  NC_000913/2518136‑2518410
                                                                                                                                            |                                                                                                                                       
aGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCGGACCTCTGGCATGCCATGCAAGCGCCCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCg                                                                                                                                      <  1:256450/143‑1 (MQ=255)
      cTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGaa                                                                                                                                <  1:797306/143‑1 (MQ=255)
       tGGGTGAGAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAAt                                                                                                                               <  2:520655/143‑1 (MQ=255)
       tGGGGGCAACATTTGGTGGAGCTAAGCGGGATCGAACCGCTGCCCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAAt                                                                                                                               <  2:112950/143‑1 (MQ=255)
        gggTGCCACATTTGGGGGAGCTAAGCGGGATCGAACCGCTGACCTCGTGCATGCCATGCAAGCGCTCCCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGGCGAGTTGACGGGCGGAATTCATATGATGCCGCCCgtcaa                                                                                                                              <  2:395549/143‑5 (MQ=18)
         ggTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAAtgt                                                                                                                             <  1:450040/143‑1 (MQ=255)
             cAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCCTCATATGAATTCCGCCCGAATGTGTCa                                                                                                                         >  1:651462/1‑143 (MQ=255)
             cAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCa                                                                                                                         >  2:382380/1‑143 (MQ=255)
              acaggttGGTCGGGGTAAGGGGCCTCGAACCGCTGACGACTTGCATGCCTTCCAAGCGCTCTTCCAGCTGCGCGTTAGCCCCGTCACGTAAAGCTTGACGCGTCGTAAGTCACGATTTTATGAATTCCGCCCGAATGTGTCaa                                                                                                                        <  2:209497/138‑1 (MQ=255)
                 aTTTGGTGGAGCGAGTCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACgg                                                                                                                     <  2:511716/143‑1 (MQ=255)
                       tggAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGCATGTGTCAACGGCAAAta                                                                                                               >  2:519960/1‑142 (MQ=255)
                         gAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCCCCGGCCCATTGc                                                                                                             >  1:628167/1‑143 (MQ=255)
                           gCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCACGAATGTGTCAACGGCAAATTGaaa                                                                                                           >  1:230038/1‑143 (MQ=255)
                            cTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAAc                                                                                                          >  1:160068/1‑143 (MQ=255)
                             tAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAACCCTTGTCGAGTTGACGGGCGGCATCAAATGAATTCCGCCCGCATGTGGCACCGGCGAACTGccagg                                                                                                         >  2:196375/1‑139 (MQ=21)
                                        gAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCaa                                                                                              <  1:536351/143‑1 (MQ=255)
                                         aaCCGCTGACCTCTTGCATGCCATGCAAGCGCTCGCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAAt                                                                                             <  2:108983/143‑1 (MQ=255)
                                                    tctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGCTTCAATCGCTGAAAAAt                                                                                  >  2:187153/1‑143 (MQ=255)
                                                    tctTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGGAAAGCTTGTCGAGTTGAGGGGCGGCATCATCTGAATTCCGCCCGAATGTGTCAACGGCAAATTTCAACACGGAGCTTCAACCGCCGAAAAAt                                                                                  >  2:136283/1‑143 (MQ=37)
                                                      ttGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCACCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGGATGTGTCCACGGCCATTTGCATAGTGTCGTTTCCCTCGCTTAACACTCa                                                                                >  1:15615/1‑143 (MQ=25)
                                                             ccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGAATTCATATGATGCCGCCCg                                                                                                                                  <  1:632074/86‑1 (MQ=25)
                                                             ccatgcAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGAATTCATATGATGCCGCCCg                                                                                                                                  >  2:632074/1‑86 (MQ=25)
                                                                      cgcTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAAGGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAGATGAAGAAttt                                                                >  1:382900/1‑143 (MQ=255)
                                                                       gcgctcCCAGCTGAGCTATAGCCCCGTCACGTAAAGCGTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAAtttt                                                               <  1:269440/140‑1 (MQ=255)
                                                                       gcTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAAtttt                                                               <  2:773628/143‑1 (MQ=255)
                                                                        ctctCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTg                                                              >  1:769139/1‑143 (MQ=255)
                                                                         tctcCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAAGCAGGCAAATGAACAATTTTgg                                                             >  2:708812/1‑143 (MQ=255)
                                                                            cccAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAAGGAACAATTTTGGaaa                                                          >  1:295989/1‑143 (MQ=255)
                                                                               aGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAAcc                                                       >  1:757847/1‑143 (MQ=255)
                                                                                        ataGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCAGCCCGAATGTGTCAACGGCACATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCAtt                                              >  1:224519/1‑143 (MQ=255)
                                                                                          ggCCCCGTCACGTAAAGCTGGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCa                                            <  2:190565/142‑1 (MQ=255)
                                                                                          aGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCa                                            <  1:493562/143‑1 (MQ=255)
                                                                                            ccccGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCagt                                          >  1:82693/1‑143 (MQ=255)
                                                                                            ccccGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAACCAGGCAAATGAAAAATTTTGGAACACCACTCGCATGCagt                                          >  2:476458/1‑143 (MQ=255)
                                                                                                        aaGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGt                              >  2:181119/1‑143 (MQ=255)
                                                                                                          gCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCCTGTCa                            >  2:76176/1‑143 (MQ=255)
                                                                                                            gtGTCGAGTTGGCGGGCGGCATCATATGAATGCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACg                          <  2:628167/142‑1 (MQ=255)
                                                                                                              gTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCCTGTCACggc                        >  1:343200/1‑142 (MQ=255)
                                                                                                                  aGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACACTTTTGGAAAACCGCTCGCATACAGTACTTATTCCTGTCACGGTTccc                    >  2:85286/1‑143 (MQ=255)
                                                                                                                  aGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTGAAATCGCTGAAAAATCAGGCAAAGGAAGAATTTTGGAAGACCGCTCGCATTAAGAGGGGATGCATGTGACGGtatac                    >  1:319581/1‑139 (MQ=255)
                                                                                                                   gTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCt                   <  1:533740/143‑1 (MQ=255)
                                                                                                                   gTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCt                   <  1:359633/143‑1 (MQ=255)
                                                                                                                   gTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCt                   >  1:702539/1‑143 (MQ=255)
                                                                                                                   gTTGACGCGCGGCATCATATGAATTCCGCCCGAATCTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCt                   <  2:424749/143‑1 (MQ=255)
                                                                                                                         gggCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGACAACTCAGGCAAATGAACAATTTTGGAAAACCGCACCCATTCAGCACGTATTCACGTCACGGTATCCTGTaacg             >  2:458979/1‑141 (MQ=255)
                                                                                                                           gcggaattcatATGATGCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGATTCCTGTAAAgtg           >  1:43071/8‑143 (MQ=255)
                                                                                                                                    atatGAATTCCGCACGAATGTGTAAACAGCAAAATGCAACGTGTAGTTTAAATCGCTGAAAAATAAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTAATGCACCGAATTCCTGAAAAGTAGAGAtata   >  2:217828/1‑143 (MQ=255)
                                                                                                                                     tatGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATaa  >  2:232627/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                       
AGTTTGCTGGGTGCAAAATTTGGTGGAGCTAAGCGGGATCGAACCGCTGACCTCTTGCATGCCATGCAAGCGCTCTCCCAGCTGAGCTATAGCCCCGTCACGTAAAGCTTGTCGAGTTGACGGGCGGCATCATATGAATTCCGCCCGAATGTGTCAACGGCAAATTGCAACGTGTAGTTTCAATCGCTGAAAAATCAGGCAAATGAACAATTTTGGAAAACCGCTCGCATTCAGTAGTTATTCATGTCACGGTTTCCTGTAAAGTGGTGTTATAA  >  NC_000913/2518136‑2518410

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