Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 468,126 Δ1 bp 39.1% coding (73/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,1260A.39.1% 30.9 / 28.3 23coding (73/339 nt)glnKnitrogen assimilation regulatory protein for GlnL, GlnE, and AmtB
Reads supporting (aligned to +/- strand):  ref base A (8/6);  new base . (3/6);  total (11/12)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

TTTCTGGCACACCGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGAT  >  CP009273/467982‑468274
                                                                                                                                                |                                                                                                                                                    
tttCTGGCACACCCCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTTATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTTTCAg                                                                                                                                                 <  1:238615/149‑1 (MQ=255)
 ttCTGGCACACCGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATTATCAAATCATTTAAGCTGGAAGACGTTCGTGAAGCGTTATCTCCCATTGGTTTCAgg                                                                                                                                                <  1:28957/149‑1 (MQ=255)
      gCACACCGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTATGACAGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGCAGCGTTATCTTCCATTGGTATTCAGGGCCt                                                                                                                                            >  1:84959/1‑149 (MQ=255)
       cacaCCGCTTGCAATACCTTCTTCGTGTAGCAGAACCCTTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGATATCTTCCATTGGTTTCAGGGCCTGa                                                                                                                                          <  1:338310/149‑1 (MQ=255)
                 gCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCg                                                                                                                                 >  1:325339/1‑149 (MQ=255)
                 gCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCATGACCGTCACCg                                                                                                                                 >  2:315034/1‑149 (MQ=255)
                                 gtAGCCGAACCTTTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGTGTTATCTTTCATTGGTTTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCggg                                                                                                                <  1:162741/149‑1 (MQ=255)
                                 gtAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGAGATGTTTTGTCTGCCATGGGTTTAAGGGCCTGACCGTCACCGAAGTGAAAGGGTTCggg                                                                                                                >  1:241153/1‑149 (MQ=255)
                                           cATTACCTAATTCTGACCGGAGGCGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGaaa                                                                                                       <  2:91426/149‑1 (MQ=255)
                                              ctcccaaTTCTGACCCGAGGGGATCTATGAAGCTGGTGAACGTGATAATTAATCCATTCAAGCTGTAAGACGTTCGTGAAGCGTTATCGTTCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAggg                                                                                                    <  2:130162/144‑1 (MQ=255)
                                                              gAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAATCATTCAAGCTGGAAGACATTCGTGCAGCGTTATCTTCCATTGGTATTCAGGGCCTGTACGGCACCGAAGTGACAGGTTTAGGGCGTCAGAAAGGGCATGCCGAGGTGTAcc                                                                                    >  2:14382/1‑149 (MQ=255)
                                                                       tATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTACCTTCCATTGGTTTCAGGGCCTGACCGTCACCGAAGTGAAAGGGTTCGGGCGACAGAAAGCGTATGCCAAGCTGTGCCGGGGGGCGGa                                                                          >  1:56817/1‑149 (MQ=255)
                                                                                        tGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCATGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTc                                                          >  1:28140/1‑149 (MQ=255)
                                                                                        tGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTTTCAGGGCCTGACCGTCCCCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACACCGTCAATTTcc                                                         >  2:295390/1‑149 (MQ=255)
                                                                                        tGAGAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCTCCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTc                                                          <  1:226154/149‑1 (MQ=255)
                                                                                             aTCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGcc                                                     >  1:319986/1‑149 (MQ=255)
                                                                                             aTCAAACCATTCAAGCTGGAAGACGTTCGTCAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGGCAATTTCCTGcc                                                     >  1:168480/1‑149 (MQ=255)
                                                                                                 aaCCATTCAAGCTGGAAGACGTTCGTGCAGCGATAGCGTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGGATACAGCGTCAATTTCCTGCCaaaa                                                 >  2:68032/1‑149 (MQ=255)
                                                                                                                 aGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAATGTAAAAATTGATGTgg                                 <  1:203762/149‑1 (MQ=255)
                                                                                                                     gTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGTAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGAt                             <  1:14382/149‑1 (MQ=255)
                                                                                                                              gCGTTATCTGCCATCGGTTTCAGTGCCGGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGAc                   <  1:106916/149‑1 (MQ=255)
                                                                                                                                      ttCCATTGGTTGCAGGGCCTGACCTGCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGa           <  1:311672/149‑1 (MQ=255)
                                                                                                                                       tCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGa           <  2:319986/149‑1 (MQ=255)
                                                                                                                                                aTTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCTAAAAGAACATCAGCAAAATGCAAGTGCTGATGACCAACTCGATGAAGTGAt  <  1:244511/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TTTCTGGCACACCGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGCCAAAAGTAAAAATTGATGTGGCGATTGCTGATGACCAACTCGATGAAGTGAT  >  CP009273/467982‑468274

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