Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,097,209 T→G 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,2090TG92.9% 152.6 / ‑0.3 71intergenic (+380/+356)ycdU/serXputative inner membrane protein/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base G (33/32);  total (37/34)
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.98e-01

TGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCAC  >  NC_000913/1097112‑1097257
                                                                                                 |                                                
tGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGa                                                >  1:539648/1‑100 (MQ=25)
 gTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAg                                               >  2:56275/1‑100 (MQ=38)
  tCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAgg                                              >  2:2867000/1‑100 (MQ=255)
   cGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAggg                                             >  2:2936082/1‑100 (MQ=255)
    gaaccccGGTCGGCGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAACAAAGCCCGTACTCTCGTACGTGCTCCTCTTTAACTACGGCGGTCAg                                               >  1:974368/4‑97 (MQ=1)
    gACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAgggg                                            >  2:90369/1‑100 (MQ=38)
    gACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAgggg                                            <  2:2844083/100‑1 (MQ=255)
     aCCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAggggg                                           >  2:1255245/1‑100 (MQ=25)
      cccccGGTCGGGTCTTCTCATCCCCCCGGTGTGTGAAATACACGAAAAAAAAGCCCGTACTTTAGTACGAGCTCTTCTTTAAATATGGCGGGGAgggggg                                          <  2:2443793/100‑1 (MQ=11)
      cccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAgggggg                                          >  1:2163159/1‑100 (MQ=25)
       ccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCgggggggggggg                                         >  1:2196090/1‑100 (MQ=16)
       ccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAggggggg                                         >  1:999117/1‑100 (MQ=18)
        cccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGa                                        >  1:938393/1‑100 (MQ=18)
        cccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGa                                        >  2:452165/1‑100 (MQ=18)
         ccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGAt                                       <  1:2936035/100‑1 (MQ=17)
          cGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGAtt                                      <  2:138786/100‑1 (MQ=17)
          cGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGACGGGGAGGGGGGGAtt                                      >  2:588061/1‑100 (MQ=12)
           ggTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTg                                     <  1:2046052/100‑1 (MQ=17)
           ggTCGGGGCTTCTCATCCCCCCGGTGGGTGAAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTg                                     <  2:3115775/100‑1 (MQ=255)
             tCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGAc                                   >  1:2839224/1‑100 (MQ=18)
             tCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGAc                                   <  2:2743875/100‑1 (MQ=17)
             tCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGAc                                   <  2:938405/100‑1 (MQ=18)
              cGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAgggggg                                          >  1:725174/1‑92 (MQ=17)
              cGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACt                                  <  2:1747269/100‑1 (MQ=18)
              cGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACt                                  >  2:1424541/1‑100 (MQ=18)
              cGGGGCTTCTCATCCCCCCGGTGTGTACAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGAc                                   <  2:1478148/99‑1 (MQ=12)
               ggggCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGCGCTCTTCTTTAAATATGGAGGGGAGGGGGGGATTGACtc                                 >  2:973585/1‑100 (MQ=255)
               ggggCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtc                                 <  1:1490878/100‑1 (MQ=18)
                gggCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtcg                                >  1:2752286/1‑100 (MQ=17)
                 ggCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtcgc                               <  2:1829064/100‑1 (MQ=17)
                  gCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtcgct                              <  1:750480/100‑1 (MQ=18)
                   cTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtcgctt                             <  2:903532/100‑1 (MQ=18)
                   cTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtcgctt                             >  2:1599819/1‑100 (MQ=18)
                   cTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtcgctt                             >  2:1525285/1‑100 (MQ=18)
                     tctcATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtcgcttcg                           >  2:3044955/1‑100 (MQ=18)
                     tctcATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGGACTCTCATACGAGCTCTTCTTTAAAGATAGCGGGGCGGGGGGGACTGGCgcgcttc                            >  2:2929094/1‑99 (MQ=11)
                      tgagaagccccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtcgcttcgc                          >  2:1009864/8‑100 (MQ=255)
                      ctcggccccccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtcg                                <  1:2544474/89‑1 (MQ=255)
                      cgcgggcccccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtcgcttc                            <  2:3354959/92‑1 (MQ=255)
                       tcATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtcgcttcgct                         <  1:1298037/100‑1 (MQ=18)
                        cATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTtcgctc                        >  1:2238806/1‑100 (MQ=18)
                         aTCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTtcgctcg                       <  2:2163190/100‑1 (MQ=18)
                         aTCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTtcgctcg                       <  2:3441295/100‑1 (MQ=17)
                          tCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTtcgctcgc                      <  2:2839268/100‑1 (MQ=18)
                          tCCCCCCGGTGTGTGCAATACACAAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACtcgcttcgc                          >  1:3297046/1‑96 (MQ=255)
                           ccccccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACtcgcttcgc                          <  1:2923529/95‑1 (MQ=32)
                           ccccccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGcc                     >  2:2923575/1‑100 (MQ=17)
                           ccccccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGcc                     >  2:2998404/1‑100 (MQ=18)
                           ccccccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGcc                     <  2:1000716/100‑1 (MQ=18)
                           ccccccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGcc                     <  2:3423984/100‑1 (MQ=18)
                            cccccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGccc                    <  1:193460/100‑1 (MQ=17)
                            cccccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGccc                    >  2:775223/1‑100 (MQ=18)
                            cccccGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGcc                     >  2:3057041/1‑99 (MQ=18)
                            cccccGGTGTGTGCAATACACAAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGccc                    <  2:3297097/100‑1 (MQ=255)
                                 ggTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCgg               <  2:1435626/100‑1 (MQ=17)
                                  gtgtgtGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCggg              <  2:2224491/100‑1 (MQ=18)
                                    gtgtgAAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCa            <  2:2180861/100‑1 (MQ=12)
                                     tgtgCAATACACAACAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAg           <  2:1951993/100‑1 (MQ=255)
                                      gtgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGc          >  2:2486877/1‑100 (MQ=18)
                                      gtgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGc          >  2:230328/1‑100 (MQ=18)
                                      gtgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGc          >  1:2485026/1‑100 (MQ=18)
                                       tgCAATACACGAAAAAAAAGCCCGTTCTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGcc         <  2:1112559/100‑1 (MQ=255)
                                       tgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGcc         <  1:1255228/100‑1 (MQ=32)
                                       tgCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGcc         >  1:194727/1‑100 (MQ=18)
                                         cAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTTACTCGCTTAGCTCGcccggggg               >  1:2958527/1‑87 (MQ=11)
                                          aaTACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGc      >  2:3203559/1‑100 (MQ=17)
                                          aaTACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGccgcccgc      >  1:100410/1‑100 (MQ=12)
                                           aTACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCt     <  2:2300462/100‑1 (MQ=17)
                                            tACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTc    <  1:2866954/100‑1 (MQ=17)
                                            tACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTc    >  1:1748462/1‑100 (MQ=17)
                                            tACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTc    <  1:1134870/100‑1 (MQ=17)
                                              cacGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGCGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCAc  >  1:2072721/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: