Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,097,144 (C)6→7 intergenic (+315/+421) 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,1381.C97.1% 105.9 / ‑3.3 34intergenic (+309/+427)ycdU/serXputative inner membrane protein/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (13/20);  total (13/21)
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

CGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCAC  >  NC_000913/1097000‑1097266
                                                                                                                                             |                                                                                                                                
cGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑t                                                                                                                                   <  2:1093007/139‑1 (MQ=32)
     taGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCTCTCTTGCTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTGCTCA‑‑Gccccc                                                                                                                             <  2:66009/138‑1 (MQ=11)
     gAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑Tccccc                                                                                                                             <  2:1458060/139‑1 (MQ=33)
       gCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑Tccccccc                                                                                                                           >  2:295233/1‑138 (MQ=17)
       gCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑Tccccccc                                                                                                                           >  2:1165979/1‑138 (MQ=18)
        ctctTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑Tcccccccg                                                                                                                          <  1:143793/139‑3 (MQ=17)
        ctctTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑Tcccccccg                                                                                                                          <  2:1039374/139‑3 (MQ=17)
        ctctTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑Tcccccccg                                                                                                                          <  1:769959/139‑3 (MQ=17)
        ctctTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑Tcccccccg                                                                                                                          <  1:630578/139‑3 (MQ=18)
          cttcttTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑Tcccccccggt                                                                                                                        <  2:745929/139‑5 (MQ=17)
           ttcttTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑Tccccc                                                                                                                             >  1:1031424/1‑133 (MQ=32)
           ttcttTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑Tccccc                                                                                                                             <  2:1031424/133‑1 (MQ=32)
               ttAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGtgtgtg                                                                                                                    <  2:1374591/139‑1 (MQ=17)
                  aaTATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCaa                                                                                                                 <  1:1444827/139‑1 (MQ=17)
                       ggCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATacac                                                                                                            >  2:994261/1‑139 (MQ=21)
                        gCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACg                                                                                                           <  2:980678/139‑1 (MQ=25)
                            tGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGaaaa                                                                                                       >  1:1175334/1‑139 (MQ=21)
                                    ggATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGccc                                                                                               >  2:750180/1‑139 (MQ=255)
                                     gATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCg                                                                                              <  1:994261/139‑1 (MQ=21)
                                       ttGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTa                                                                                            <  1:464306/139‑1 (MQ=255)
                                          aCTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGaaaaaa                                                                                                     >  1:1009831/1‑127 (MQ=255)
                                          aCTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGaaaaaa                                                                                                     <  2:1009831/127‑1 (MQ=25)
                                                   gctcgcCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTTGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAg                                                                                <  1:362577/139‑1 (MQ=255)
                                                     tcgcCCTGCGGGCAGCCCGCTTACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGct                                                                              <  1:192419/139‑1 (MQ=255)
                                                     tcgcCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGct                                                                              <  2:1036373/139‑1 (MQ=25)
                                                     tcgcCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGct                                                                              <  1:295233/139‑1 (MQ=25)
                                                     tcgcCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGct                                                                              <  1:1165979/139‑1 (MQ=255)
                                                     tcgcCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGct                                                                              <  1:750180/139‑1 (MQ=255)
                                                      cgcCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCAACCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGct                                                                              <  2:462662/139‑1 (MQ=255)
                                                        cccTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCtctt                                                                           >  2:566417/1‑139 (MQ=255)
                                                                 cAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAata                                                                  <  2:1175334/139‑1 (MQ=25)
                                                                   gCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGaaaaaa                                                                                                     <  2:1201860/102‑1 (MQ=255)
                                                                   gCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGaaaaaa                                                                                                     >  1:1201860/1‑102 (MQ=255)
                                                                         tCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑Tcccccccgg                                                                                                                         >  1:1363388/1‑72 (MQ=21)
                                                                               cGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAggggg                                                    >  1:523936/1‑139 (MQ=21)
                                                                                 ttCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGctct                                                                            >  1:1403463/1‑113 (MQ=255)
                                                                                 ttCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGctct                                                                            <  2:1403463/113‑1 (MQ=255)
                                                                                                      tGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGct                                                                              <  2:468850/90‑1 (MQ=255)
                                                                                                      tGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGct                                                                              >  2:515755/1‑90 (MQ=255)
                                                                                                      tGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGct                                                                              >  1:468850/1‑90 (MQ=255)
                                                                                                      tGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGct                                                                              <  1:515755/90‑1 (MQ=255)
                                                                                                          gCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGt                                                           <  2:745468/105‑1 (MQ=255)
                                                                                                          gCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGt                                                           >  1:745468/1‑105 (MQ=255)
                                                                                                               ttGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAg                    <  1:566417/139‑1 (MQ=25)
                                                                                                                  tCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGccc                 >  1:902998/1‑139 (MQ=25)
                                                                                                                          cGGTCGGGGCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTg         <  2:902998/139‑1 (MQ=17)
                                                                                                                                 ggCTTCTCA‑‑TCCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCAc  >  2:1406453/1‑139 (MQ=18)
                                                                                                                                             |                                                                                                                                
CGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGACCCCCGGTCGGGGCTTCTCA‑‑TCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGCCCGCTCACTGCGTTCAC  >  NC_000913/1097000‑1097266

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: