Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,657,830 C→A 13.8% R50L (CGA→CTA)  gadW ← transcriptional activator of gadA and gadBC; repressor of gadX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,657,8300CA13.8% 80.9 / 4.6 29R50L (CGA→CTA) gadWtranscriptional activator of gadA and gadBC; repressor of gadX
Reads supporting (aligned to +/- strand):  ref base C (13/12);  new base A (2/2);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

TGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAGGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATC  >  CP009273/3657585‑3658092
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                      
tGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGTGACTGGTTAATTAACAAACGTTCTTTGCTGTGTAAGCATAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTATCAGCTCATCGATATCCAGTCGTCTCGTATCCGATGAACAAAACGCaaaa                                                                                                                                                                                                                                           >  1:268830/1‑275 (MQ=255)
                 aGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACACAAAATCATCTGCCAAAttttt                                                                                                                                                                                                                          >  1:200460/1‑275 (MQ=255)
                       tttttAAGANAAAAACTCGATACACATGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCcag                                                                                                                                                                                                                    <  2:200453/272‑1 (MQ=255)
                           tATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATATTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAACCGCAAAATCGTCTGCCAAATTTTTCTCcagcagc                                                                                                                                                                                                                >  1:251959/1‑275 (MQ=255)
                                 aaaaCATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAgtcgt                                                                                                                                                                                                                                                                      <  1:196599/215‑1 (MQ=255)
                                  aaaaATAGAGAAAGAGGAAAGGTAGAGCATATGGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAAGCGCCGTCACCAGCGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATAc                                                                                                                                                                                                         <  2:367593/275‑1 (MQ=255)
                                                                               tGAATTTACACAGGATTCATTGAAGCTATCAAAAATCGGCGTCACCAGCGGTATGGGGGATGACTGGTTAATGAACAAACGGTCTTTGCTGTGTAACACTAAAGTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAACCTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGa                                                                                                                                                            <  1:57075/275‑1 (MQ=255)
                                                                                       ccctGGATTTTTTGAAGCTATCAAAAATAGCCGTCACCAGCGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAacg                                                                                                                                                     <  1:36974/274‑1 (MQ=255)
                                                                                                                    gCCGTCAACAGCGGCATTGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCATACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATAGTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAGGATCACCGAGCAGACATGAGt                                                                                                                       <  1:312311/275‑1 (MQ=255)
                                                                                                                     ccGTCACCAGCGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGATGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCt                                                                                                                                                                                                                                     >  1:180370/1‑164 (MQ=255)
                                                                                                                         cACCAGCGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCa                                                                                                                                                                                                                                 <  2:180363/164‑1 (MQ=255)
                                                                                                                                            gACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAGGATCAc                                                                                                                                      >  1:267580/1‑236 (MQ=255)
                                                                                                                                                 gTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGATGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTGGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCGAGCAGGAGAATACTCTCGTTATGCAGAGATATTTTTTGCTGTTCAAGATAAATAGCGAATGAACGAGGAATGAGGATCACCAAGCAAACATGAGTCATGATAATACCTTATATTTCATACTgag                                                                                          >  1:269667/1‑273 (MQ=255)
                                                                                                                                                  ttaattaaTAAGTGGTTTTTGATTNGTAACCCTCAATTGCGTGGTAACTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAGGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGa                                                                                         <  2:307027/275‑1 (MQ=255)
                                                                                                                                                                         gtgtAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGc                                                                                                                                                                                              >  2:421577/1‑151 (MQ=255)
                                                                                                                                                                           gtAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTAGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAAttttt                                                                                                                                                                                                                          >  1:234611/1‑121 (MQ=255)
                                                                                                                                                                             aaCCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCg                                                                                                                                                                                          <  1:421592/151‑1 (MQ=255)
                                                                                                                                                                               cccTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAGGATCACCGAGCAGACATGAGTCATGAtt                                                                                                                <  2:174441/223‑1 (MQ=255)
                                                                                                                                                                               cccTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTAGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTcc                                                                                                                                                                                                                      <  2:234602/121‑1 (MQ=255)
                                                                                                                                                                                    aaTTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAAttttt                                                                                                                                                                                                                          >  1:151596/1‑112 (MQ=255)
                                                                                                                                                                                        gCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAACACGCAAAATCGTCTGCCAAATTTTTCTcc                                                                                                                                                                                                                      <  2:151593/112‑1 (MQ=255)
                                                                                                                                                                                         cGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAGGATCACCGAGCa                                                                                                                                >  2:432281/1‑197 (MQ=255)
                                                                                                                                                                                             gTAGCTTACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCAGCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAGGATCACCGAGCAGACa                                                                                                                            <  1:432296/197‑1 (MQ=255)
                                                                                                                                                                                                        atatTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGaa                                                                                                                                                           >  1:147582/1‑155 (MQ=255)
                                                                                                                                                                                                            tttGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGaa                                                                                                                                                       <  2:147579/155‑1 (MQ=255)
                                                                                                                                                                                                                   gtaaTGGCAAACTGTCAGCTCATCGATATCCAGTAGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAacg                                                                                                                                                     >  2:70956/1‑150 (MQ=255)
                                                                                                                                                                                                                       tGGCAAACTGTCAGCTCATCGATATCCAGTAGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAAcgacga                                                                                                                                                 <  1:70956/150‑1 (MQ=255)
                                                                                                                                                                                                                            aaCTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAGGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATAAAATGATGCCAGACTGTTTCTTATTGCAAAAACTGCCAGATTTggaaag               >  2:176586/1‑271 (MQ=255)
                                                                                                                                                                                                                                          cGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTATGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAGGATCACCGAGCAGACATGAGTAATGATTATCCATTATATTTCATACTTCGATTATTTCAATTTCAGTATAATAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCAAAAATTATAAGGACCAACATTTATc  >  1:210618/1‑274 (MQ=255)
                                                                                                                                                                                                                                               tCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAGGATCAc                                                                                                                                      >  1:192871/1‑137 (MQ=255)
                                                                                                                                                                                                                                                   gtcgtcGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAGGATCACCGAg                                                                                                                                  <  2:192864/137‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                      
TGAAAAGTAAGGGGATCAGTTCTTTCTTATGAGAAAACATCGATAAACAGGAAAGGTAGAGCATATTGCTCAGTATCGGTGAATTTACCCCGGATTCATTGAAGCTATCAAAAATCGCCGTCACCAGCGGCATGGGGGGTGACTGGTTAATTAACAAACGGTCTTTGCTGTGTAACCCTAAATTGCGTGGTAGCTGACGAATATTTTGTAAGTAATGGCAAACTGTCAGCTCATCGATATCCAGTCGTCGCGTATCCGGTGAACAAAACGCAAAATCGTCTGCCAAATTTTTCTCCAGCAGCAGAATACTCTCGTTATGCAGCGATATTTTTTGCTGTTCATGATAAATATCGAATGAACGACGAATGAGGATCACCGAGCAGACATGAGTCATGATTATCCCTTATATTTCATACTGCGATTATTTCAATTTCAGTATACTAATGAAATGATGCCAGACTGTTTCTTAGCGCAAAAACTGCCAGATTTGGTAGGACCAAAATTTATC  >  CP009273/3657585‑3658092

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