Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,664,594 G→A L9F (CTC→TTC)  gadW ← DNA‑binding transcriptional dual regulator GadW

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,664,5940GA90.5% 52.4 / ‑1.2 23L9F (CTC→TTC) gadWDNA‑binding transcriptional dual regulator GadW
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (9/10);  total (11/12)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-02

GCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAGGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGATAATAAGTAAATGCAGCACGAATATATTTTCGCACAGCGTATAGCTTATGTTTATAAAAAAATGGCT  >  NC_000913/3664354‑3664830
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
gCATGGGGGGTGACTGGTTAAGTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAg                                                                                                                                                                                                                                     <  1:450685/251‑1 (MQ=255)
gCATGGGGGGGGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAAGCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTAGGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAg                                                                                                                                                                                                                                     <  2:136031/251‑1 (MQ=255)
               gtttaattatCAAACGGTCTGTGCTTGGGAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAATTAATGGCAAACTGTTAGCTCATCGATATCCAGTCTTCGCGTAGCCGGTGAACAAAACGCAAAATCGTCTGCTAAATTTTTCTCCAGCAGCAGAATTCTCTCGTTATGCAGCGAGATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCGGACATGAGTCAt                                                                                                                                                                                                                       <  2:286977/248‑1 (MQ=255)
                 ttaattaaCAAAAGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGACTGAAGATAACCGAGCAGACATGAGTCATgtt                                                                                                                                                                                                                    >  1:291114/1‑249 (MQ=255)
                     ttaaCAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGTTTATc                                                                                                                                                                                                                >  2:185277/1‑251 (MQ=255)
                          aaaCGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAGGATCACCGAGCAGACATGAGTCATGATTATCCCTta                                                                                                                                                                                                           <  1:96629/251‑1 (MQ=255)
                                  tttGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCCTAGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGAAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGGGTCATGATTATCCCTTATATTTCAt                                                                                                                                                                                                   >  1:447134/1‑251 (MQ=255)
                                                                                    aaTGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGtt                                                                                                                                                 >  1:344010/1‑251 (MQ=255)
                                                                                                         cGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAg                                                                                                                            <  1:154312/251‑1 (MQ=255)
                                                                                                                 aggcgtcgCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCACGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTaatgaaat                                                                                                                                                               <  2:99886/205‑1 (MQ=255)
                                                                                                                 aGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTaatgaaat                                                                                                                                                               >  1:99886/1‑208 (MQ=255)
                                                                                                                                         aaCGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAACTTATCCTCAAagta                                                                                           <  2:241921/251‑1 (MQ=255)
                                                                                                                                              aaaaTCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGAtt                                                                                                                                                                                                                   >  1:73916/1‑127 (MQ=255)
                                                                                                                                              aaaaTCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGAtt                                                                                                                                                                                                                   <  2:73916/127‑1 (MQ=255)
                                                                                                                                                             aTTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAAtt                                                                       >  1:371960/1‑251 (MQ=255)
                                                                                                                                                                      cagcaccagAATACTATCGTTATGCAGCGATATTTTGTGATGTTCATGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCATACATGAGTCATGATTCTCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAACCTGCCAGATTTTGTAGGACCAAAATTTCTCCACAAAGTAGTTTTGCATAACAGTTAAGTGAtaataac                                                              <  2:215478/251‑2 (MQ=255)
                                                                                                                                                                                       cGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAGGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGATAATAAGTAAATGCAGCACGAata                                             <  2:164678/251‑1 (MQ=255)
                                                                                                                                                                                                      ttttttGCTGTTCATGATACATATCGAATGAACGACGAATGACGTTCACCGTGCAGACATGAGCAATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTACTGAAATGATGCCAGACTGTTTCTCAGCGTAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGCTAATAAGTAAATGCAGCACGAATATATTTTCGCAAAGCg                              >  2:272940/1‑251 (MQ=255)
                                                                                                                                                                                                                   aTGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGATAATAAGTAAATGCAGCACGAATATATTTTCGCACAGCGTATAGCTTATGtt                 <  1:285586/251‑1 (MQ=255)
                                                                                                                                                                                                                   aTGATAAATATCGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGATAATAAGTAAATGCAGCACGAATATATTTTCGCACAGCGTATAGCTTATGtt                 >  1:259417/1‑251 (MQ=255)
                                                                                                                                                                                                                              cGAATGAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGATAATAAGTAAATGCAGCACGAATATATTTTCGCACAGCGTATAGCTTATGTTTATAAAAAAAt      >  1:194965/1‑251 (MQ=255)
                                                                                                                                                                                                                                 aTGAACGACGAATGACGAACACCGAGCAGCCATGAGTCGTGATTATACCTTATATTTCACACTGCGATGATTTCACTTTCAAAATACTACTGAAAAGCTAACAGACTGATTCTTAGCGCAAAGACTGCCAGATTTGACGGAGCCAAAATTTAGCCACACAGAACTTTTGCATAAAAGTTAATTGATAATAAGTATGTACAGCACGAATATATTTTCGCACAGCGTATAGCTTAACTTTATAAAAAACTGGc   >  2:401872/1‑251 (MQ=255)
                                                                                                                                                                                                                                  agAACGACGAATGAAGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGATAATAAGTAAATGCAGCACGAATATATTTTCGCACAGCGTATAGCTTATGTTTATAAAAAAATGGCt  <  1:151093/250‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                             
GCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAGGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATCCACAAAGTAGTTTTGCATAACAGTTAATTGATAATAAGTAAATGCAGCACGAATATATTTTCGCACAGCGTATAGCTTATGTTTATAAAAAAATGGCT  >  NC_000913/3664354‑3664830

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: