Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,011,323 +T intergenic (‑258/+192) PP_0871 ← / ← SBRG_Ecoli_galM glycine betaine/carnitine/choline ABC transporter permease/SBRG_Ecoli_galM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,011,3231.T97.4% 126.4 / ‑3.0 39intergenic (‑258/+192)PP_0871/SBRG_Ecoli_galMglycine betaine/carnitine/choline ABC transporter permease/SBRG_Ecoli_galM
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base T (17/21);  total (18/21)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GCTGATGAAACCAGTCCCCAACCGAAGCGGGTTGCAGGGGAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCG  >  NC_002947/1011180‑1011457
                                                                                                                                                 |                                                                                                                                      
gCTGATGAAACCAGTCCCCAACCGAAGCGGGTTGCAGGGGAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCaa                                                                                                                                   <  1:391291/150‑1 (MQ=255)
 cTGATTAAACCAGTCCCCAACCGAAGCGGGTTGCAGGGGATTGGACGTCTCCGTGGGCGTAAAGGTTCCAATCTGGGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGCTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCaaa                                                                                                                                  <  1:85361/150‑1 (MQ=255)
  tgatgaAACCAGTCCCCAACCGAAGCGGGTTGCAGGGGAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTCTTGGGCGATTAAAACAGCGACTTGATGTCTCAAat                                                                                                                                 >  2:337639/1‑150 (MQ=255)
         aCCAGTCCCCAACCGAAGCGGGTTGCAGGGGAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCttt                                                                                                                          >  1:543184/1‑150 (MQ=255)
            atctcccAACCGAAGCGGGTTGCAGGGGAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTaa                                                                                                                       <  1:1245906/146‑1 (MQ=255)
                cccAACCGCAGCGGGTTGCAGGGGAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGGATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATCTCTCAAATCATCTTTTAATGat                                                                                                                   >  2:582973/1‑150 (MQ=255)
                      cGAAGCGGGTTGCAGGGGAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATaataat                                                                                                             >  1:918561/1‑150 (MQ=255)
                                  cAGGGGAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATAtt                                                                                                 <  1:1268943/150‑1 (MQ=255)
                                   aGGGGAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATta                                                                                                     <  2:646218/145‑1 (MQ=255)
                                   aGGGGAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATta                                                                                                     >  1:646218/1‑145 (MQ=255)
                                     gggAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGa                                                                                                                    <  2:187310/128‑1 (MQ=255)
                                     gggAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGa                                                                                                                    >  1:187310/1‑128 (MQ=255)
                                       gAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCcgc                                                                                            <  2:1279265/150‑1 (MQ=255)
                                         gTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGcga                                                                                          <  2:918561/150‑1 (MQ=255)
                                             aCGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGaa                                                                                      >  2:775635/1‑150 (MQ=255)
                                                  tcCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAAc                                                                                 >  1:860024/1‑150 (MQ=255)
                                                    cGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAg                                                                               >  2:1290941/1‑150 (MQ=255)
                                                             aaaGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCcg                                                                                             <  1:1261591/127‑1 (MQ=255)
                                                             aaaGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCcg                                                                                             >  2:1261591/1‑127 (MQ=255)
                                                                    cccATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATaa                                                               >  1:615323/1‑150 (MQ=255)
                                                                     ccATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAAc                                                              >  1:490770/1‑150 (MQ=255)
                                                                         cTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATa                                                          <  1:775635/150‑1 (MQ=255)
                                                                           gAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGATGTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGc                                                        <  1:337639/150‑1 (MQ=255)
                                                                              gcggcATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAAtt                                                                                                            >  2:657077/1‑95 (MQ=255)
                                                                              gcggcATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAAtt                                                                                                            <  1:657077/95‑1 (MQ=255)
                                                                              gcggcATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGaa                                                     >  2:1175530/1‑150 (MQ=255)
                                                                                ggcATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGaaaa                                                   <  2:543184/150‑1 (MQ=255)
                                                                                    tttGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGt                                               >  2:833825/1‑150 (MQ=255)
                                                                                       ggCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGgc                                            <  2:615323/150‑1 (MQ=255)
                                                                                            gtgtTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAgt                                       <  1:833825/150‑1 (MQ=255)
                                                                                            gtgtTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAgt                                       <  2:860024/150‑1 (MQ=255)
                                                                                                tAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATCACCATAGCGAAAATAGTGGCGCAGTGTaa                                   <  1:881580/150‑1 (MQ=255)
                                                                                                 aaTCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTCAg                                  <  1:1290941/150‑1 (MQ=255)
                                                                                                                ggTTTTTGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATAtt                   >  1:738114/1‑150 (MQ=255)
                                                                                                                   ttttGGGCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAg                <  2:738114/150‑1 (MQ=255)
                                                                                                                        ggCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTgg                                             >  1:99034/1‑116 (MQ=255)
                                                                                                                        ggCGAGAAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTgg                                             <  2:99034/116‑1 (MQ=255)
                                                                                                                           gagaAAAACAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTCGCAATATTGAGTTGCaaaa        >  1:66038/1‑149 (MQ=255)
                                                                                                                                 aaCAGCGACTTGATGTCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCg  <  1:1175530/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                      
GCTGATGAAACCAGTCCCCAACCGAAGCGGGTTGCAGGGGAGTGGACGTCTCCGTGGGCGTAAAGGTTCCCATCTGAGGCGGCATTTGGCCAGTGTTAATCGCTGGGTTGCGGGTTTTTGGGCGAGAAAAACAGCGACTTGATGCTCAAATCATCTTTTAATGATAATAATTCTCATTATATTGCCGCGACGAAGCAACAGCAATGCTTACGCATAACCATAGCGAAAATAGTGGCGCAGTGTAAGGTTGTTGTGAATATTGAGTTGCAAATATGTCG  >  NC_002947/1011180‑1011457

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