Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 19,634 (C)5→6 intergenic (+14/+177) nhaR → / ← insB1 DNA‑binding transcriptional activator NhaR/IS1 family protein InsB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091319,6341.C95.0% 57.1 / ‑2.8 20intergenic (+14/+177)nhaR/insB1DNA‑binding transcriptional activator NhaR/IS1 family protein InsB
Reads supporting (aligned to +/- strand):  ref base . (1/0);  new base C (7/12);  total (8/12)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

TTGATGATGCCGCTTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGTGCGTTAATCGGCAGCTCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACATTACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTAT  >  NC_000913/19385‑19814
                                                                                                                                                                                                                                                          |                                                                                                                                                                                     
gtgatgaTGCCGCGTTGTTGAAAGCTTTTGGTGAGCTGCAAAATGCAATCTCCGTTGACCCAACGTTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTATAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATGCTGCGCTTATTAGTCCAGCGGCGCGTTAATCGGCAACTcaccc                                                                                                                                                                                        <  2:152441/249‑1 (MQ=255)
               tgatgaAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATCCGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCTCCCCCCAAAGTTAAGGTgggg                                                                                                                                                                       >  2:94948/1‑251 (MQ=255)
                   aaaaGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCTACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTTGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAACCGGCACCTCCCCCCAAAGTTAAGGTGTGGgaga                                                                                                                                                                   <  2:163588/250‑1 (MQ=255)
                                    tgcacctttcactCTTCGTTGCCCCAACGCGTTATGCCTTTGACTTTGATGCCGGAAAAACTGGTGTAGTAATTGTACGCGTCGAGTATTTGATGGTATAGTACCATTGTATTTTTGCTGAGCGTATGATTCAGCAGCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCTCTTTTTAGGCCAGAGGCGCGTTAATCGGCACTTCCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACAt                                                                                                                                                   <  2:209940/238‑1 (MQ=255)
                                     gCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCTCCCCCCCAAGTTAAGGTGGGGGCGATAGATTAGTTGTACAta                                                                                                                                                  >  1:68260/1‑249 (MQ=255)
                                                       tGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCTCCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCgg                                                                                                                                <  1:94948/250‑1 (MQ=255)
                                                            cAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCTCCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCatt                                                                                                                          <  1:85700/251‑1 (MQ=255)
                                                                                 ttgaTGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCTCCCCCCAAAGTTAAGGTGGGGGATATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACAt                                                                                                       <  1:251904/246‑1 (MQ=255)
                                                                                     tGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCTCCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTc                                                                                                 >  1:266122/1‑251 (MQ=255)
                                                                                            aaaaCTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCAGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTCGGACAGCGGCGAGTTAATCTGCACCTCCCCCCAAAGTTAAGGTGGGGGAGATAGATTATTTGTACATAACCCCGATTTTGACTCGGCTCATTATTGTCCCGATTTAGACAATATTTCCATATGt                                                                                          >  2:201590/1‑251 (MQ=255)
                                                                                                    cGTAGAAATTGTGCGCGTCGAGAATTTTATGGATGAGTACCATGCTATTTTTGCTGAGCGTATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGTATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCTCCCCCCAAAGTTAAGGTTGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGATCATTATTTGCCTGCTTGAGACATTATTTCCATATGTACGCggg                                                                                   <  1:79690/250‑1 (MQ=255)
                                                                                                                          aaTGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGAATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATGCGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCTCCCCCCAAAGTTAAGGTAGGGGAGAGAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCAATATGGACGCGGGCGAATAAATAGAGGAATCTGAtt                                                            >  2:210806/1‑251 (MQ=255)
                                                                                                                                            cATGCTATTTTTTCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTAGTCTGCGTTTTTTAGTCCAGCGGCGCGTTAATCGGCACCTCCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATg                                          <  2:166177/251‑1 (MQ=255)
                                                                                                                                                        gCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCACTACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCTCCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACATAACAACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGCTGCTGAAAAGagt                               >  2:134627/1‑250 (MQ=255)
                                                                                                                                                                    atTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCTCCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCtttt                  >  1:12504/1‑248 (MQ=255)
                                                                                                                                                                               ggcggTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCTGCACCTCCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCttttaccatactgga       <  1:353701/251‑15 (MQ=255)
                                                                                                                                                                                gcggTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCTCCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCttttaccatactggat      <  1:201590/251‑16 (MQ=255)
                                                                                                                                                                                 cggTACATCGAACTTGCAATACGGATTAGTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCTCCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGGAAAGAGTAGTAATTGCttttaccatactggatt     <  2:119880/251‑17 (MQ=255)
                                                                                                                                                                                   gTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGGCCAGCGGCGCGTTAATCGGCACCTCCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGTGATGCTGAAAAGAGTAGTAATTGCttttaccatactggatttc   <  1:134206/251‑19 (MQ=255)
                                                                                                                                                                                    tACAGCGATTCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCTCCCCCCAAAGTTAAGGTGGGGGAGATAGCTTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGAGGATGCTGAAAAGAGGAGTAAtagatcgtaagatcgtcgtggagg  >  2:224728/1‑228 (MQ=255)
                                                                                                                                                                                    tACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCTCCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACATAACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTATTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAAt                         <  1:224728/228‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                     
TTGATGATGCCGCTTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGTGCGTTAATCGGCAGCTCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACATTACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTAT  >  NC_000913/19385‑19814

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