Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 89,570 A→C intergenic (+538/‑64) cra → / → mraZ DNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091389,5700AC96.3% 86.5 / NA 27intergenic (+538/‑64)cra/mraZDNA‑binding transcriptional dual regulator Cra/DNA‑binding transcriptional repressor MraZ
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (19/7);  total (20/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTTCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCAGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTATCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATC  >  NC_000913/89327‑89801
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       
gTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGt                                                                                                                                                                                                                                  <  2:40408/251‑1 (MQ=255)
                          aTGCAGCTCGCGTCGCAAACTGACACTTTATCTTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGGAGTTGTgg                                                                                                                                                                                                        >  2:199025/1‑251 (MQ=255)
                                cTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCaa                                                                                                                                                                                                   >  2:16227/1‑250 (MQ=255)
                                 tCGCGTCGCAAACTGACACTTTATATTTGCTGTGTAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCaa                                                                                                                                                                                                   >  1:148837/1‑249 (MQ=255)
                                     gTCGCACACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGGGGGGCAAAGTggg                                                                                                                                                                                             >  2:241634/1‑251 (MQ=255)
                                        gCAAACTGAAACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAAt                                                                                                                                                                                          >  2:266521/1‑251 (MQ=255)
                                                         atTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGc                                                                                                                                                                          <  1:163402/250‑1 (MQ=255)
                                                                      aaTAGTGAGTCATTTTAAAACGGTGATGACGTAGAGGGATTTTTTCTTACAGATATTCATAACGCTAATTTGCTGCGCAAGTTGGACGACAAATACACACCGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGATTTTCCTCCGCTCCGTAAACTCACTTCAGTGGGACATTATGGGGCAAAATGGGAATAAGGGGTGAGGCTGGCCTGTTCCCGGGAGc                                                                                                                                                            >  1:19687/1‑251 (MQ=255)
                                                                       aTAGTGAGTCATTTTACAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGCAAAAGAGGGTATAAGGGGTGAGGCTGGCATGTTCCGGGGAGc                                                                                                                                                            >  1:192941/1‑250 (MQ=255)
                                                                        tAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCa                                                                                                                                                           >  2:191793/1‑250 (MQ=255)
                                                                                         aCGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGAc                                                                                                                                         >  1:47734/1‑251 (MQ=255)
                                                                                                       gAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCtt                                                                                                                           >  1:231731/1‑251 (MQ=255)
                                                                                                                  tCTTACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAATGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCt                                                                                                                 >  2:209220/1‑250 (MQ=255)
                                                                                                                     tACAGCTATTCATAACGTTAATTTGCTGCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGTCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTa                                                                                                                <  2:218029/248‑1 (MQ=255)
                                                                                                                                                      gTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCg                                                                            >  1:186496/1‑251 (MQ=255)
                                                                                                                                                         ggACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTGCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCAGGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCTGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGt                                                                          <  2:10079/250‑1 (MQ=255)
                                                                                                                                                             gTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGTTCAAAt                                                                     <  1:219046/251‑1 (MQ=255)
                                                                                                                                                              tAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCCCGGTAAACTCCTTTCAGTTGGAAATTGTGGGGAAAAGTGGGAATAAGTGGTTAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGTGCGCTTATCAGTGCCTACCCGTTACCGGGCACCGCTGCTTGAGAACGCTGCAGGTccactg                                                                    >  2:45856/1‑246 (MQ=255)
                                                                                                                                                                           gCTGATATTCGCCGTAAACATCGGGTTTATTACCTCGGTCTGCCTTGTGACTGGCTTGACAAGCTTTTACTCAGCCCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGCGGCTGGAATGTTCCGGGGAGCAACGTTAGTCAATCTCGCCAGCAAAGGGCGCTTATCCGTGCCTACCCGTTCCCGTGAACCGCTGCGTGAGAACGCTGCCGGTCAAATGGTTCGCACCAttt                                                       >  2:227655/1‑250 (MQ=255)
                                                                                                                                                                             tGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCAGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTCTCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGAc                                                     >  1:128770/1‑251 (MQ=255)
                                                                                                                                                                                         tAAACATCGGGTTTTTTACCTCGGTATGCCTTGCGACTGGCTTGAAAAGCTTTTCCTCCCCTCCGTAAACTCCTTTAAGTGGGTAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGGCAAACTCGACATCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTCATTGAGAACGATGCCCGTCAAATGGTTTGCACCATTGACATTCATCAcacg                                         >  2:89596/1‑248 (MQ=255)
                                                                                                                                                                                               tCGGGTTTTTTACCTCGGTATGCTTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGTGAAGGGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCATCGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGATAGCTCTGTCTGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCct                                    <  1:116527/250‑1 (MQ=255)
                                                                                                                                                                                                 gggTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCctgct                                 <  2:75625/251‑1 (MQ=255)
                                                                                                                                                                                                                         gtgACAGGCTTGAAAAGTTTTTCATCCGCTCCGTAAAATCCTTTCAGAGGGAAATTGTGGGGCAAAGTGGGCATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAAAGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACCGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCCGCTCCTTTACCCCCTACCTGCATgcg          >  2:261918/1‑248 (MQ=255)
                                                                                                                                                                                                                             cTGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAAtt     >  1:186251/1‑251 (MQ=255)
                                                                                                                                                                                                                              tGGCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAAtt     >  1:252392/1‑250 (MQ=255)
                                                                                                                                                                                                                                gCTTGACAAGCTTTTCCTCCGCTCCGTAAACTCCTTTCAGTGGGACATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTACCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTCATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATc  >  1:122865/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                       
GTTAACCGGGGAAGATATGTCCTAAAATGCCGCTCGCGTCGCAAACTGACACTTTATATTTGCTGTGGAAAATAGTGAGTCATTTTAAAACGGTGATGACGATGAGGGATTTTTTCTTACAGCTATTCATAACGTTAATTTGCTTCGCACGTTGGACGTAAAATAAACAACGCTGATATTAGCCGTAAACATCGGGTTTTTTACCTCGGTATGCCTTGTGACTGGCTTGACAAGCTTTTCCTCAGCTCCGTAAACTCCTTTCAGTGGGAAATTGTGGGGCAAAGTGGGAATAAGGGGTGAGGCTGGCATGTTCCGGGGAGCAACGTTAGTCAATCTCGACAGCAAAGGGCGCTTATCAGTGCCTACCCGTTATCGGGAACAGCTGCTTGAGAACGCTGCCGGTCAAATGGTTTGCACCATTGACATTTATCACCCGTGCCTGCTGCTTTACCCCCTGCCTGAATGGGAAATTATC  >  NC_000913/89327‑89801

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑