Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,066,951 C→T 9.7% S319S (TCG→TCA)  gatA ← Glutamyl‑tRNA(Gln) amidotransferase subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,066,9510CT 9.7% 123.0 / 3.5 41S319S (TCG→TCA) gatAGlutamyl‑tRNA(Gln) amidotransferase subunit A
Reads supporting (aligned to +/- strand):  ref base C (29/8);  new base T (2/2);  total (31/10)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

CGACAGGGCGTAGGTACCGACCATGATGCGCCGCTGTACTTCAACGCCGAACCCTTCGCCGCGGGAGCGCTTGTACAGGTCGGTGAGGTCTTTGGGCTCTTCACAGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGACCACTGCGCCGAGCTTTTCCAGCTCTTTGACGCT  >  ACB122/1066816‑1067070
                                                                                                                                       |                                                                                                                       
cgACAGGGCGTAGGTACCGACCATGATGCGCCGCTGTACTTCAACGCCGAACCCTTCGCCGCGGGAGCGCTTGTACAGGTCGGTGAGGTCTTTGGGCTCTTCACAGCGGTAGCCGAAGCGTACGCCGTCGAAAcgc                                                                                                                         >  4:46967/1‑136 (MQ=255)
         gTAGGTACCGACCATGATGCGCCGCTGTACTTCAACGCCGAACCCTTCGCCGCGGGAGCGCTTGTACAGGTCGGTGAGGTCTTTGGGCTCTTCACAGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTg                                                                                                                >  5:27342/1‑136 (MQ=255)
           aGGTACCGACCATGATGCGCCGCTGTACTTCAACGCCGAACCCTTCGCCGCGGGAGCGCTTGTACAGGTCGGTGAGGTCTTTGGGCTCTTCACAGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTgg                                                                                                               >  2:49892/1‑135 (MQ=255)
                   aCCATGATGCGCCGCTGTACTTCAACGCCGCACCCTTCGCCGCGGGAGCGGTTGTACAGGTCGGTGAGGTCTTTGGGCTCTTCACAGCGGTCGCCGAAGCGTACGCCGTCGAACCGCGACCGGTTGGAAGAGGCtt                                                                                                      >  3:92575/1‑136 (MQ=255)
                   aCCATGATGCGCCGCTGTACTTCAACGCCGAACCCTTCGCCGCGGGAGCGCTTGTACAGGTCGGTGAGGTCTTTGGGCTCTTCACAGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCtt                                                                                                      >  2:120474/1‑136 (MQ=255)
                    ccATGATGCGCCGCTGTACTTCAACGCCGAACCCTTCGCCGCGGGAGCGCTTGTACAGGTCGGTGAGGTCTTTGGGCTCTTCACAGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCtt                                                                                                      >  7:23041/1‑135 (MQ=255)
                           gcgcCGCTGTACTTCAACGCCGAACCCTTCGCCGCGGGAGCGCTTGTACAGGTCGGTGAGGTCTTTGGGCTCTTCACAGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCg                                                                                                >  3:113081/1‑134 (MQ=255)
                               cgcTGTACTTCAACGCCGAACCCTTCGCCGCGGGAGCGCTTGTACAGGTCGGTGAGGTCTTTGGGCTCTTCACAGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGa                                                                                          >  8:65459/1‑136 (MQ=255)
                                              ccGAACCCTTCGCCGCGGGAGCGCTTGTACAGGTCGGTGAGGTCTTTGGGCTCTTCACAGCGCTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCg                                                                           >  7:116394/1‑136 (MQ=255)
                                                 aaCCCTTCGCCGCGGGAGCGCTTGTACAGGTCGGTGAGGTCTTTGGGCTCTTCACAGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGaa                                                                        >  8:110164/1‑136 (MQ=255)
                                                        cgccggGGGAGCGCTTGTCCTGGGCGGTGAGGTGTGTGGGCTCTTCCCAGCGGTAGCCGAAGCGTACGCCGGCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCAt                                                                  <  6:27342/135‑1 (MQ=255)
                                                         gccgcGGGAGCGCTTGTACAGGTCGGTGAGGTCTTTGGGCTCTTCACAGCGGTAGCCGAAGCGTACGCCGTCGCACCGCGACAGGTTGGCAGAGGCGTCCGCCGGAGCGGTCACGTAGTAGGCCGGTATGGCAtgc                                                                >  2:97278/1‑136 (MQ=255)
                                                           cgcgGGAGCGCTTGTACAGGTCGGTGAGGTCTTTGGGCTCTTCACAGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCAtgctg                                                              >  7:197884/1‑136 (MQ=255)
                                                                         tACAGGTCGGTGAGGTCTTTGGGCTCTTCACAGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAgg                                                 <  2:10066/135‑1 (MQ=255)
                                                                         gcctggTCGGTGAGGTCTTTGGCATCTTCACAGCGGTCGCCTATGGGTCCGGCGTCGAACCGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCTATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGc                                                <  3:23170/132‑1 (MQ=255)
                                                                                cGGTGAGGGCTTTGGGCGCTTCACAGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTc                                         <  3:160400/136‑1 (MQ=255)
                                                                                 ggTGAGGTCTTTGGGCTCTTCACAGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTcc                                        >  2:62206/1‑136 (MQ=255)
                                                                                 ggTGAGGTCTTTGGGCTCTTCACAGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCGGATTTcc                                        >  5:97329/1‑136 (MQ=255)
                                                                                     aGGTCTTTGGGCTCTTCACAGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGa                                    >  4:103277/1‑136 (MQ=255)
                                                                                      ggTCTTTGGGCTCTTCACAGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGAc                                   >  8:163809/1‑136 (MQ=255)
                                                                                           ttGGGCTCTTCACAGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGACCACTg                              >  4:144158/1‑136 (MQ=255)
                                                                                             gggCTCTTCACAGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGACCACTgcg                            >  2:30641/1‑136 (MQ=255)
                                                                                             gggATCTTCCAAGCGGTAGCGGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGACCACTgcg                            <  4:127718/136‑1 (MQ=255)
                                                                                                    tCACAGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGACCACTGCGCCGAGCt                     >  4:133395/1‑136 (MQ=255)
                                                                                                     cacaGCGGTAGCCGAAGCGTACGCCGTCGCAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGACCACTGCGCCGAGCtt                    >  2:6376/1‑136 (MQ=255)
                                                                                                     cacaGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGACCACTGCGCCGAGCtt                    >  8:200757/1‑136 (MQ=255)
                                                                                                     cacaGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGACCACTGCGCCGAGCtt                    >  2:105605/1‑136 (MQ=255)
                                                                                                     cacaGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGACCACTGCGCCGAGCtt                    >  2:29293/1‑136 (MQ=255)
                                                                                                        aGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGACCACTGCGCCGAGCTTTTc                 <  1:30641/136‑1 (MQ=255)
                                                                                                                ccGAAGCGTACGCCGTCGAAACGTGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGACCACTGCGCCGAGCTTTTc                 <  4:37398/128‑1 (MQ=255)
                                                                                                                ccGAAGCGTACGCCGTCGAAACGTGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGACCACTGCGCCGAGCTTTTc                 >  3:37398/1‑128 (MQ=255)
                                                                                                                ccGAAGCGTACGCCGTCGAAACGTGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGACCACTGCGCCGAGCTTTTc                 >  3:37399/1‑128 (MQ=255)
                                                                                                                ccGAAGCGTACGCCGTCGAAACGTGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGACCACTGCGCCGAGCTTTTc                 <  4:37399/128‑1 (MQ=255)
                                                                                                                ccGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGCTCACGTAGTAGGCCGGAATGGCATGCTGCATCTTCGGCAGGCTGATTTCCATGACAACTGCGACGAGCTTTTCCAGCTCtt         >  3:187383/1‑136 (MQ=255)
                                                                                                                ccGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGACCACTGCGCCGAGCTTTTCCAGCTCtt         >  7:124306/1‑136 (MQ=255)
                                                                                                                  gAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGACCACTGCGCCGAGCTTTTCCAGCTCTTTg       <  8:124306/136‑1 (MQ=255)
                                                                                                                  gAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGAGCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGACCACTGCGCCGAGCTTTTCCAGCTCTTTg       >  5:52109/1‑136 (MQ=255)
                                                                                                                       gTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGACCACTGCGCCGAGCTTTTCCAGCTCTTTGACGCt  >  6:180328/1‑136 (MQ=255)
                                                                                                                       gTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGACCACTGCGCCGAGCTTTTCCAGCTCTTTGACGCt  >  8:187671/1‑136 (MQ=255)
                                                                                                                       gTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGACCACTGCGCCGAGCTTTTCCAGCTCTTTGACGCt  >  1:99851/1‑136 (MQ=255)
                                                                                                                                  aaaCGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGa                                                                                          <  6:180169/37‑1 (MQ=255)
                                                                                                                                  aaaCGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGa                                                                                          >  5:180169/1‑37 (MQ=255)
                                                                                                                                       |                                                                                                                       
CGACAGGGCGTAGGTACCGACCATGATGCGCCGCTGTACTTCAACGCCGAACCCTTCGCCGCGGGAGCGCTTGTACAGGTCGGTGAGGTCTTTGGGCTCTTCACAGCGGTAGCCGAAGCGTACGCCGTCGAAACGCGACAGGTTGGAAGAGGCTTCAGCCGGAGCGATCACGTAGTAGGCCGGAATGGCATGCTGCATGTTCGGCAGGCTGATTTCCTTGACCACTGCGCCGAGCTTTTCCAGCTCTTTGACGCT  >  ACB122/1066816‑1067070

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑