Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,348,853 A→G 16.0% K257R (AAA→AGA)  gltB → glutamate synthase, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,348,8530AG16.0% 44.9 / 3.2 25K257R (AAA→AGA) gltBglutamate synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base A (10/11);  new base G (2/2);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.04e-01

GCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGT  >  CP009273/3348711‑3348989
                                                                                                                                              |                                                                                                                                        
gCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCa                                                                                                                                    >  1:83999/1‑149 (MQ=255)
          tGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAGATTCCAGCCACCGCtt                                                                                                                          >  2:112454/1‑149 (MQ=255)
          tGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGAGTGCGCGTACTTATAGATTCCAGACACCGCtt                                                                                                                          >  2:243559/1‑149 (MQ=255)
                accgGCGCGTCTCCCCTACCCCCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACTCCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCt                                                                                                                    <  1:128921/149‑1 (MQ=255)
                                   cacCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCAc                                                                                                 >  2:268522/1‑149 (MQ=255)
                                            gcgcTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCCCGACGCCGCACCGTTCGTCa                                                                                        >  2:73561/1‑149 (MQ=255)
                                             ccctcgccgccggcGCAACCGCTCCGCGACCTCGCGCATAACGGTCAAATCAACCCCATCACCGGTAACCGCCAATGGGCGCCCGCGCCCACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCaa                                                                                       <  1:8704/144‑1 (MQ=255)
                                                             aaCCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCAGCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACCCCCAACCCTTCGTCAACGAAACCGGCTCTGAc                                                                       >  1:47498/1‑149 (MQ=255)
                                                                cGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTcc                                                                    >  2:46428/1‑149 (MQ=255)
                                                                     cGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAAAGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCCTATCCCTGACTTGCACGACGCCGCACCGTTCGTCAACGAAACCAGCTCTGACTCCAGTTc                                                               >  2:249009/1‑149 (MQ=255)
                                                                       tgaTCTGGCGCATAACGGTGAAATCCACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCATTTcta                                                             <  1:237861/147‑3 (MQ=255)
                                                                       ccaccTGGCGTACAACGGTGAACTCAACACCACCACCGGCAACCGCCACTGGGCTCGTGCGCGTACTTCTAAATGCCAGCCATCGCTTATCCCTGTCCTGCACGACGCCGCACCGTTCGTCAACGAAACCTGCTCTGACTCCAGTTCGa                                                             <  2:115274/145‑1 (MQ=255)
                                                                                   tAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGc                                                 <  2:278655/149‑1 (MQ=255)
                                                                                        tggAAATCGACCCCATCACCGGTAACCGCCAATGGGCGCGGGCGCGTACTTATAAATTCCAGACACCGCTTATCCCAGAAGCGCGGTATTCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCCGTTCCATGGATAACATGCTGGaa                                            <  1:76605/147‑1 (MQ=255)
                                                                                        gTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGaa                                            <  2:177568/149‑1 (MQ=255)
                                                                                           aaaTCAACACCATCACCGGGAACCGCCAATGGGCGCGTTCGCGTACTTATAAATACCAGTGATCCCGTACCCCTGCCCTGAACGACGCCGCACCGTCCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAActg                                         <  2:83999/149‑1 (MQ=255)
                                                                                             aTCAACCCCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGCCTCCAGTTCGATGGATAACATGCTGGAActgct                                       <  2:55738/149‑1 (MQ=255)
                                                                                                aaCACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAGATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATTGATAACATGCTGGAActgctgct                                    <  1:243559/149‑1 (MQ=255)
                                                                                                aaCACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAGATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAActgctgct                                    <  1:112454/149‑1 (MQ=255)
                                                                                                                 ccGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGGTGGata                   >  1:23183/1‑149 (MQ=255)
                                                                                                                   gCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATatc                 >  1:92282/1‑149 (MQ=255)
                                                                                                                           gcgcgtgcgcgtACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAATTGCTGCTGGCAGGCGGGATGGATATCATCcgcgg         >  1:120400/1‑148 (MQ=255)
                                                                                                                           gcgcgtgcgcgtACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCgcgc         >  1:257648/1‑149 (MQ=255)
                                                                                                                                  ctcgtACTAATGAATTCCAGACACCGCTTATCCCGGACCTGCACGCCGCCGCACCGTTCGTCAACGACACCGGCTCTGGCTCCAGTTCGATGGATAACATGCTGTAACTGCGGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGt  <  1:335266/147‑1 (MQ=255)
                                                                                                                                  cgcgtACTTATAAATTCCAGGCACCGCTTATCCCTGAGCGGCACGACGCCGCACCGTGCGTCAACGGAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTAGAACTGCTGCTGGCAGGCGGGATGGATATCAGCCGCGCCATGCGt  <  1:249009/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                        
GCCATTTGCCTGTTCCACCAGCGCTTCTCCACTAACACCGTACCGCGCTGGCCGCTGGCGCAACCGTTCCGCTATCTGGCGCATAACGGTGAAATCAACACCATCACCGGTAACCGCCAATGGGCGCGTGCGCGTACTTATAAATTCCAGACACCGCTTATCCCTGACCTGCACGACGCCGCACCGTTCGTCAACGAAACCGGCTCTGACTCCAGTTCGATGGATAACATGCTGGAACTGCTGCTGGCAGGCGGGATGGATATCATCCGCGCCATGCGT  >  CP009273/3348711‑3348989

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: