Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009273 468,093 Δ1 bp 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.97.3% 159.3 / ‑2.3 37coding (39/339 nt)glnKnitrogen assimilation regulatory protein for GlnL, GlnE, and AmtB
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base . (15/21);  minor base T (1/0);  total (16/21)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

TTCTTAACTTCCTGCTCTCTTTCTCGTTTTTCATTTCTGGCACACCGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGC  >  CP009273/467949‑468222
                                                                                                                                               |                                                                                                                                  
ttcttAACTTCCTGCTCTCTTTCTCGTTTTTCATTTCTGGCACACCGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGAc                                                                                                                              <  1:249233/149‑1 (MQ=255)
          ccTGCTCTCTTTCTCGTTTTTCATTTCTGGCACACCGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAg                                                                                                                    <  1:424908/149‑1 (MQ=255)
           cTGCTCTCTTTCTCGTTTTTCATTTCTGGCACACCGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCg                                                                                                                         <  2:575854/143‑1 (MQ=255)
            tGCTCTCTTTCTCGTTTTTCATTTCTGGCACACCGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACggtggt                                                                                                                        >  2:475973/1‑138 (MQ=255)
                      ctcGTTTTTCATTTCTGGCACACCGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCCCGCTGAAGACGTCCGTGAAGCGTTAt                                                                                                              <  2:412354/143‑1 (MQ=255)
                      ctcGTTTTTCATTTCTGGCACACCGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCa                                                                                                        <  1:331218/149‑1 (MQ=255)
                        cGTTTTTCATTTTTGGCACACCGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCAtt                                                                                                      <  1:6190/149‑1 (MQ=255)
                        cGTTTTTCATTTCTGGCACACCGCTTGCAATACCTTCGTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCAtt                                                                                                      <  1:67346/149‑1 (MQ=255)
                         gTTTTTCATTTCTGGCACACCGCTTGCAATACCTTCTTCGTGTAGCAGAACCCTTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGGTAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCtt                                                                                                           >  2:287866/1‑143 (MQ=255)
                         gTTTTTCATTTCTGGCACACCGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTg                                                                                                     <  1:543056/149‑1 (MQ=255)
                             ttCATTTCTGGCACACCGCTTGCAATACCTTCTTCGTGTAGCAGCACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCAt                                                                                                       <  2:541894/143‑1 (MQ=255)
                              tCATTTCTGGCACACCGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCAtt                                                                                                      <  2:471561/143‑1 (MQ=255)
                                       gCACACCGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATATATGAAGCTGGTGACCGTGGTAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTGATCTTCCCTTGGTATTCAg                                                                                             >  2:343070/1‑143 (MQ=255)
                                             cGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTc                                                                                 <  1:276161/149‑1 (MQ=255)
                                              gCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGGGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAATCGTTATCTTCCATTGGTATTCAGGGCCTGa                                                                                      <  2:690582/143‑1 (MQ=255)
                                                   cAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATCCAGGGCCTGACCGGCACCGaa                                                                           >  1:763218/1‑149 (MQ=255)
                                                      tACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTg                                                                        >  1:729209/1‑149 (MQ=255)
                                                       aCCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAAACGTTCGTGAAGCGGTATCTTCCATTGAGATTCAGGGCCTGACGGTCGCCg                                                                             >  2:262915/1‑143 (MQ=255)
                                                            cttcGTGTAGCAGAACCATTACCGAATTCTTACCGGGGGGGATCTATGAAGCTGGTGACCGTGATAATCGAAGCATTCAAGCTGAAGCCGTTCTTGAAGCGTTATTTTCCGTTGGTATTCAGGGCCTGACCGCCCCCGAAGTGAAAGGt                                                                  >  1:64160/1‑149 (MQ=255)
                                                             ttcGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGtt                                                                 >  1:552194/1‑149 (MQ=255)
                                                                  caaGCAGAACCATTACGGAATTTTGACCGCAGGGGATCTATGAAGCTTGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGt                                                                  <  2:38778/141‑1 (MQ=255)
                                                                     gCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGt                                                         >  1:406595/1‑149 (MQ=255)
                                                                       aGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCTTGATAATCAAACCGTTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCgg                                                             >  2:463782/1‑143 (MQ=255)
                                                                                  cGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGTTGAAGACGTTCGTGAAGCGTTATCTTCCACTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAg                                                  >  2:523890/1‑143 (MQ=255)
                                                                                   gAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAgg                                                 <  2:701480/143‑1 (MQ=255)
                                                                                            ccGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGAGTCAGAAAGGGCATGCGGa                                        >  2:407200/1‑143 (MQ=255)
                                                                                              tgAGGCGATCTATGGCGCTGGTGCCCGTGATGATCAAACCAGTCAAGCTGAAGACGTGCGTGAAGCGTTATCTTCCATTTGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGACAGGGCATGCCGAGc                                      <  2:2008/142‑1 (MQ=255)
                                                                                                gtgggATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCgg                              <  1:535000/147‑1 (MQ=255)
                                                                                                  gggATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCGTTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCgggg                            >  1:268360/1‑149 (MQ=255)
                                                                                                    gATCTATGAAGCTTGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTAcc                                <  2:371621/143‑1 (MQ=255)
                                                                                                       cTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCggg                             <  2:492186/143‑1 (MQ=255)
                                                                                                                    gACCCTGATAATCAACCCATTCAAGCTGAAGACGTCCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGCAACCGAAAAGTAAGGTTTCGGGAGGAACAGAGGGCACGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCaa          <  1:188592/149‑1 (MQ=255)
                                                                                                                     aCCGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCAGTGGTATTCAGGGCCTGACCGTCACGGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGc               >  2:59429/1‑143 (MQ=255)
                                                                                                                       cGTGATAATCAAACCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGt             >  2:232703/1‑143 (MQ=255)
                                                                                                                         tGATAATCAAACCATTCCAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTcc     <  1:391367/149‑1 (MQ=255)
                                                                                                                          gATCATCAATTCATTCAAGCTGAAGACGTTCGGGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGACAGGTTTCGGGCGTCAGAAAGGGCATGCGGATCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCt    <  1:158896/149‑1 (MQ=255)
                                                                                                                                  aaCCATTCAAGCTGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGc  >  2:452272/1‑143 (MQ=255)
                                                                                                                                               |                                                                                                                                  
TTCTTAACTTCCTGCTCTCTTTCTCGTTTTTCATTTCTGGCACACCGCTTGCAATACCTTCTTCGTGTAGCAGAACCATTACCGAATTCTGACCGGAGGGGATCTATGAAGCTGGTGACCGTGATAATCAAACCATTCAAGCTGGAAGACGTTCGTGAAGCGTTATCTTCCATTGGTATTCAGGGCCTGACCGTCACCGAAGTGAAAGGTTTCGGGCGTCAGAAAGGGCATGCCGAGCTGTACCGGGGGGCGGAATACAGCGTCAATTTCCTGC  >  CP009273/467949‑468222

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