Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 468,093 Δ1 bp 33.3% coding (40/339 nt) glnK → nitrogen assimilation regulatory protein for GlnL, GlnE, and AmtB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273468,0920G.33.3% 39.6 / 23.7 21coding (39/339 nt)glnKnitrogen assimilation regulatory protein for GlnL, GlnE, and AmtB
Reads supporting (aligned to +/- strand):  ref base G (7/6);  new base . (4/3);  total (11/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TGCTCTCTTTCTCGTTTTTCATTTCTGGCACACCGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAA  >  CP009273/467961‑468232
                                                                                                                                   |                                                                                                                                            
tGCTCTCTTTCTCGTTTTTCATTTCTGGCACACCGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGc                                                                                                                             >  1:241725/1‑149 (MQ=255)
      ctTTCTTGTTGTTCCTTTCTGGCACACCGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCt                                                                                                                      <  2:315555/149‑1 (MQ=255)
               ttttCATTTCTGGCACACCGCTTGCCATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTCACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTGGt                                                                                                             <  1:1003/149‑1 (MQ=255)
                     cttCTGGCACACCGTTTGCAATACCGTCTTCGTGTAGCGGAAACTTTCCCAAATTCGGCCCGGGGGTGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCGTGAAGACGTTCGTGGAGCGTTATCTTCCATTGGTATTCa                                                                                                        <  1:85337/148‑1 (MQ=255)
                      ttccgggccacCGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAg                                                                                                       <  1:270001/141‑1 (MQ=255)
                        cTGGCACACCGCTTGCAATACCTTCTTCGTGGAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGc                                                                                                    <  1:228757/149‑1 (MQ=255)
                                    ttGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTGGTGAAGCGTTCTCTTCCATTGGTATTCAGGGCCTGACCGTccc                                                                                         >  1:200823/1‑147 (MQ=255)
                                         aTACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGGGGGGATCTATGAAGCTGGTGACCGGGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAg                                                                                    <  1:99571/149‑1 (MQ=255)
                                               ttgtcGTGTAGCAGAACCATTTCCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAg                                                                              <  1:135539/146‑1 (MQ=255)
                                               tcctcGTGTAGCAGACTCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAg                                                                              <  1:239793/146‑1 (MQ=255)
                                                  tcGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTCATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGGGAAAGGtt                                                                           >  1:304218/1‑149 (MQ=255)
                                                   cGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGttt                                                                          <  2:270110/149‑1 (MQ=255)
                                                         gCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCg                                                                    >  1:96724/1‑149 (MQ=255)
                                                         gCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCg                                                                    >  2:112736/1‑149 (MQ=255)
                                                                       gAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATg                                                      <  2:241725/149‑1 (MQ=255)
                                                                                      gggATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCggg                                       >  2:108680/1‑149 (MQ=255)
                                                                                         aTCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGAGAAGGGTACCGGACGTCAGAAAGGGCATGCCGAGCTGTACCggggggg                                   >  2:333463/1‑148 (MQ=255)
                                                                                                                 aaTCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGCCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGCGCGGAATACAGCGACAATTTCCTGcc           >  1:29780/1‑149 (MQ=255)
                                                                                                                     aaaCCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCaaaa       >  1:304934/1‑149 (MQ=255)
                                                                                                                      aaCCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAATGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGCCAATTTCCTGCCaaaa       >  2:201576/1‑149 (MQ=255)
                                                                                                                          aTTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACGGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTaaa  >  1:96763/1‑149 (MQ=255)
                                                                                                                                   |                                                                                                                                            
TGCTCTCTTTCTCGTTTTTCATTTCTGGCACACCGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAA  >  CP009273/467961‑468232

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