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

CTGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACG‑C‑CCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTAC  >  NC_000913/2824336‑2824599
                                                                                                                                         |                                                                                                                                
cTGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACG‑c‑c                                                                                                                               >  2:474748/1‑139 (MQ=255)
    cagtcaTCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACG‑c‑ccccc                                                                                                                           <  1:229169/139‑1 (MQ=255)
      gtcaTCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACG‑c‑cccccct                                                                                                                         <  1:474748/139‑3 (MQ=255)
           cacaACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCC‑CCCCCTCTTGt                                                                                                                     >  2:205645/1‑139 (MQ=255)
            acaACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCC‑CCCCCTCTTGTa                                                                                                                    <  1:277743/139‑1 (MQ=255)
                 ctctGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCC‑CCCCCTCTTGTAAATAg                                                                                                               >  2:615917/1‑139 (MQ=255)
                 ctctGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACG‑c‑ccccc                                                                                                                           >  1:427726/1‑126 (MQ=255)
                 ctctGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACG‑c‑ccccc                                                                                                                           >  1:447636/1‑126 (MQ=255)
                 ctctGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACG‑c‑ccccc                                                                                                                           <  2:447636/126‑1 (MQ=255)
                 ctctGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACG‑c‑ccccc                                                                                                                           <  2:427726/126‑1 (MQ=255)
                                                         tCGTCAGGAAGCTGGTGCTGGCCTTGCAATTTTACGAGGAGGATTCAGAAACAAGCTGATTAGCCAGAGGGAAGCTCACG‑CTCCCCCCCTTGTAAATAGTTACT‑TACTCGCGCCAGCGCCACGGCTTGTCCTAACTGcc                                                                      >  2:419058/1‑139 (MQ=255)
                                                               ggAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCC‑CCCCCTCTTGTAAATAGTTACT‑TACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCg                                                                <  1:742800/139‑1 (MQ=255)
                                                                        tgtTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCC‑CCCCCTCTTGTAAATAGTTACT‑TACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGt                                                       >  2:653082/1‑139 (MQ=255)
                                                                          ttGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACG‑c‑ccc                                                                                                                             >  2:718213/1‑67 (MQ=255)
                                                                          ttGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACG‑c‑ccc                                                                                                                             <  1:718213/67‑1 (MQ=255)
                                                                                  gCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCC‑CCCCCTCTTGTAAATAGTTACT‑TACTCgcgccagcgcc                                                                                         >  2:142769/1‑95 (MQ=255)
                                                                                  gCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCC‑CCCCCTCTTGTAAATAGTTACT‑TACTCgcgccagcgcc                                                                                         <  1:142769/95‑1 (MQ=255)
                                                                                     aTTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCC‑CCCCCTCTTGTAAATAGTTACT‑TACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCtgtg                                          <  1:713897/139‑1 (MQ=255)
                                                                                                   aTTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACGCC‑CCCCCTCTTGTAAATAGTTACT‑TACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGa                            <  1:414363/139‑1 (MQ=255)
                                                                                                                        ccAGAGGGAAGCTCACGCC‑CCCCCTCTTGTAAATAGTTACT‑TACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGAc       >  1:815787/1‑139 (MQ=255)
                                                                                                                        ccAGAGGGAAGCTCACGCC‑CCCCCTCTTGTAAATAGTTACT‑TACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGAc       <  1:653082/139‑1 (MQ=255)
                                                                                                                             gggAAGCTCACGCC‑CCCCCTCTTGTAAATAGTTACT‑TACTCGCGCCAGCGCCACGGCTTGTCCGAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTAc  <  1:615917/139‑1 (MQ=255)
                                                                                                                                 aGCTCACGCC‑CCCCCTCTTGTAAATAGTTACT‑TACTCgcgc                                                                                                >  2:157305/1‑41 (MQ=255)
                                                                                                                                 aGCTCACGCC‑CCCCCTCTTGTAAATAGTTACT‑TACTCgcgc                                                                                                <  1:157305/41‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                
CTGTCAGTCATCACAACCTCTGTCAGTCAAAAATTTATGCCGGACAAGATAGCACTTTCGTCAGGAAGATGGTGTTGGCGTTGCAATTTTACGAGGAGGATTCAGAAAAAAGCTGATTAGCCAGAGGGAAGCTCACG‑C‑CCCCCTCTTGTAAATAGCTACTGTACTCGCGCCAGCGCCACGGCTTGTCCTAACTGCCAGACCGCCATTGCGTAATGGGTGCTGTGGTTGTAACGGGTGATGGTGTAGAAGTTCGGCAGACCGTAC  >  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: ‑