Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,182,864 A→T 16.0% intergenic (+25/+69) tadA → / ← mltF tRNA adenosine(34) deaminase TadA/membrane‑bound lytic murein transglycosylase F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,182,8640AT16.0% 59.6 / 5.8 25intergenic (+25/+69)tadA/mltFtRNA adenosine(34) deaminase TadA/membrane‑bound lytic murein transglycosylase F
Reads supporting (aligned to +/- strand):  ref base A (11/10);  new base T (2/2);  total (13/12)
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.72e-01

GCCAGGGGCAGTTCTTTGGCCAAGGCTTTCTCAACCACCGGGTGATGGTGGAGGGCGGGGTGCTGGCGCAAGAGTGCGGGCAGATCCTCAGCGACTTCTTCAAGGCCCGCCGGGCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGTTCTGCGCAAATCTTGTAGGCGCGGGTTTGCCCGCGAAGAGGTCCTGACAGGCTTACATCGGCGACGACTTATCCGCTGGCTTGTCCTTGTTGACCCCGGGCACATGCAGGTTGCCCTC  >  NC_002947/1182715‑1182998
                                                                                                                                                     |                                                                                                                                      
gCCAGGGGCAGTTCTTTGGCCAAGGCTTTCTCAACCACCGGGTGATGGTGGAGGGCGGGGTGCTGGCGCAAGAGTGCGGGCAGATCCTCAGCGACTTCTTCAAGGCCCGCCGGGCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCaa                                                                                                                                        <  1:72506/150‑1 (MQ=255)
 ccAGGGGCAGTTCTTTGGCCAAGGCTTTCTCAACCACCGGGTGATGGTGGAGGGCGGGGTGCTGGCGCAAGAGTGCGGGCAGATCCTCAGCGACTTCTTCAAGGCCCGCCGGGCCAAGGCTTGAGCTGTACCGGCCTCTTCGCGGGCAAg                                                                                                                                       >  2:296088/1‑150 (MQ=255)
                    cAAGGCTTTCTCAACCACCGGGTGATGGTGGAGGGCGGGGTGCTGGCGCAAGAGTGCGGGCAGATCCTCAGCGACTTCTTCAAGGCCCGCCGGGCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGTTCTg                                                                                                                    >  1:375033/1‑150 (MQ=255)
                     aaGGCTTTCTCAACCACCGGGTGATGGTGGAGGGCGGGGTGCTGGCGCAAGAGTGCGGGCAGATCCTCAGCGACTTCTTCAAGGCCCGCCGGGCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGTTCTgc                                                                                                                   <  2:492472/150‑1 (MQ=255)
                                        ggTGATGGTGGAGGGCGGGGTGCTGGCGCAAGAGTGCGGGCAGATCCTCAGCGACTTCTTCAAGGCCCGCCGGGCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGTTCTGCGCAAATCTTGTAGGCGCgg                                                                                                >  1:590119/1‑150 (MQ=255)
                                                tggAGGGCGGGGTGCTGGCGCAAGAGTGCGGGCAGATCCTCAGCGACTTCTTCAAGGCCCGCCGGGCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGTTCTGCGCAAATCTTGTAGGCGCGGGTTTGccc                                                                                        >  2:123264/1‑150 (MQ=255)
                                                   aGGGCGGGGTGCTGGCGCAAGAGTGCGGGCAGATCCTCAGCGACTTCTTCAAGGCCCGCCGGGCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCAAACCCGCGCCTACAAGATTTGCGCAGAACCTGTAGGAGCGGGCTTGCCcgcg                                                                                     <  1:601391/150‑1 (MQ=21)
                                                       cGGGGTGCTGGCGCAAGAGTGCGGGCAGATCCTCAGCGACTTCTTCAAGGCCCGCCGGGCCAAGACTTGAGCTGCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGTTCTGCGCAAATCTTGTAGGCGCGGGTTTGCCCGCGAaga                                                                                 <  2:333733/150‑1 (MQ=255)
                                                                  cgcAAGAGTGCGGGCAGATCCTCAGCGACTTCTTCAAGGCCCGCCGGGCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGTTCTGCGCAAATCTTGTAGGCGCGGGTTTGCCCGCGAAGAGGTCCTGACAg                                                                      <  2:516737/150‑1 (MQ=255)
                                                                   gcAAGAGTGCGGGCAGATCCTCAGCGACTTCTTCAAGGCCCGCCGGGCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTAAAGGTTGTGCGCAAATCTTGTAGGCGCGGGTTTGCCCGCGAAGAGATCCTGCCAgg                                                                     >  1:183338/1‑150 (MQ=255)
                                                                              ggCAGATCCTCAGCGACTTCTTCAAGGCCCGCCGGGCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGTTCTGCGCAAATCTTGTAGGCGCGGGTTTGCCCGCGAAGAGGTCCTGACAGGCTTACATCGGc                                                          <  1:111470/150‑1 (MQ=255)
                                                                              ggCAGATCCTCAGCGACTTCTTCAAGGCCCGCCGGGCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGTTCTGCGCAAATCTTGTAGGCGCGGGTTTGCCCGCGAAGAGGTCCTGACAGGCTTACATCGGc                                                          <  1:111435/150‑1 (MQ=255)
                                                                                       tCAGCGACTTCTTCAAGGCCCGCCGGGCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCAAGCCCGCTc                                                                                                                                <  2:363304/71‑1 (MQ=255)
                                                                                       tCAGCGACTTCTTCAAGGCCCGCCGGGCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCAAGCCCGCTc                                                                                                                                >  1:363304/1‑71 (MQ=255)
                                                                                       tCAGCGACTTCTTCAAGGCCCGCCGGGCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGTTCTGCGCAAATCTTGTAGGCGCGGGTTTGCCCGCGAAGAGGTCCTGACAGGCTTACATCGGCGACGACTTa                                                 >  1:491656/1‑150 (MQ=255)
                                                                                                tcttcAAGGCCCGCCGGGCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGTTCTGCGCAAATCTTGTAGGCGCGGGTTTGCCCGCGAAGAGGTCCTGACAGACTTACATCGGCGACGACTTATCCGCTGGc                                        <  2:375033/150‑1 (MQ=255)
                                                                                                 cttcAAGGCCCGCCGGGCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGTTCTGCGCAAATCTTGTAGGCGCGGGTTTGCCCGCGAagag                                                                                >  1:448830/1‑109 (MQ=255)
                                                                                                 cttcAAGGCCCGCCGGGCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGTTCTGCGCAAATCTTGTAGGCGCGGGTTTGCCCGCGAagag                                                                                <  2:448830/109‑1 (MQ=255)
                                                                                                  ttcAAGGCCCGCCGGGCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCTTGCCCGcgaa                                                                                                                               >  1:373701/1‑58 (MQ=21)
                                                                                                  ttcAAGGCCCGCCGGGCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCTTGCCCGcgaa                                                                                                                               <  2:373701/61‑4 (MQ=21)
                                                                                                      attcccGCCGGGCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGTTCTGCGCAAATATTGTAGGCGCGGGTTTGCCCGCGAAGAGGTCCTGACAGGCTTACATCGGCGACGACTTATCCGCTGGCTTGTcc                                  >  1:547523/4‑150 (MQ=255)
                                                                                                                 gCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGTTCTGCGCAAATCTTGTAGGCGCGGGTTTGCCCGCGAAGAGGTCCTGACAGGCTTACATCGGCGACGACTTATCCGCTGGCTTGTCCTTGTTGAcccc                       <  2:590119/150‑1 (MQ=255)
                                                                                                                                   ccGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGTTCTGCGCAAATCTTGTAGGCGCGGGTTTGCCCGCGAAGAGGTCCTGACAGGCTTACATCGGCGACGACTTATCCGCTGGCTTGTCCTTGTTGACCCCGGGCACATGCAGGTTGcc     >  2:548282/1‑150 (MQ=255)
                                                                                                                                      gCCTCTTCGCGGGCAAGCCCGCTCCTACAGGTTCTGCGCAAATCTTGTAGGCGCGGGTTTGCCCGCGAAGAGGTCCTGACAGGCTTACATCGGCGACGACTTATCCGCTGGCTTGTCCTTGTTGACCCCGGGCACATGCAGGTTGCCCTc  >  2:608128/1‑150 (MQ=255)
                                                                                                                                        gtaggagcgGGCTTGCCCGCTCCTACAGGTTCTGCGCAAATCTTGTAGGCGCGGGTTTGCCCGCGAAGAGGTCCTGACAGGCTTACATCGGCGACGACTTATCCGCTGGCTTGTCCTTGTTGACCCCGGGCACATGCAGGTTGCCCTc  >  1:492915/7‑148 (MQ=255)
                                                                                                                                        gtaggagcgGGCTTGCCCGCTCCTACAGGTTCTGCGCAAATCTTGTAGGCGCGGGTTTGCCCGCGAAGAGGTCCTGACAGGCTTACATCGGCGACGACTTATCCGCTGGCTTGTCCTTGTTGACCCCGGGCACATGCAGGTTGCCCTc  <  2:492915/142‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                      
GCCAGGGGCAGTTCTTTGGCCAAGGCTTTCTCAACCACCGGGTGATGGTGGAGGGCGGGGTGCTGGCGCAAGAGTGCGGGCAGATCCTCAGCGACTTCTTCAAGGCCCGCCGGGCCAAGGCTTGAGCTGCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACAGGTTCTGCGCAAATCTTGTAGGCGCGGGTTTGCCCGCGAAGAGGTCCTGACAGGCTTACATCGGCGACGACTTATCCGCTGGCTTGTCCTTGTTGACCCCGGGCACATGCAGGTTGCCCTC  >  NC_002947/1182715‑1182998

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