Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,489,089 A→T 9.2% A325A (GCT→GCA)  yjgR ← DUF853 domain‑containing protein YjgR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,489,0890AT 9.2% 95.7 / 4.2 44A325A (GCT→GCA) yjgRDUF853 domain‑containing protein YjgR
Reads supporting (aligned to +/- strand):  ref base A (14/25);  new base T (2/2);  total (17/27)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-01

CTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGTA  >  NC_000913/4488952‑4489231
                                                                                                                                         |                                                                                                                                              
cTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCg                                                                                                                                           <  1:217050/143‑1 (MQ=255)
    gtttCGCCTCCAGAAAAGAGTGCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAATCGCCGGATTGGCCCTCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCaaa                                                                                                                                       <  2:517849/142‑1 (MQ=255)
         gCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATCGGCCCGCATGCTGCGCGCCCCTGCCTGCACTGCTTTCCGACCTTTCGGCGTAACAGCCCGCGAACCacg                                                                                                                                  >  2:238512/1‑141 (MQ=255)
              cAGAAAAACGCTCAGCGCCTCCCCGGTCCCCAGTTCCTCAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGaa                                                                                                                             <  2:57107/143‑1 (MQ=255)
               aGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTCTTGGCGTGAAAGCCCGCAGAGCCGGTTGAAc                                                                                                                            >  2:474415/1‑143 (MQ=255)
                         tCAACGCCCCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTGTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCg                                                                                                                  <  2:179023/143‑1 (MQ=255)
                           cccGCCTCCCCGCTCCCCAGTTCCTCATTCGCCTTTTCGGTACCCAACGCCCGATTGTCCCTTTTGTTTTGCGCCGCTGCCTTCACTGCTTTCTGCCCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAg                                                                                                                <  2:180318/141‑1 (MQ=255)
                           atcgccgcGTCGGTGCCCAGTTCCTGAATCGCCTTGTCGGTATCAACCGGCGGCTTGGCCCGCATGGTTTGAGCCGCTGCCTTCACTGCTTTTTGATCTTTTGGCGTAGAAGCCCGCAAAGCATGTTGAACGCGATTACCgcg                                                                                                                <  2:180293/141‑3 (MQ=255)
                             cgtcgcGCTGGTGCCGCGTTCCGGCTTCGCCTTTTCGGTATCAACCGCCGGATGGGTCCGCTTGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTACAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCt                                                                                                              <  2:192953/138‑1 (MQ=255)
                               ccgcGCCTGTGCCTGGGTCATGAGTCGCTTTTCCGTTATCAAACGCCGGGTTGGCCCGTATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTgt                                                                                                            <  2:20655/140‑2 (MQ=255)
                                 tCGCCGGGGCCCAGTGCCTGAATCGCCGTTTCGGGATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGccc                                                                                                          <  1:6325/143‑1 (MQ=255)
                                   gCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAATGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGa                                                                                                        >  1:217981/1‑143 (MQ=255)
                                   gCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAATGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGa                                                                                                        >  1:502420/1‑143 (MQ=255)
                                    ccGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCCTGTTGAACGCGATTACCGAGCTGCCCGAg                                                                                                       >  2:12771/1‑143 (MQ=255)
                                       gTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGCTTACCGACCTGCCCGAGcac                                                                                                    >  1:43823/1‑143 (MQ=255)
                                                tCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCgg                                                                                           <  1:371149/143‑1 (MQ=255)
                                                     acgcGCCGTTTCGGGATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAatat                                                                                      <  2:487427/140‑1 (MQ=255)
                                                             tttCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACggg                                                                              <  1:272053/143‑1 (MQ=255)
                                                                    gtCAGACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCGTCACTGCTCTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCg                                                                       <  1:76962/142‑1 (MQ=255)
                                                                     tCACACCCCGGGTTGACCCGCATGGTTTGCGCCGCTGCCTTCCCTGCTTTCTCATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGa                                                                      <  2:401043/143‑1 (MQ=255)
                                                                       aaaCGCCGGATTGGCCGGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGaaa                                                                    <  1:24804/143‑1 (MQ=255)
                                                                            ccGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAAc                                                               <  1:392847/143‑1 (MQ=255)
                                                                             cGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTACAAGCCCGCAAAGCATGTTGAACGCGATTACCGCGCTGCCCGAGCCCATTGTCTGGGATATCAGACGGGTTGTGCGAAACGAAcc                                                              >  1:17540/1‑143 (MQ=255)
                                                                                 ttGGCCCGCATGGTTGGCGCCGCTGCCTGCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGAc                                                          <  1:187967/143‑1 (MQ=255)
                                                                                    gCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGCTCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGCTTACCGAGCTGCCCGAGCACCTTATCCGGGATGTCCGACGGGTTTTGCGAGACGACCCCGACGcc                                                       >  2:487141/1‑143 (MQ=255)
                                                                                    ccccGCATGGTTTGCGCCGCTGCCTTCACTCCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGcc                                                       <  2:8990/142‑1 (MQ=255)
                                                                                             gTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCttt                                              >  1:249046/1‑143 (MQ=255)
                                                                                              tttGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAACGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCtttt                                             >  2:466126/1‑143 (MQ=255)
                                                                                              tttGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCtttt                                             >  2:311928/1‑143 (MQ=255)
                                                                                                tGCGCCGCTGCCTTTACTCCGTTCTGATCGTGTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGa                                           <  2:487561/143‑1 (MQ=255)
                                                                                                        tGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATaa                                   <  1:369726/143‑1 (MQ=255)
                                                                                                            ttCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCt                               >  1:246688/1‑143 (MQ=255)
                                                                                                              cACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGATTTGCCCTAGCTCATTATCCGTAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTa                             <  1:51611/143‑1 (MQ=255)
                                                                                                               aCTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGCACTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGACAAGTACAATAGGCCCCCGCCTATGGACCGACTAAGcataat                            >  2:500979/1‑138 (MQ=255)
                                                                                                                    tttCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCt                       >  1:217982/1‑143 (MQ=255)
                                                                                                                       cggATCCTCTGCCGTGAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCt                    <  2:267375/141‑1 (MQ=255)
                                                                                                                       cTGATCTTTTGTCTTAAAAGCCCGCAAAGTATGTTGAACGCGATTGCCGAGCTGGCCGAGCACATGATCCGGAATAGCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCt                    <  1:238512/143‑1 (MQ=255)
                                                                                                                        tGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTCTAACCTGCTc                   >  2:383449/1‑143 (MQ=255)
                                                                                                                        tGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTc                   >  2:359799/1‑143 (MQ=255)
                                                                                                                          aTCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATCTAAGACGGGCTTCGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGGATAAGCCTTATCACCTGCTCac                 >  1:238041/1‑142 (MQ=255)
                                                                                                                            cTTTTGGCGTAACAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATc               <  1:446130/143‑1 (MQ=255)
                                                                                                                            cTTTTGGCGCAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATc               <  1:247191/143‑1 (MQ=255)
                                                                                                                               ttGGCGTAAATGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTa            <  2:217981/143‑1 (MQ=255)
                                                                                                                               ttGGCGTAAATGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTa            <  2:502420/143‑1 (MQ=255)
                                                                                                                               gtggctgaaaaGCCCGCAAAGCATGTTGAACGCGATTCCCGAGCTGCGCGAGCACATTATCCGGAATATCAGACGTGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTa            <  2:514430/136‑1 (MQ=255)
                                                                                                                                    ggcaaaGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCcag       <  1:150434/140‑1 (MQ=255)
                                                                                                                                      aaaaGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCcagcc     >  1:426309/1‑142 (MQ=255)
                                                                                                                                         aGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGTa  <  1:386651/143‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                              
CTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGTA  >  NC_000913/4488952‑4489231

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