Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,259,611 G→A A131A (GCC→GCT zur ← DNA‑binding transcriptional repressor Zur

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,259,6110GA95.8% 57.5 / NA 24A131A (GCC→GCTzurDNA‑binding transcriptional repressor Zur
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (8/15);  minor base C (1/0);  total (9/15)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

ATCAGAAAGATCAGGCAGAAAAAGAGGTCGTGGATTGGGAAACCCGCAATGAATATCGCTGGTAATTAATCCCTCCTGCCCGACGTGTACAAGGATGTACGCCCTCTTAACGCGGTTTCTTTTTCACCTGCACAGAGTGATCATGCTGGCACTGTTCAGGATGACGACACGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGTCTTCCACGCCTTCTGCACACTCTTCTTTCACTGCGCCGCAGCGATCGCAAATAAACATGGCTGACGTATGGGTGGGCTGATCGAACAGATGACAGAGCACATAACTGTTGGTGGATTCCACCTTATGCACAAAACCTTGCTCAAGCAGAAAATCCAGCGCGCGATAAACCGTTGGCGGCTTGGCTTGCGGTTCAGCTTCGCGCAGT  >  NC_000913/4259382‑4259854
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
aTCAGAAAGATCAGGCAGAAAAAGAGGTCGTGGATTGGGAAACCCGCAATGAATATCGCTGGTAATTAATCCCTCCTGCCCGACGTGTACAAGGATGTACGCCCTCTTAACGCGGTTTCTTTTTCACCTGCACAGAGTGATCATGCTGGCACTGTTCAGGATGACGACACGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGAGCAAACCCCATTTTTgccgn                                                                                                                                                                                                                                 <  1:274788/250‑2 (MQ=255)
                                       aaaCCCGCAATGAATAACGCTGGTAATTAAGCCCTCCTGCCCGACGTGTACAAGGATGTACGCCCTCTTAACGCGGTTTCTTTTTCACCTGCACAGAGTGATCATGCTGGCACTGTTCAGGATGACGACACGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGAGCAAACCCCATTTTTGCCGCCATCGTATGCATAATGTCTTCCACGCCTTCTGCACActct                                                                                                                                                                                         <  1:196512/251‑1 (MQ=255)
                                          cccGCAATGAATATCGCTGGTAATTAATCCCTCCTGCCCGACGTGTACAAGGATGTACGCCCTCTTAACGCGGTTTCTTTTTCACCTGCACAGAGTGATCATGCTGGCACTGTTCAGGATGACGACACGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGAGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGTCTTcc                                                                                                                                                                                                           >  1:385190/1‑230 (MQ=255)
                                                            ggTAATTAATCCCTCCTGCCCGACGTGTACAAGGATGTACGCCCTCGTAACGCGGTTTCTTTTTCACCTTCACAGAGTGATCATGCTGGCACTGTTCAGGATTACGACACGCTTCCACTTCTACACATTCCGCACAGATCCCATGTGCTTCAATCACATTATGCCGCAGAGCAAACCCCATTTTTGCCGCCAGCGTATTCATAATGTCGTCCACGCCTTCTGCACACTCTTCTTTCACTGCGCCGCAGCg                                                                                                                                                                     <  1:253817/250‑1 (MQ=255)
                                                                                      gtACAAGGATGTACGCCCTCTTAACGCGGTTTCTTTTTCACCTGCACAGAGTGATCATGCTGGCACTGTTCAGAATGACGACACGCTTCCACTTCTACACATGCAGCACAGCGCCAATGTGCTTCAATCACATTATGCCGCAGAGCAATCCCCATTTTTGCCGCCAGCGTATGCATAATGTCTTCCACGCCTTCTGCACACTCTTCTTTCACTGCGCCGCAGCGATCGCAAATAAACATGGCTGACGtctg                                                                                                                                          >  1:431772/1‑248 (MQ=255)
                                                                                            gtaTGTACGCCCTCTTAACGCGGTTTCTTTTTCACCTGCACAGAGTGATCATGATGGAACTGTTAAGGATGACGACAAGCTTCCCATTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGAGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGTCTTCCACGCCTTCTGCACACTCTTCTTTCACTGCGCCGCAGCGATCGCAAATAAACATGGCTGACGTAtgggtggg                                                                                                                                    >  1:416580/3‑251 (MQ=255)
                                                                                                  aCGCCCTCTTAACGCGGTTTCTTTGTCGCCTGCACAGAGTGAGCGTGCGGGCGCTGTTCAGGATGACGACTCGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATTCCGCAGAGCAAACGCCATTGTTGCCGCCAGCGTATGCATAATGTCTTCCACGCCTTCGGCGGACTCTTCTTTCACTGCGCCGCAGCGATCGCAAATAAACATGGCTGACGTATGGGTGGGCTGAt                                                                                                                               <  1:254002/250‑1 (MQ=255)
                                                                                                    gCCCTCTTAACGCGGTTTCTTTTTCACCTGCACAGAGTGATCATGCTGGCACTGTTCAGGATGACGACACGCTTCCACTTCTACACATGCCGCACAGAACCCATGTGCTTCAATCACATTATGACGCAGAGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGTCTTCCACGCCTTCTGCACACTCTTCTTTCACTGCGCCGCAGCGATCGCAAATAAACATGGCTGACGTATGGGTGGGCTGATCg                                                                                                                             >  1:187840/1‑250 (MQ=255)
                                                                                                      ccTCTTAACGCGGTTTCTTTTTCACCTGCACAGAGTGATCATGCTGGCACTGTTCAGGATGACGACACGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGAGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGTCTTCCAc                                                                                                                                                                                                         >  1:252632/1‑172 (MQ=255)
                                                                                                               gcgGTTTCTTTTTCACCGGCCCAGAGTGATCATGCTGGCACTGTTCAGGATGACGACACGCTGCCGCTTCTACACATGCCGCACAGAGCCCATGTGCTGCAAGCACATTATGCCGCAGAGCAAACCCCATTTTGGCCTCCAGCTTAGGCATAATGTCTGCCGCGGCTTCTTCACACTCGTCTATCACTgc                                                                                                                                                                              <  1:367889/190‑1 (MQ=255)
                                                                                                                        ttttCTCCGGCACAGAGTGATCATGCTGGCACTGTTCAGGATGACGACACGCTTCCACTTCTACACATGCCGCACATAGCTCATGTGCTCCAAACACACTAGGCCGCAGAGCAAACCCCATTTTTGCCGCCAGCGTATGCGTAATGTCTACCACGCCTTCTGCACACTCGTCTTGCACTGCGCCGCAGCg                                                                                                                                                                     <  1:299423/190‑1 (MQ=255)
                                                                                                                                                     cACTGTTCAGGATGACGACACGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGAGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGGCTTCCACGCCTTCTGCACACTCTTCGGTCACTGCGCCGCAGCGATCGCAAATAAACATGGCTGACGTAtgggtgg                                                                                                                                     <  1:385064/193‑1 (MQ=255)
                                                                                                                                                       cTGTTCAGGATGACGACACGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGAGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGTCTTCCACGCCTTCTGCACACTCTTCTTTCACTGCGCCGCAGCGATCGCAAATAAACATGGCTGACGTATGGGTGGGCTGATCGAACAGATGACAGAGCACATAACTGTTGGTGGATTCCACCTTATGCACAAAAc                                                                         >  1:29078/1‑251 (MQ=255)
                                                                                                                                                       cTGTTCAGGATGACGACACGCTTCCACTTCTACACATGCAGCACAGAGCACATGTGCTTCACTAACATTATGCAGTAGAGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGTCTTCCAc                                                                                                                                                                                                         >  1:323915/1‑123 (MQ=255)
                                                                                                                                                       cTGTTCAGGATGACGACACGCTTCCACTNCTACACAGGCNGCANNGAGCNCNTNTGCTTCAATCACATTATGCCGCAGAGCAANCCCCATTTTTGCCGCCAGCGTGTNNTTNNNGTNNTCCAc                                                                                                                                                                                                         <  2:323915/123‑1 (MQ=255)
                                                                                                                                                             aGGATGACGACACGCTTCCACGTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGAGCAAACCCCAttttt                                                                                                                                                                                                                                      <  1:409603/88‑1 (MQ=255)
                                                                                                                                                              ngATGACGACACGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCATAGCAAACCCCATTTTTGCCGCCATCGTATg                                                                                                                                                                                                                        >  1:262806/2‑101 (MQ=255)
                                                                                                                                                              ggATGACGACACGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAACCACAGTATGCCGCAGAGCAAACCCCATCTTTgccgactgcgtcgg                                                                                                                                                                                                                        <  2:262806/101‑11 (MQ=255)
                                                                                                                                                                                     tACACATGCCGCACAGAGCCCATGTTCTTCAATCACATTATGCCGCAGAGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGTCTGCCACGCCTTCGGCACACTCTTCTTTCACTGCGCCGCAGCGATCGCAAATAAACATGGCtn                                                                                                                                                 <  1:48299/149‑2 (MQ=255)
                                                                                                                                                                                            gCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGAGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGTCTGCCACGCCTTCTGCACACTCTTCTTGCACTGCGCCGCAGCGATCGCAAATAAACATGGCTGACGTATGGGTGGGCGGATCGAACAGATGACAGAGCACATAACTGTTGGTGGATTCCACCTTATGCACAAAACCTTGCTCAAGCAGAAAATCCAGCGCGCGAt                                           <  1:393978/244‑1 (MQ=255)
                                                                                                                                                                                                                  cAATCACATTATGCCGCAGAGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGTCTTCCACGCCTTCTGCACACTCTTCTTTCACTGCGCCGCAGCGATCGCAAATAAACATGGCTGACTTATGGGTGGGCTGATCGAACAGATGACAGAGCACATAACTGTTGGTGGATTCCACCTTATGcac                                                                              <  1:288982/187‑1 (MQ=255)
                                                                                                                                                                                                                        caTTATGCCGCAGAGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGTCTTCCACGCCTTCTGCACACTCTTCTTTCACTGCGCCGCAGCGATCGCAAATAAACATGGCTGACGTATGGGTGGGCTGATCGAACAGATGACAGAGCACATAACTGTTGGTGGATTCCACCTTATGCACAAAACCTTg                                                                     <  1:326545/190‑1 (MQ=255)
                                                                                                                                                                                                                        caTTANNCCNNNGCNNAAACCCCATTTTTGCCGCCAGCGNATGCATAATGTCTTCCACGCCTTCTGCACNCNCNTCTTNNACTNCGCCGCAGCGNTCGCAAATAAACATGGCTGACGTATGGGTGGGCTGATCGAACAGATGANANAGCACATAACTGTTGGTGGCTTCAACCTTATNCACAAAACCTTg                                                                     >  2:326545/1‑190 (MQ=255)
                                                                                                                                                                                                                              gCCGCAGAGCAAACCCCATTTTTGCCGACAGCGTATGCATAATGTCTGCCACGCCGTCTGCACACTCTTCTTTCACTGCGCCGCAGCGATCGCAAATAAACAATGCTGACGTATGGGTGGGCTGATCTAACAGATGACAGAGCACAGAACGGTTGGTGGATTCCACTTTAGGCACAAAACCTTGCTGAAGCAGAAAAGCCAGCGCGCGAGAAACCGTTGGCGGCGTGGCGTGCGGTTCAGCTTCGCGCAGt  <  1:445089/251‑1 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   
ATCAGAAAGATCAGGCAGAAAAAGAGGTCGTGGATTGGGAAACCCGCAATGAATATCGCTGGTAATTAATCCCTCCTGCCCGACGTGTACAAGGATGTACGCCCTCTTAACGCGGTTTCTTTTTCACCTGCACAGAGTGATCATGCTGGCACTGTTCAGGATGACGACACGCTTCCACTTCTACACATGCCGCACAGAGCCCATGTGCTTCAATCACATTATGCCGCAGGGCAAACCCCATTTTTGCCGCCAGCGTATGCATAATGTCTTCCACGCCTTCTGCACACTCTTCTTTCACTGCGCCGCAGCGATCGCAAATAAACATGGCTGACGTATGGGTGGGCTGATCGAACAGATGACAGAGCACATAACTGTTGGTGGATTCCACCTTATGCACAAAACCTTGCTCAAGCAGAAAATCCAGCGCGCGATAAACCGTTGGCGGCTTGGCTTGCGGTTCAGCTTCGCGCAGT  >  NC_000913/4259382‑4259854

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