Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 496,458 C→A 17.0% S220I (AGC→ATC)  gltC ← glutamate synthase operon transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP000730496,4580CA17.0% 74.3 / 6.7 28S220I (AGC→ATC) gltCglutamate synthase operon transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (13/10);  new base A (2/3);  total (15/13)
Fisher's exact test for biased strand distribution p-value = 6.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.14e-01

AGCACATATCGTTCTTTTTAATTCAGTATGATCTAATTTTATATCTATCCATGATTTAGATTCTGGTAAATGTATATTTTGTGATGAAATGATGTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGCTACTGATTGTGTATACATCTGATTTAGTAATCTTGCGCAATTGTTTTTTTACAGTTTCGGCAAATGGTGCCAAGCAATAAATATGACTATGCTCAAACTGAATTAATGGTGGGTGTGTCGCCATCGTAATTGGATCGTC  >  CP000730/496322‑496597
                                                                                                                                        |                                                                                                                                           
aGCACATATCGTTCTTTTTAATTCAGTATGATCTAATTTTATATCTATCCATGATTTAGATTCTGGTAAATGTATATTTTGTGATGAAATGATGTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGCTACt                                                                                                                                         >  2:108927/1‑141 (MQ=255)
aGCACATATCGTTCTTTTTAATTCAGTATGATCTAATTTTATATCTATCCATGATTTAGATTCTGGTAAATGTATATTTTGTGATGAAATGATGTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGCTACt                                                                                                                                         >  2:130915/1‑141 (MQ=255)
         cGTTCTTTTTAATTCAGTATGATCTAATTTTATATCTATCCATGATTTAGATTCTGGTAAATGTATATTTTGTGATGAAATGATGTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGCTACTGATAGTGta                                                                                                                                >  2:240054/1‑141 (MQ=255)
                tttAATTCAGTATGATCTAATTTTATATCTATCCATGATTTAGATTCTGGTAAATGTATATTTTGTGATGAAATGATGTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGCTACTGATTGTGTATACATCt                                                                                                                         >  1:105668/1‑141 (MQ=255)
                      tCAGTATGATCTAATTTTATATCTATCCATGATTTAGATTCTGGTAAATGTATATTTTGTGATGAAATGATGTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGCTACGTATTGTGTATACATCTGATTTa                                                                                                                   <  2:377540/141‑1 (MQ=255)
                       cAGTATGATCTAATTTTATATCTATCCATGATTTAGATTCTGGTAAATGTATATTTTGTGATGAAATGATGTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGCTACTGATTGTGTATACATCTGATTTAg                                                                                                                  >  1:101530/1‑141 (MQ=255)
                       cAGTATGATCTAATTTTATATCTATCCATGATTTAGATTCTGGTAAATGTATATTTTGTGATGAAATGATGTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGATACTGATTGTGTATACATCTGATTTAg                                                                                                                  <  2:157511/141‑1 (MQ=255)
                           aTGATCTAATTTTATATCTATCCATGATTTAGATTCTGGTAAATGTATATTTTGTGATGAAATGATGTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTCGTTGGCTACTGATTGTGTATACATCTGATTTAGTAAt                                                                                                              <  1:272498/141‑1 (MQ=255)
                                  aaTTTTATATCTATCCATGATTTAGATTCTGGTAAATGTATATTTTGTGATGAAATGATGTAAACTTCTTTTTGCCAAAGGAGATACTGCGCAAGTGGTTGGCTACTGCTTGTGTATACATCTGATTTAGTAATCTTgcgc                                                                                                       >  1:93542/1‑141 (MQ=255)
                                    ttttATATCTATCCATGATTTAGATTCTGGTAAATGTATATTTTGTGATGAAATGATGTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGCTACTGATTGTGTATACATCTGATTTAGTAATCTTGCGCaa                                                                                                     >  2:197490/1‑141 (MQ=255)
                                       tataTCTATCCATGATTTAGATTCTGGTAAATGTATATTTTGTGATGAAATGATGTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGCTACTGATTGTGTATACATCTGATTTAGTAATCTTGCGCAATTg                                                                                                  <  1:59843/141‑1 (MQ=255)
                                        atatCTATCCATGATTTAGATTCTGGTAAATGTATATTTTGTGATGAAATGATGTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGATACTGATTGTGTATACATCTGATTTAGTAATCTTGCGCAATTGt                                                                                                 <  2:162806/141‑1 (MQ=255)
                                           tctatcCATGATTTAGATTCTGGTAAATGTATATTTTGTGATGAAATGATGTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGCTACTGATTGTGTATACATCTGATTTAGTAATCTTGCGCAATTGtttt                                                                                              >  1:32310/1‑141 (MQ=255)
                                            ctatcCATGATTTAGATTCTGGTAAATGTATATTTTGTGATGAAATGATGTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGCTACTGATTGTGTATACATCTGATTTAGTAATCTTGCGCAATTGttttt                                                                                             >  2:212817/1‑141 (MQ=255)
                                             tatcCATGATTTAGATTCTGGTAAATGTATATTTTGTGATGAAATGATGTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGCTACTGATTGTGTATACATCTGATTTAGTAATCTTGCGCAATTGtttttt                                                                                            >  2:157434/1‑141 (MQ=255)
                                               tcCATGATTTAGATTCTGGTAAATGTATATTTTGTGATGAAATGATGTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGCTACTGATTGTGTATACATCTGATTTAGTAATCTTGCGCAATTGTTTTTTTa                                                                                          <  1:346592/141‑1 (MQ=255)
                                                      tttAGATTCTGGTAAATGTATATTTTGTGATGAAATGATGTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGCTACTGATTGTGTATACATCTGATTTAGTAATCTTGCGCAATTGTTTTTTTACAGTTTc                                                                                   >  1:19050/1‑141 (MQ=255)
                                                              cTGGTAAATGTATATTTTGTGATGAAATGATGTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGCTACTGATTGTGTATACATCTGATCTAGTAATCTTGCGCAATTGTTTTTTTACAGTTTCGGCAAATg                                                                           >  2:360784/1‑141 (MQ=255)
                                                                                         tgatgTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGCTACTGATTGTGTATACATCTGATTTAGTAATCTTGCGCAATTGTTTTTTTACAGTTTCGGCAAATGGTGCCAAGCAATAAATATGACTATGCt                                                <  1:77310/141‑1 (MQ=255)
                                                                                          gatgTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGATACTGATTGTGTATACATCTGATTTAGTAATCTTGCGCAATTGTTTTTTTACAGTTTCGGCAAATGGTGCCAAGCAATAAATATGACTATGCTc                                               <  2:24354/141‑1 (MQ=255)
                                                                                           atgTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGCTACTGATTGTGTATACATCTGATTTAGTAATCTTGCGCAATTGTTTTTTTACAGTTTCGGCAAATGGTGCCAAGCAATAAATATGACTATGCTCa                                              <  1:303651/141‑1 (MQ=255)
                                                                                              tAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGCTACTGATTGTGTATACATCTGATTTAGTAATCTTGCGCAATTGTTTTTTTACAGTTTCGGCAAATGGTGCCAAGCAATAAATATGACTATGCTCAAAc                                           >  1:89351/1‑141 (MQ=255)
                                                                                                       ttttGACGAAGGAGATACTGCGCAAGTGGTTGGCTACTGATTGTGTATACATCTGATTTAGTAATCTTGCGCAATTGTTTTTTTACAGTTTCGGCAAATGGTGCCAAGCAATAAATATGACTATGCTCAAACTGAATTAAt                                  <  1:333537/141‑1 (MQ=255)
                                                                                                                  gagaTACTGCGCAAGTGGTTGGCTACTGATTGTGTATACATCTGATTTAGTAATCTTGCGCAATTGTTTTTTTACAGTTTCGGCAAATGGTGCCAAGCAATAAATATGACTATGCTCAAACTGAATTAATGGTGGgtgtgt                       <  1:108927/141‑1 (MQ=255)
                                                                                                                        cTGCGCAAGTGGTTGGCTACTGATTGTGTATACATCTGATTTAGTAATCTTGCGCAATTGTTTTTTTACAGTTTCGGCAAATGGTGCCAAGCAATAAATATGACTATGCTCAAACTGAATTAATGGTGGGTGTGTCGCCAt                 <  2:37589/141‑1 (MQ=255)
                                                                                                                                    ttGGATACTGATTGTGTATACATCTGATTTAGTAATCTTGCGCAATTGTTTTTTTACAGTTTCGGCAAATGGTGCCAAGCAATAAATATGACTATGCTCAAACTGAATTCATGGTGGGTGTGTCGCCATCGTAATTGGATc     >  2:128021/1‑141 (MQ=255)
                                                                                                                                    ttGGATACTGATTGTGTATACATCTGATTTAGTAATCTTGCGCAATTGTTTTTTTACAGTTTCGGCAAATGGTGCCAAGCAATAAATATGACTATGCTCAAACTGAATTAATGGTGGGTGTGTCGCCATCGTAATTGGATc     >  2:151792/1‑141 (MQ=255)
                                                                                                                                       gCTACTGATTGTGTATACATCTGATTTAGTAATCTTGCGCAATTGTTTTTTTACAGTTTCGGCAAATGGTGCCAAGCAATAAATATGACTATGCTCAAACTGAATTAATGGTGGGTGTGTCGCCATCGTAATTGGATCGTc  <  1:360784/141‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                           
AGCACATATCGTTCTTTTTAATTCAGTATGATCTAATTTTATATCTATCCATGATTTAGATTCTGGTAAATGTATATTTTGTGATGAAATGATGTAACCTTCTTTTTGACGAAGGAGATACTGCGCAAGTGGTTGGCTACTGATTGTGTATACATCTGATTTAGTAATCTTGCGCAATTGTTTTTTTACAGTTTCGGCAAATGGTGCCAAGCAATAAATATGACTATGCTCAAACTGAATTAATGGTGGGTGTGTCGCCATCGTAATTGGATCGTC  >  CP000730/496322‑496597

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