Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,824,478 (C)6→7 intergenic (‑132/+13) pncC ← / ← mltB nicotinamide‑nucleotide amidohydrolase, NMN amidohydrolase/membrane‑bound lytic murein transglycosylase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,824,4721.C93.3% 44.7 / ‑2.7 15intergenic (‑126/+19)pncC/mltBnicotinamide‑nucleotide amidohydrolase, NMN amidohydrolase/membrane‑bound lytic murein transglycosylase B
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (6/8);  total (7/8)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

CTGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTAC  >  NC_000913/2824336‑2824599
                                                                                                                                         |                                                                                                                                
cTGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGcc                                                                                                                               >  2:474748/1‑139 (MQ=255)
    cagtcaTCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGcccccc                                                                                                                           <  1:229169/139‑1 (MQ=255)
      gtcaTCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGccccccct                                                                                                                         <  1:474748/139‑3 (MQ=255)
           cacaACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGt                                                                                                                     >  2:205645/1‑139 (MQ=255)
            acaACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTa                                                                                                                    <  1:277743/139‑1 (MQ=255)
                 ctctGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAg                                                                                                               >  2:615917/1‑139 (MQ=255)
                 ctctGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGcccccc                                                                                                                           >  1:427726/1‑126 (MQ=255)
                 ctctGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGcccccc                                                                                                                           >  1:447636/1‑126 (MQ=255)
                 ctctGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGcccccc                                                                                                                           <  2:447636/126‑1 (MQ=255)
                 ctctGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGcccccc                                                                                                                           <  2:427726/126‑1 (MQ=255)
                                                         tCGTCAGGAAGCTGGTGCTGGCCTTGCAATTTTACGAGGAGGATTCAGAAACAAGCTGATTAGCCAGAGGGAAGCTCACGCTCCCCCCCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                                                                      >  2:419058/1‑139 (MQ=255)
                                                               ggAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCg                                                                <  1:742800/139‑1 (MQ=255)
                                                                        tgtTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGt                                                       >  2:653082/1‑139 (MQ=255)
                                                                          ttGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGcccc                                                                                                                             >  2:718213/1‑67 (MQ=255)
                                                                          ttGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGcccc                                                                                                                             <  1:718213/67‑1 (MQ=255)
                                                                                  gCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGTTACTTACTCgcgccagcgcc                                                                                         >  2:142769/1‑95 (MQ=255)
                                                                                  gCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGTTACTTACTCgcgccagcgcc                                                                                         <  1:142769/95‑1 (MQ=255)
                                                                                     aTTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCtgtg                                          <  1:713897/139‑1 (MQ=255)
                                                                                                   aTTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGa                            <  1:414363/139‑1 (MQ=255)
                                                                                                                        ccAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGAc       >  1:815787/1‑139 (MQ=255)
                                                                                                                        ccAGAGGGAAGCTCACGCCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGAc       <  1:653082/139‑1 (MQ=255)
                                                                                                                             gggAAGCTCACGCCCCCCCTCTTGTAAATAGTTACTTACTCGCGCCAGCGCCACGGCTTGTCCGAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTAc  <  1:615917/139‑1 (MQ=255)
                                                                                                                                 aGCTCACGCCCCCCCTCTTGTAAATAGTTACTTACTCgcgc                                                                                                >  2:157305/1‑41 (MQ=255)
                                                                                                                                 aGCTCACGCCCCCCCTCTTGTAAATAGTTACTTACTCgcgc                                                                                                <  1:157305/41‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                
CTGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCCCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTAC  >  NC_000913/2824336‑2824599

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