| 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 | ||
| * | CP000730 | 496,458 | 0 | C | A | 17.0% | 74.3 / 6.7 | 28 | S220I (AGC→ATC) | gltC | glutamate 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: ‑ |