Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 19,627 G→C intergenic (+7/+184) 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,6270GC95.5% 47.7 / ‑6.7 22intergenic (+7/+184)nhaR/insB1DNA‑binding transcriptional activator NhaR/IS1 family protein InsB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (9/12);  minor base A (0/1);  total (9/13)
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.97e-01

AGTTTGATGATGCCGCTTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGTGCGTTAATCGGCAGCTCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACATTACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTAT  >  NC_000913/19382‑19814
                                                                                                                                                                                                                                                     |                                                                                                                                                                                             
aGTTTGATGATGCCGCGTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGCGCGTTAATCGGCACCTccc                                                                                                                                                                                          >  1:253904/1‑251 (MQ=255)
   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)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                             
AGTTTGATGATGCCGCTTTGATGAAAGCTTTTGGTGCGATGCACAATGCAATCTTCGTTGCCCCAACGCTTTATGCATATGACTTTTATGCCGATAAAACTGTCGTAGAAATTGGTCGCGTCGAGAATGTGATGGAAGAGTACCATGCTATTTTTGCTGAGCGGATGATTCAGCACCCGGCGGTACAGCGAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGTGCGTTAATCGGCAGCTCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACATTACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTAT  >  NC_000913/19382‑19814

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