Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,184,588 G→C 17.1% intergenic (‑204/‑208) 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,5880GC17.1% 93.3 / 14.8 41intergenic (‑204/‑208)mltF/purLmembrane‑bound lytic murein transglycosylase F/phosphoribosylformylglycinamidine synthase
Reads supporting (aligned to +/- strand):  ref base G (20/14);  new base C (4/3);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.04e-01

GTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCA  >  NC_002947/1184445‑1184736
                                                                                                                                               |                                                                                                                                                      
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                                                                      <  1:237503/139‑1 (MQ=255)
                                                                                     aGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACctc                                                                      >  2:237503/1‑139 (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                                                          <  3:99319/110‑1 (MQ=255)
                                                                                                                              agaTCGAGCGCCGCCCGCGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGAcc                                                          >  4:99319/1‑110 (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)
                                                                                                                                               ggcggcGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCa  <  4:318916/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                      
GTAGCGGGTTGTGTGGTTGGGGGAGGGGGGATTGCATCAGGGCTGCCAGGTGGCCGGCTTGAGATATTCAGCGCCTGTGAGATCGAGCGCCGCCAGGTGTCCGCCTTAAGATATTCAGCGCCTATGAGATCGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCAGAATGTGTCGAGCCGGACACCTGCG‑‑GCGCTCGATCTCGTCTCAAGCCGGACACCTCTGAAAACTGACCCCGAAGTCCGGTTCCGCCGCCCGGCAGGCGTGGCTTAGAGCAGGTGCCGAAAGCCA  >  NC_002947/1184445‑1184736

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: ‑