Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,757,688 T→A 26.0% intergenic (+31/+33) lpp → / ← ldtE murein lipoprotein/murein L,D‑transpeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,757,6880TA26.0% 38.3 / 11.6 23intergenic (+31/+33)lpp/ldtEmurein lipoprotein/murein L,D‑transpeptidase
Reads supporting (aligned to +/- strand):  ref base T (9/8);  new base A (3/3);  total (12/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

TTGACCAGCTGAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGAAACCGG  >  NC_000913/1757542‑1757831
                                                                                                                                                  |                                                                                                                                               
ttGACCAGCTGAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACAtt                                                                                                                                                <  1:469405/148‑1 (MQ=255)
 tGACCAGCTGAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATtg                                                                                                                                               <  2:208050/148‑1 (MQ=255)
       gCTGAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACAATGTGCGcc                                                                                                                                         <  2:666451/148‑1 (MQ=255)
        cTGAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACAATGTGCGCCa                                                                                                                                        <  2:589488/148‑1 (MQ=255)
              aaCGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCAttttttt                                                                                                                                 <  3:665071/149‑1 (MQ=255)
                       aaCGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCg                                                                                                                        >  3:407483/1‑149 (MQ=255)
                            aaTGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTAc                                                                                                                   >  3:430346/1‑149 (MQ=255)
                             ttGCGTTCCGACGTTCAGGCTGCTAAAGATGATGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGTAGTGAAAAATGGCGCACAATGTGCGCCa                                                                                                                                        >  3:357290/2‑127 (MQ=255)
                             aTGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACAATGTGCGCCa                                                                                                                                        <  4:357290/127‑1 (MQ=255)
                                                    tAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCAttttttt                                                                                                                                 >  1:28315/1‑111 (MQ=255)
                                                    tAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCAttttttt                                                                                                                                 <  2:28315/111‑1 (MQ=255)
                                                         aTGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTccc                                                                                      <  4:430346/149‑1 (MQ=255)
                                                                               gCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGctc                                                                                >  4:896490/1‑133 (MQ=255)
                                                                               gCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGctc                                                                                <  3:896490/133‑1 (MQ=255)
                                                                                   cTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCt                                                             >  2:894534/1‑148 (MQ=255)
                                                                                       acaacaTGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGAct                                                         >  1:593440/1‑148 (MQ=255)
                                                                                                       aTACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTg                                        >  3:831943/1‑149 (MQ=255)
                                                                                                                 tAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCg                               <  2:593440/148‑1 (MQ=255)
                                                                                                                                           gcgcACAATGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGAAAcc    >  4:13889/1‑149 (MQ=255)
                                                                                                                                           gcgcACAATGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGAAAc     >  2:589401/1‑148 (MQ=255)
                                                                                                                                             gcACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGAAACCgg  >  4:632591/1‑149 (MQ=255)
                                                                                                                                              cacaTTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGAAACCgg  <  1:894534/148‑1 (MQ=255)
                                                                                                                                                caTTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCggg                             <  2:895349/119‑1 (MQ=255)
                                                                                                                                                caTTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCggg                             >  1:895349/1‑119 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
TTGACCAGCTGAGCAACGACGTGAACGCAATGCGTTCCGACGTTCAGGCTGCTAAAGATGACGCAGCTCGTGCTAACCAGCGTCTGGACAACATGGCTACTAAATACCGCAAGTAATAGTACCTGTGAAGTGAAAAATGGCGCACATTGTGCGCCATTTTTTTTGTCTGCCGTTTACCGCTACTGCGTCACGCGTAACATATTCCCTTGCTCTGGTTCACCATTCTGCGCTGACTCTACTGAAGGCGCATTGCTGGCTGCGGGAGTTGCTCCACTGCTCACCGAAACCGG  >  NC_000913/1757542‑1757831

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