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,2090TG58.4% 4.1 / 49.0 48intergenic (+380/+356)ycdU/serXputative inner membrane protein/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base T (10/10);  new base G (17/11);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 5.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCAC  >  NC_000913/1097112‑1097257
                                                                                                 |                                                
tGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGa                                                >  1:2757249/1‑100 (MQ=255)
   cGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAAAACACGAAAAAAAAGCCCGTACTTTCGTACGCGATCTTCTTTACATATGGCGGGGAg                                               >  2:714635/1‑98 (MQ=11)
    gACCCCCGGTCGGGGCTTCTCATCCCCCGGGTGTGTGCAATACAAGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAgggg                                            <  2:2530784/100‑1 (MQ=255)
    gACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAgggg                                            >  2:2108402/1‑100 (MQ=255)
    gACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAgggg                                            >  2:3268313/1‑100 (MQ=255)
      cccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAgggggg                                          <  2:82156/100‑1 (MQ=33)
      cccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAgggggg                                          >  2:1393426/1‑100 (MQ=25)
      cccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAgggggg                                          >  2:1832309/1‑100 (MQ=25)
      cccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAgggggg                                          <  2:736848/100‑1 (MQ=25)
       ccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAggggggg                                         <  2:3074741/100‑1 (MQ=18)
         ccGGTCGGGGTTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGa                                        <  2:2617123/99‑1 (MQ=12)
           ggTCGGGGCTGCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTg                                     <  2:790312/100‑1 (MQ=255)
            gagagggCTTCGGAGCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTg                                     <  1:3117391/95‑1 (MQ=11)
             tCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGAc                                   <  1:2541246/100‑1 (MQ=32)
             tCGGGGCTGCTCATCCCCCCGGTGTGTGAAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGAc                                   <  1:853991/100‑1 (MQ=255)
              cGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTg                                     <  1:3029968/97‑1 (MQ=32)
               ggggCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGTGGGGGGGATTg                                     >  1:137297/1‑96 (MQ=16)
               ggggCTTCTCATCCCCCCGGTGTGTCCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtc                                 <  2:2146157/100‑1 (MQ=12)
                   cTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACtcgctt                             >  1:1115077/1‑100 (MQ=33)
                   cTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtcgctt                             >  2:2129750/1‑100 (MQ=18)
                     tctcATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtcgctt                             >  1:747823/1‑98 (MQ=18)
                     tcgcATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtcgcttc                            <  2:1588211/96‑1 (MQ=255)
                       tcATCCCCCGGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtcgcttcgct                         <  2:2940406/100‑1 (MQ=12)
                        cATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACtcgcttcgc                          >  1:1887463/1‑98 (MQ=33)
                        cATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACtcgcatatcgg                        >  1:2846337/1‑93 (MQ=255)
                        cATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTtcgctc                        >  2:1157717/1‑100 (MQ=33)
                         aTCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTtcgctcg                       >  2:123476/1‑100 (MQ=18)
                         aTCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCGCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTtcgctcg                       >  1:1639578/1‑100 (MQ=255)
                          tCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTtcgctcgc                      <  1:1410544/100‑1 (MQ=33)
                            cccccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTtcgctcgc                      >  2:1250308/1‑98 (MQ=33)
                            cccccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGccc                    >  2:1172806/1‑100 (MQ=33)
                            cccccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGccc                    >  2:670414/1‑100 (MQ=33)
                            cccccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGGTCGcc                     >  1:259518/1‑99 (MQ=12)
                             ccccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCt                   >  2:951402/1‑100 (MQ=33)
                              cccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTg                  >  2:101867/1‑100 (MQ=18)
                               ccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGc                 <  2:2443275/100‑1 (MQ=32)
                                  gtgtgtGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCggg              <  2:923055/100‑1 (MQ=33)
                                   tgtgtgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGc             <  2:2252112/100‑1 (MQ=32)
                                   tgtgtgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGc             <  1:2129735/100‑1 (MQ=18)
                                   tgtgtgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGc             >  1:103267/1‑100 (MQ=18)
                                   gatctgCAATACACGAAAAAAAAGCCCGTTCTTTCATACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCgg               <  1:2284379/94‑1 (MQ=255)
                                    gtgtgCAATACAGGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCa            >  2:2220620/1‑100 (MQ=255)
                                      gtgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGc          >  2:3266779/1‑100 (MQ=33)
                                      gtgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGc          >  2:1926191/1‑100 (MQ=33)
                                      gtgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGc          <  2:2265092/100‑1 (MQ=18)
                                        gaaaTACACGAAACAAAAGCCCGTACTTTCGAACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCggggatccg        >  2:2284395/3‑94 (MQ=11)
                                             acacGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCa   <  1:3071866/100‑1 (MQ=31)
                                              cacGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCAc  >  2:3079438/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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: