Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,350,594 G→A 19.1% T837T (ACG→ACA gltB → glutamate synthase, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,350,5940GA19.1% 36.0 / 4.6 21T837T (ACG→ACAgltBglutamate synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base G (6/11);  new base A (2/2);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.38e-01

CCACGCCTACAACCCGGACGTGGTGCGCACGCTGCAACAAGCGGTACAAAGCGGCGAGTACAGCGACTATCAGGAATACGCGAAGCTGGTTAATGAGCGTCCGGCAACCACGCTGCGCGATCTGCTGGCAATTACGCCGGGTGAAAACGCGGTCAACATTGCTGATGTTGAACCGGCAAGCGAACTGTTTAAACGCTTTGATACCGCCGCGATGTCTATCGGCGCGTTAAGCCCGGAAGCCCACGAGGCGCTGGCGGAAGCGATGAACAGCATCGGCGGTAATT  >  CP009273/3350459‑3350742
                                                                                                                                       |                                                                                                                                                    
ccACGCCTACAACCCGGCCGTGGGACGCACGCTGCAACAAGCGGTACACAGCGGCAAGAACAGCGACTATCAGGAATACGCGAAGCGGGTTAATCGGCGTCCGGCAACCACCCTGCGCCACCTGCTGGCAATTACGCCGGGTGa                                                                                                                                              >  1:110570/1‑144 (MQ=255)
ccACGCCTACAACCCGGACGTGGTGCGCACGCTGCAACAAACGGAACAAAGCAGCGAGAACAGCGACTATCAGGAATACGCGAAGCTGGTTAATGAGCGTCCGGCAACCACGCTGCGCGATCTGCTGGCAATTACACCGGGTGAAAAcg                                                                                                                                         >  2:20197/1‑149 (MQ=255)
         cAACCCGGACGTGGTGCGCACGCTGCAACAAGCGGTACAAAGCGGCGAGTACAGCGACTATCAGGAATACGCGAAGCTGGTTAATGAGCGTCCGGCAACCACGCTGCGCGATCTGCCGGCAATTACGCCGGGTGAAAACGCGGt                                                                                                                                     <  1:207406/144‑1 (MQ=255)
             ccGGACGTGGTGCGCACGCTGCAACAAGCGGTACAAAGCGGCGAGTACAGCGACTATCAGGAATACGCGAAGCTGGTTAATGAGCGTCCGGCAACCACGCTGCGCGATCTGCTGGCAATTACACCGGGTGAAAACGCGGTCAACATTGc                                                                                                                            >  2:278952/1‑149 (MQ=255)
               ggACGTGGTGCGCACGCTGCAACAAGCGGTACAAAGCGGCGAGTACAGCGACTATCAGGAATACGCGAAGCTGGTTAATGAGCGTCCGGCAACCACGCTGCGCGATCTGCTGGCAATTACGCCGGGTGAAAACGCGGTCAACAt                                                                                                                               <  1:172202/144‑1 (MQ=255)
                        gcgcACGCTGCAACAAGCGGTACAAAGCGGCGAGTACAGCGACTATCAGGAATACGCGAAGCTGGTTAATGAGCGTCCGGCAACCACGCTGCGCGATCTGCTGGCAATCCCGCCGTGTGAAAACGCGGTCAACATTGCTGATGt                                                                                                                      <  1:106571/144‑1 (MQ=255)
                                    acaaGCGGTACAAAGCGGCGAGTACAGCGACTATCAGGAATACGCGAAGCTGGTTAATGAGCGTCCGGCAACCACGCTGCGCGATCTGCTGGCAATTACGCCGGGTGAAAACGCGGTCAACATTGCTGATGTTGAACCGGCAAg                                                                                                          <  1:251179/144‑1 (MQ=255)
                                                           aCAGCGACTATCAGGAATACGCGAAGCTGGTTAATGAGCGTCCGGCAACCACGCTGCGCGATCTGCTGGCAATTACGCCGGGTGAAAACGCGGTCAACATTGCTGATGTTGAACCGGCAAGTGAACTGTTTAAACGCTTTGATAccgcc                                                                              >  2:186283/1‑149 (MQ=255)
                                                                     ccAGGAATACGCGAAGCTGGTTAATGAGCGTCCGGCAACCACGCTGCGCGATCTGCTGGCAATTACGCCGGGTGAAAACGCGGTCAACATTGCTGATGTTGAATTGGTAATTGAACTGTTTAAACGCTTTGATACCGCCGCGATGTCTa                                                                    <  2:83881/148‑1 (MQ=255)
                                                                           aTACGCGAAGCTGGTTAATGAGCGTCCGGCAACCACGCTGCGCGATCTGCTGGCAATTACGCCGGGTGAAAACGCGGTCAACATTGCTGATGTTGAACCGGCAAGCGAACTGTTTAAACGCTTTGATACCGCCGCGATGTCTAt                                                                   >  1:113313/1‑144 (MQ=255)
                                                                               gcgAAGCTGGTTAATTAGCGTCCGGCAACCACGCTGCGCGATCTGCTGGCAATTACACCGGGTGAAAACGCGGTCAACATTGCTGATGTTGAACCGGCAAGCGAACTGTTTAAACGCTTTGATACCGCCGCGATGTCTATCGgc                                                               <  1:255463/144‑1 (MQ=255)
                                                                                     cTGGTTAATGAGCGTCCGGCAACCACGCTGCGCGGTCTGCTGGCAATTACGCCGGGTGAAAACGCGGTCAACATTGCTGATGTTGAACCGGCAAGCGAACTGTTTAAACGCTTTGATACCGCCGCGATGTCTATCGGCGCGTTAAGccc                                                    >  2:50239/1‑149 (MQ=255)
                                                                                          tAATGAGCGCCCGGCAACCACGCCGCGCGATCTGCTGGCAATTACACCGGGTGAAAACGCGGTCAACATTGCTGATGTTGAACCGGCAAGCGAACTGTTTAAACGCTTTGATACCGCCGCGATGTCTATCGGCGCGTTAAGccc                                                    <  1:18555/144‑1 (MQ=255)
                                                                                             ggAGCCTCCGGCAACCCCGCTGCGCGATCTGCTGGCAATTACGCCGGGTGAAAACGCGGTCAACATTTCTGATGTTGAACCGGCAAGCGAACTGTTTAAACGCTTTGATACCGCCGCGATGTCTATCGGCGCGTTAAGCCCGGa                                                 <  1:5277/143‑1 (MQ=255)
                                                                                                               gCTGCGCGATCTGCTGGCAATTACGCCGGGTGAAAACGCGGTTAACATTGCTGATGTTGAACCGGCAAGCGAACTGTTTAAACGCTTTGATACCGCTGCGATGTTTATCGGTGCGTTAAGCCCGGAAGCCCACGAGGCGCTGGCGGAAg                          <  2:104057/149‑1 (MQ=255)
                                                                                                                     cgATCTGCTGGCAATTACGCCGGGTGAAAACGCGGTCAACATTGCTGATGTTGAACCGGCAAGCGAACTGTTTAAACGCTTTGATACCGCCGCGATGTCTATCGGCGCGTTAAGCCCGGAAGCCCACGAGGCGCTGGCGGAAGc                         <  1:46117/144‑1 (MQ=255)
                                                                                                                      gATCTGCTGGCAATTACGCCGGGTGCAAACGCGGTCAACATTGCTGATGTTGAACCGGCAAGCGAACTGTTTAAACGCTTTGATACCGCCGCGATGTCTATTGACGCGTTAATCCCGGAAGCCGACGAGGCGCTGGCGGAAGCGCTGaa                   >  2:57931/1‑149 (MQ=255)
                                                                                                                      gATCTGCTGGCAATTACGCCGGGTGAAAACGCGGTCAACATTGCTGATGTTGAACCGGCAAGCGAACTGTTTAAACGCTTTGATACCGCCGCGATGTCTATCGGCGCGTTAAGCCCGGAAGCCCACGAGGCGCTGGCGGAAGCg                        <  1:235214/144‑1 (MQ=255)
                                                                                                                      gATCGGCTGGCAATTACGCCGGGTGAAAACGCGGTCAACATTGCTGATGTTGAACCGGCAAGCGAACTGTTTAAACGCTTTGATACCGCCGCGATGTCTTTCGGGGCGTTAAGCCCGGAAGCCCACGAGGCGCTGGCGGAAGCGATGaa                   <  2:183989/149‑1 (MQ=255)
                                                                                                                                   ttACGCCGGGTGAAGACGCGGTCAACATTGCTGATGTTGAACCGGCAAGCGAACTGTTTAGCTGCTCTGATACCGGCGCGATGTCTATCGGCGCGTTAAGCCCGGAAGCCGACGAGGCGCTGGCGGAAGCGATGAACAGCATCGGCGGt      <  2:123332/149‑1 (MQ=255)
                                                                                                                                    tACGCCGGGTGAAAACGCGGTCAACATTGCTGATGTTGAACCGACAAACGAAATATTTAAACGCTTTGATACCGCCGCGATGTCTATCGGCGCGTTAAGCCCGGAAGCCCACGAGGCGCTGGCGGAAGCGATGAACAGCATCGGCGGTa     >  2:259629/1‑149 (MQ=255)
                                                                                                                                       gCCGGGTGAAAACGCGGTCAACATTGCTGATGTTGAACCGGCAAGCGAACTGTTTAAACGCTTTGATACCGCCGCGATGTCTATCGGCGCGTTAAGCCCGGAAGCCCACGAGGCGCTGGCGGAAGCGATGAACAGCATCGGCGGTAAtt  >  2:135156/1‑149 (MQ=255)
                                                                                                                                       |                                                                                                                                                    
CCACGCCTACAACCCGGACGTGGTGCGCACGCTGCAACAAGCGGTACAAAGCGGCGAGTACAGCGACTATCAGGAATACGCGAAGCTGGTTAATGAGCGTCCGGCAACCACGCTGCGCGATCTGCTGGCAATTACGCCGGGTGAAAACGCGGTCAACATTGCTGATGTTGAACCGGCAAGCGAACTGTTTAAACGCTTTGATACCGCCGCGATGTCTATCGGCGCGTTAAGCCCGGAAGCCCACGAGGCGCTGGCGGAAGCGATGAACAGCATCGGCGGTAATT  >  CP009273/3350459‑3350742

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: