Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,405,758 G→A intergenic (‑7/‑221) mcaS ← / → smrA small regulatory RNA McaS/DNA endonuclease SmrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,405,7580GA95.5% 62.8 / ‑5.1 22intergenic (‑7/‑221)mcaS/smrAsmall regulatory RNA McaS/DNA endonuclease SmrA
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (8/13);  total (9/13)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGACCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGACTCTACAGTACACACAGCAGTGCATCCGCGTCTTAAATCCGGCATTGTCTCCTCTGCGCCGGTGACTGTGTTATAACGGTTTTAACCTTTAGTTGCCAATTTTCTTCAGTGACAGATTTCACGAAAATATAAGTTACATTATTAATATCGTGAATGAATAATCATGCATAAGTATTTTGCTTAAAATATCGGCAATATTTGGAACTTATTACTGGAAATTTGGGTAATACGTTGTTGGACCGACCCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTG  >  NC_000913/1405535‑1406002
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
aaTTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGACCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGACTCTACAGTACACACAGCAGTGCATCCGCGTCTTAAATCCGGCATTGTCTCCTCTGCGCCGGTGACTGTATTATAACGGTTTTAACCTTTAGTTGcc                                                                                                                                                                                                                           <  1:539376/251‑1 (MQ=255)
    ttACCGAAAGCGGCCTATTATTTCATTCAGCGGCGTGATCGTCGCCAGACCCCACTGACCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAATGTTAAAACTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGACTCTACAGTTCACACAGCAGTGCATCCGCGTCTTAAACCCGGCATTGTCTCCTCTGCGCCGGTGACTGTATTATAACGGTTTTAACCTTTAGTTGCCAAt                                                                                                                                                                                                                        <  1:27308/250‑1 (MQ=255)
               ggCCTATTATTGGATTCAGCGGCGTGATAGTCGCCAGACACCACTGACCGCGTCATTAATCACGCAATACCGGCGAGAATGCGGCTATCTGCAAATTTAAAACTGCATAAAAAAAAAGAGTCTGTCGACAACCGCCAGACTCTACAGTACACACAGCAGTGCATCCGCGTCTTAAATCCGGCATTGTCTCCTCTGCGCCGGTGACTGTATTATAACGGTTTTAGCCTTTAGTTGCCAATTTTCTTCAGTGa                                                                                                                                                                                                            <  2:461562/251‑1 (MQ=255)
                 ccTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGACCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGACTCTACAGTACACACAGCAGTGCATCCGCGTCTTAAATCCGGCATTGTCTCCTCTGCGCCGGTGACTGTATTATAACGGTTTTAACCTTTAGTTGCCAATTTTCTTCAGTGACa                                                                                                                                                                                                          <  1:164847/251‑1 (MQ=255)
                        tttGATTCAGCGGCGTGATAGTCGCCAGACACCACTGACCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGACTCTACAGTACACACAGCAGTGCATCCGCGTCTTAAATCCGGCATTGTCTCCTCTGCGCCGGTGACTGTATTATAACGGTTTTAACCTTTAGTTGCCAATTTTCTTCAGTGACAGATTTc                                                                                                                                                                                                    >  1:173763/1‑250 (MQ=255)
                         ttGATTCAGCGGCGTGATAGTCGCCAGACACCACTGACCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGACTCTACAGTACACACAGCAGTGCATCCGCGTCTTAAATCCGGCATTGTCTCCTCTGCGCCGGTGACTGTATTATAACGGTTTTAACCTTTAGTTGCCAATTTTCTTCAGTGACAGATTTCa                                                                                                                                                                                                   <  2:173763/250‑1 (MQ=255)
                                aGCGGCGTGATAGTCGCCAGACACCACTGACCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGACTCTACAGTACACACAGCAGTGCATCCGCGTCTTAAATCCGGCATTGTCTCCTCTGCGCCGGTGACTGTATTATAACGGTTTTAACCTTTAGTTGCCAATTTTCTTCAGTGACAGATTTCACGAAAata                                                                                                                                                                                           <  1:16621/251‑1 (MQ=255)
                                              cGCCAGACACCACTGACCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGACTCTACAGTACACACAGCAGTGCATCCGCGTCTTAAATCCGGCATTGTCTCCTCTGCGCCGGTGACTGTATTATAACGGTTTTAACCTTTAGTTGCCAATTTTCTTCCGTGACAGATTTCACGAAAATATAAGTTACattatt                                                                                                                                                                             >  1:219323/1‑251 (MQ=255)
                                                 cAGACACCACTGACCGCGTCATTAATCACGCACTTCCGGCGAGAATGCGGCTATCTGCAAAGTTACAACTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGCCTCTACAGTACACACAGCAGTGCATCCGCGTCTTAAATCCGGCATTGTCTCATCTGCGCCGGTGACTGTATTATAACGGTTTTAACCTTTAGTTGCCATTTTCCTTCACTGCCAGATTTCACGAAAATTTAATTTACAtaattaat                                                                                                                                                                          >  1:532634/1‑251 (MQ=255)
                                                                                aaTTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGACAGCCGCCAGACTCTACAGTACACACAGCAGTGCATCCGCGTCTTAAATCCGGCATTGTCTCCTCTGCGCCGGTGACTGTATTATAACGGTTTTAACCTTTAGTTGCCAATTTTCTTCAGTGACAGATTTCAAGAAAATATAAGTTACATTATTAATATCGTGAATGAATAATCATGCATAAGTAttt                                                                                                                                           <  1:360070/251‑1 (MQ=255)
                                                                                                    tATCTGCACAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGACTCTACAGTACACACAGCAGTGCATCCGCGTCTTAAATCCGGCATTGTCTCCTCTGCGCCGGTGACTGTATTATAACGGTTTTAACCTTTAGTTGCCAATTTTCTTCAGTGACAGATTTCACGAAAATATAAGTTACATTATTAATATCGTGAATGAATAATCATGCATAAGTATTTTGCTTAAAATATCGGCAat                                                                                                                        >  1:200993/1‑250 (MQ=255)
                                                                                                                   aaCTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGACTCTACAGTACACACAGCAGTGCATCCGCGTCTTAAATCCGGCATTGTCTCCTCTGCGCCGGTGACTGTATTATAACGGTTTTAACCTTTAGTTGCCAATTTTCTTCAGTGACAGATTTCACGAAAATATAAGTTACATTATTACTATCGTGAATGAATAATCATGCATAAGTATTTTGCTTAAAATATCGGCAATATTTGGAACTTATTAc                                                                                                        >  2:567583/1‑251 (MQ=255)
                                                                                                                    aCTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGACTCTACAGTACACACAGCAGTGCATCCGCGTCTTAAATCCGGCATTGTCTCCTCTGCGCCGGTGACTGTATTATAACGGTTTTAACCTTTAGTTGCCAATTTTCTTCAGTGACAGATTTCACGAAAATATAAGTTACATTATTAATATCGTGAATGAATAATCATGCATAAGTATTTTGCTTAAAATATCGGCAATATTTGGAACTTATTACt                                                                                                       <  2:152228/251‑1 (MQ=255)
                                                                                                                                        cTGTCGACATCCGCCAGACTCTACAGTACACACAGCAGTGCATCCGCGTCTTAAATCCGGCATTGTCTCCTCTGCGCCGGTGACTGTATTATAACGGTTTTAACCTTTAGTTGCCAATTTTCTTCAGTGACAGATTTCACGAAAATATAAGTTACATTATTAATATCGTGAATGAATAATCATGCATAAGTATTTTGCTTAAAATATCGGCAATATTTGGAACTTATTACTGGACATTTGGGTAATACgtt                                                                                   >  1:117368/1‑251 (MQ=255)
                                                                                                                                                     ccAGACTCTACAGTACACACAGCAGTGCAGCCGCGTCTTAAATCCGGCATTGTCTCCTCTGCGCCGGTGACTGTATTATAACGGTTTTAACCTTTAGTTGCCAATTTTCTTCAGTGACAGATTTCACGAAAATATAAGTTACATTATTAATATCGTGAATGAATAATCATGCATAAGTATTTTGCTTAAAATATCGGCAATATTTGGAACTTATTACTGGAAATTTGGGTAATACGTTGTTGGACCGAccc                                                                      <  2:200993/251‑1 (MQ=255)
                                                                                                                                                                              tGCATCCGCCTCCTAAATCCGGACTAGTCTCCTCTGAGCCGGCGATTGTGTTATAACGGTTTTAACCTTTAGCTGCCAATTTCATTCAGTGTCCGATTTACCGAAAATATCAGTTACATTATTAATATCGTGAATGAATAATCATGCATAAGTATTTTGCTTATAATATCGGCAATATTTGGAACTTATTACTGGAACTTTGGGTAATACGTTGTTGGACCGACCCGGTCAGGTTATCATATcgctgtctt                                             >  2:598303/1‑245 (MQ=255)
                                                                                                                                                                                      cggcTTAAATCCGGCATTGTCGACTCTGCACCTGTGACTGTATTAAAACGGTTTTAACCTGTAGTTGCCAATTTTCTTCAGTGACATATTTCACGAAAATAGAAGTTTCATTATTAATATCGTGAATGAATAATCATGCATAAGTATATTTCTTAAAATATCGGCAGTATTTGGAACTTATTACTGGAAATTTGGGTAATACTTTGTTGGACCGCCCCGGTCTGGGTATTATATCGCGATCTTAATTgt                                       <  2:356414/246‑2 (MQ=255)
                                                                                                                                                                                              aTCCGGCATTGTCTCCTCTGCGCCGGTGACTGTATTATAACGGTTTTAACCTTTAGTTGCCAATTTTCTTCAGTGACAGATTTCACGAAAATATAAGTTACAGTATTAATATCGTGAATGAATAATCATGCATAAGTATTTTGCTTAAAATATCGGCAATATTTGGAACTTATTACTGGAAATTTGGTTAATACGTTGTTGGACCGACCCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGt                             <  2:607954/251‑1 (MQ=255)
                                                                                                                                                                                                                    ccGGTGACTGTATTATAACGGTTTTAACCTTTAGTTGCCAATTTTCTTCAGTGACAGATTTCACGAAAATATAAGTTACATTATTAATATCGTGAATGAATAATCATGCATAAGTATTTTGCTTAAAATATCGGCAATATTTg                                                                                                                   >  2:448730/1‑143 (MQ=255)
                                                                                                                                                                                                                    ccGGTGACTGTATTATAACGGTTTTAACCTTTAGTTGCCAATTTTCTTCAGTGACAGATTTCACGAAAATATAAGTTACATTATTAATATCGTGAATGAATAATCATGCATAAGTATTTTGCTTAAAATATCGGCAATATTTg                                                                                                                   <  1:448730/143‑1 (MQ=255)
                                                                                                                                                                                                                        tGACTGTATTATAACGGTTTTAACCTTTAGTTGCCAATTTTCTTCAGTGACAGATTTCACGAAAATATAAGTTACATTATTAATATCGTGAATGAATAATCATGCATAAGTATTTTGCTTAAAATATCGGCAATATTTGGAACTTATTACTGGAAATTTGGGTAATACGTTGTTGGACCGACCCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCt   >  2:332462/1‑251 (MQ=255)
                                                                                                                                                                                                                         gACTGTATTATAACGGTTTTAACCTTTAGTTGCCAATTTTCTTCAGTGACAGATTTCACGAAAATATAAGTTACATTATTAATATCGTGAATGAATAATCATGCATAAGTATTTTGCTTAAAATATCGGCAATATTTGGAACTTATTACTGGAAATTTGGGTAATACGTTGTTGGACCGACCCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTg  <  2:117368/251‑1 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
AATTTTACCGAAAGCGGCCTATTATTTGATTCAGCGGCGTGATAGTCGCCAGACACCACTGACCGCGTCATTAATCACGCAATTCCGGCGAGAATGCGGCTATCTGCAAAGTTAAAACTGCATAAAAAAATAGAGTCTGTCGACATCCGCCAGACTCTACAGTACACACAGCAGTGCATCCGCGTCTTAAATCCGGCATTGTCTCCTCTGCGCCGGTGACTGTGTTATAACGGTTTTAACCTTTAGTTGCCAATTTTCTTCAGTGACAGATTTCACGAAAATATAAGTTACATTATTAATATCGTGAATGAATAATCATGCATAAGTATTTTGCTTAAAATATCGGCAATATTTGGAACTTATTACTGGAAATTTGGGTAATACGTTGTTGGACCGACCCGGTCTGGTTATCATATCGCGCTCTTAATTGCGGGAGGATGTAACATGAACCTTGACGACAAATCGCTG  >  NC_000913/1405535‑1406002

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