Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,097,209 T→G 100% intergenic (+380/+356) ycdU → / ← serX putative inner membrane protein/tRNA‑Ser

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,097,2090TG90.4% 74.3 / ‑1.6 42intergenic (+380/+356)ycdU/serXputative inner membrane protein/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base G (16/20);  total (19/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

TGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCAC  >  NC_000913/1097112‑1097257
                                                                                                 |                                                
tGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGa                                                >  1:1310460/1‑100 (MQ=255)
 gTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAACCCGTACTTTAATAAAAACTCTTCTTTAAATATGGCGGGGa                                                >  1:277621/1‑99 (MQ=11)
  tCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAgg                                              <  2:2220706/100‑1 (MQ=255)
   cGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAgg                                              >  2:1576338/1‑99 (MQ=38)
     aCCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAggggg                                           <  1:2155925/100‑1 (MQ=25)
      cccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAgggggg                                          >  1:2473229/1‑100 (MQ=25)
      cccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAgggggg                                          <  2:1963778/100‑1 (MQ=25)
        cccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGa                                        >  1:2108310/1‑100 (MQ=33)
        cccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAggggggg                                         >  1:2100116/1‑99 (MQ=18)
          cGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGAtt                                      >  1:2156894/1‑100 (MQ=18)
           ggTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTg                                     <  1:1663973/100‑1 (MQ=17)
           ggTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTg                                     <  2:2610149/100‑1 (MQ=18)
             tCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGAc                                   <  2:2853236/100‑1 (MQ=32)
              cGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACt                                  <  1:483662/100‑1 (MQ=17)
              cGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACt                                  <  1:2897272/100‑1 (MQ=17)
               ggggCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGCGGGGGGGATTGAc                                   >  2:536756/1‑98 (MQ=255)
                ggtcTTCTCGTCCCCCCGGTTTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGAtt                                      <  2:1008302/91‑1 (MQ=11)
                gggCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtc                                 >  1:2668282/1‑99 (MQ=18)
                  gCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtcgct                              <  1:1121858/100‑1 (MQ=17)
                  gCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtcgct                              >  2:158184/1‑100 (MQ=18)
                      ctcATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtcgcttcgc                          <  1:1509515/100‑1 (MQ=17)
                      cgcATCCACCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtcgct                              <  2:2432193/94‑1 (MQ=255)
                       tcATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACtcgcttcgct                         <  1:1639544/100‑1 (MQ=32)
                        cATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTtcgctc                        >  2:2798435/1‑100 (MQ=18)
                        cATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCGTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTtcgctc                        >  2:1940604/1‑100 (MQ=12)
                        cATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACCAGCTCTTCTTTAAATATGGCGGCGAGGGGGGGATTGACtaacttcgcta                        >  2:1639525/1‑99 (MQ=11)
                         ttCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTtcgctcg                       <  1:1146874/99‑1 (MQ=17)
                          tCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTtcgctcgc                      >  2:1056629/1‑100 (MQ=18)
                          tCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTtcgctcga                      >  2:1146863/1‑99 (MQ=18)
                            cccccGGTGTGTGCAATACACGAAAAAAAAGCCCGTTCTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGccc                    <  1:965849/100‑1 (MQ=255)
                            cccccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGccc                    >  1:2937449/1‑100 (MQ=32)
                            cccccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGccc                    <  1:430839/100‑1 (MQ=18)
                             acccGGTGTGTGCAATACACGAAAATAAAGCTCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGTAGGGGGGTATTGACtcgcttcg                           <  2:1001242/91‑1 (MQ=11)
                                   tttgtgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGc             <  1:1092097/98‑1 (MQ=17)
                                   tgtgtgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGc             <  2:2504963/100‑1 (MQ=17)
                                   tgtgtgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGc             <  1:2798467/100‑1 (MQ=17)
                                     tgtgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAg           <  1:1972420/100‑1 (MQ=17)
                                      gtgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGc          >  1:417706/1‑100 (MQ=18)
                                       tgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAg           <  1:1204699/98‑1 (MQ=17)
                                        gCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGccc        >  2:654788/1‑100 (MQ=18)
                                         ccgatcacGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCg       <  2:571359/95‑1 (MQ=11)
                                              cacGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCAc  >  1:1902975/1‑100 (MQ=11)
                                                                                                 |                                                
TGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCAC  >  NC_000913/1097112‑1097257

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