Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,184,585 G→C 17.5% intergenic (‑201/‑211) mltF ← / → purL membrane‑bound lytic murein transglycosylase F/phosphoribosylformylglycinamidine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,184,5850GC17.5% 81.8 / 15.3 40intergenic (‑201/‑211)mltF/purLmembrane‑bound lytic murein transglycosylase F/phosphoribosylformylglycinamidine synthase
Reads supporting (aligned to +/- strand):  ref base G (20/13);  new base C (4/3);  total (24/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.01e-01

GTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAG  >  NC_002947/1184445‑1184733
                                                                                                                                            |                                                                                                                                                      
gtcgccggtggtGGGGTTGGCGGATGGGTGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGATATATTAAGCGCCTGAGAGATCGAGCGCCGAGATGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCgcgggcgg                                                                                                                                                 <  8:104602/142‑1 (MQ=255)
gTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCgcgggc                                                                                                                                                   <  5:50512/146‑1 (MQ=255)
           cgtgggggggggAGGGGCGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGCTCGAGCGCCGCCCGCGCGGCGCTCGATCTca                                                                                                                                     <  5:353904/142‑1 (MQ=37)
            gtgGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTcac                                                                                                                                    >  7:333581/1‑149 (MQ=255)
                       aGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCGCCCAGgcgctgcat                                                                                                                         >  8:352215/1‑144 (MQ=255)
                                tGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAg                                                                                                                 <  8:333581/148‑1 (MQ=255)
                                  cATCAGGGCTGCGAGGTGGCCGGCTTGAGATATTCAGCGCCTGTTATAGCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCg                                                                                                              <  6:359557/149‑1 (MQ=255)
                                   aTCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCgg                                                                                                             >  8:304892/1‑149 (MQ=255)
                                      aGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGaca                                                                                                          >  8:337928/1‑149 (MQ=255)
                                        ggCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGAtc                                                                                                                                        <  3:222338/117‑1 (MQ=255)
                                        ggCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCCCGCGCGGCGCTCGAtc                                                                                                                                        >  4:222338/1‑117 (MQ=255)
                                           tGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTgc                                                                                                     >  1:49768/1‑149 (MQ=255)
                                           tGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTgc                                                                                                     >  6:187512/1‑149 (MQ=255)
                                              cAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTgcg‑‑gc                                                                                                <  7:337928/149‑1 (MQ=255)
                                                    gCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGa                                                                                          >  5:69479/1‑149 (MQ=255)
                                                            gagaTATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGtct                                                                                  >  2:358131/1‑149 (MQ=255)
                                                            gagaTATTCAGCGCCTGCGAGACCGAGCGCCGCCAGGTGCCCGCCTTAAGATATTCAGCCCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGtc                                                                                   <  8:60369/148‑1 (MQ=255)
                                                                tatTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAg                                                                              >  2:118168/1‑149 (MQ=255)
                                                                tatTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAg                                                                              <  6:69479/149‑1 (MQ=255)
                                                                    cAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCgg                                                                          >  6:259765/1‑149 (MQ=255)
                                                                        gcCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCACCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGacac                                                                      <  5:259765/149‑1 (MQ=255)
                                                                          cTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACct                                                                    >  4:118032/1‑149 (MQ=255)
                                                                                 aTCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACCTCTGaaaa                                                             >  1:2455/1‑149 (MQ=255)
                                                                                     aGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACctc                                                                   >  2:237503/1‑139 (MQ=255)
                                                                                     aGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACctc                                                                   <  1:237503/139‑1 (MQ=255)
                                                                                      gcgcCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGCCACCTGCG‑‑GCGCTCGCTCTCGTCTCAAGCCggtc                                                                        >  4:33922/1‑131 (MQ=255)
                                                                                            ccAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGCTCTCACAGGCGCAGAATGTGTCGAGACGGACACCTGCG‑‑GCGCTCTATCTCGTCTCAAGCCGGACACCTCTGAAAACTTACCCCGaa                                                  <  7:207091/149‑1 (MQ=255)
                                                                                                 tgtCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCg                                             >  7:30294/1‑149 (MQ=255)
                                                                                                  gtCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCgg                                            <  1:180171/149‑1 (MQ=255)
                                                                                                        cTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTccgc                                      >  5:382513/1‑149 (MQ=255)
                                                                                                         ttAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTccgcc                                     >  3:318916/1‑149 (MQ=255)
                                                                                                              atatTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCGCAGCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCgccgcc                                  <  1:118168/149‑1 (MQ=255)
                                                                                                                         ctgtGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCATAGGCGCTGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGc                     >  3:386219/4‑149 (MQ=37)
                                                                                                                         ctgtGAGATCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGCCACCTCTGAAACCTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTgc                      >  7:272048/4‑147 (MQ=37)
                                                                                                                            tGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCt                    >  5:170918/1‑147 (MQ=255)
                                                                                                                            tGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTa                  <  8:30294/149‑1 (MQ=255)
                                                                                                                              agaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGAcc                                                       >  4:99319/1‑110 (MQ=255)
                                                                                                                              agaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGAcc                                                       <  3:99319/110‑1 (MQ=255)
                                                                                                                                   gAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAg            <  5:187512/148‑1 (MQ=255)
                                                                                                                                    aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGt          >  6:349909/1‑149 (MQ=255)
                                                                                                                                    aGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGt          <  1:358131/149‑1 (MQ=255)
                                                                                                                                        ccgcGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGTCCACCAGCG‑‑GCGCTAGATATCGTCGCAAGCCGGGCACCTCTGCAAACGGCCCCCGCAGTCCGGTTCCGCCGCCCGGCAGGCGTGGGTTAGAGCAGGTGcc       >  1:87107/1‑148 (MQ=255)
                                                                                                                                           cgcgGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCTGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGaaa   <  7:352215/149‑1 (MQ=255)
                                                                                                                                            gcgggcggCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAg  <  7:304892/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                      
GTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAG  >  NC_002947/1184445‑1184733

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑