Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 1,720,951 T→G 56.1% D90A (GAC→GCC)  plsB ← Glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24401,720,9510TG56.1% 1.8 / 14.4 26D90A (GAC→GCC) plsBGlycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (8/3);  new base G (3/12);  total (11/15)
Fisher's exact test for biased strand distribution p-value = 1.49e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.87e-01

GGGGTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACGGCCAGGTCGGTGAGCGCGGGCGAGG  >  KT2440/1720819‑1721080
                                                                                                                                    |                                                                                                                                 
ggggTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCgg                                                                                                                                <  7:211219/136‑1 (MQ=255)
 gggTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCgg                                                                                                                                <  5:203642/135‑1 (MQ=255)
  ggTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGTCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCggg                                                                                                                               <  3:155874/135‑1 (MQ=255)
          cccAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCGTTTGTCCTGGCGGCCGAGCCAGTCCGGGGCGGGGGTCAggg                                                                                                                       >  1:188909/1‑134 (MQ=255)
                acacGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGGCAGGTAgaagaa                                                                                                                >  8:190058/1‑136 (MQ=255)
                     ctcaccggATTGATCTGCACATGTTATTCGGCATGCTGGCTGACGGCCGCCACCCGGCGCTCGCGCGTGGGGGGGGCGCCGCTTTTGTCCGGGGGGCCGAGCCAGTCCGGGGCGGGGGGCAGGTATAACAAcgcg                                                                                                            <  4:17390/129‑1 (MQ=255)
                     cTGACAGGAATGATCTGCGCATCTTCTTCGGCTTGGTGGCTGACGGCAGCAACCAGGCGCTCGAGCGGGGGGGGGGCGCCGTTTTTGTCCTGGCGGCCGAGCCAGTCGGGGGGGGGGGTCAGTTAGAAGAAcgcg                                                                                                            <  4:18784/135‑1 (MQ=255)
                     cTGACAGGAATGAGCTGCGCATCTTCTGCGGCATGCTGGCGGCCGGCAGCCACCTGGCGCTCGAGCGGGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGGCGGGGGTCAGGTAGAAGAAcgcg                                                                                                            <  4:251879/135‑1 (MQ=255)
                       gACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCACTCCCGGTCGGGGGTCAGGTAGAAGAACgcggc                                                                                                          >  4:258032/1‑135 (MQ=255)
                                 aTCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGc                                                                                                >  4:239520/1‑135 (MQ=255)
                                     gcgcATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCTGTCCGGGGCGGGGGTCAGGTAGAAGAAcgcg                                                                                                            <  6:108830/119‑1 (MQ=255)
                                     gcgcATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGGCGGGGGGCAGGTAGAAGAAcgcg                                                                                                            >  5:108830/1‑119 (MQ=255)
                                       gcATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTTTCCTGGCGGCCGAGCCACTCCGGGGAGCGGGGCAGGGATAAGAAACCCGGGGGGTACTGCAGgggg                                                                                          >  2:128366/1‑135 (MQ=255)
                                       gcATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCCGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGc                                                                                         >  8:65512/1‑136 (MQ=255)
                                             tcttcGGCATGCGGGCGGCCGGCAGCCACCAGGCGCTCGAGGGTGGGCGGGGCGCCGCTTTTGTCCGGGCGGCGGAGCCAGGCCGGGGCGGGGGACAGGTACAAGACCGCGGCGGGTTCCTGCAGTGGGCCGACg                                                                                    <  7:65512/135‑1 (MQ=255)
                                                   gCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTAcc                                                                              >  1:5850/1‑135 (MQ=255)
                                                                 ggCAGCCACCCGGCGCGCGAGCGGGGGGGGGGCGCCGTTTTTGGCGTGGCGGCGGAGCCGGTCGGGTGCGGGGGGCTGGTAGAAGACCGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGgcg                                                                <  1:58624/135‑1 (MQ=255)
                                                                     gCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGc                                                            >  4:105889/1‑135 (MQ=255)
                                                                      cccccaGGCGCGCGCGCGGGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCTGGCCGGGGCGGGGGGCAGGTAGACCACCGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCa                                                           <  4:183134/132‑1 (MQ=255)
                                                                      ccacccGGCGCGGGCGGGGGGGGGGGGCGCCGCTTTTGGCCGGGCGGCCGCGCCAGGCCGGGGCGGGGGGCTGGTAGAAGACCGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCa                                                           <  3:199321/135‑1 (MQ=255)
                                                                      ccaccaGGCGCTCGAGCGGGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGGCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCGGACGGCTACCGGCAGCACCGGGCGCGGCa                                                           <  4:95473/135‑1 (MQ=255)
                                                                                          ggcggggcgCCGCTGTTGGCCTGGCGGCCGAGCCAGTCCGGGGCGGGGGGCAGGTAGACGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGt                                      <  8:13086/136‑1 (MQ=255)
                                                                                                  gccgcTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCCGCCCCGCCTTGGTGCCCTCGGGGTCGAg                               >  6:40144/1‑135 (MQ=255)
                                                                                                                   ggcCGAGCCGGGGCGGGGGGGGGGTCCGTTATATCACCGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCGTGGTGCACTCGTGGTCGAGCACGGCCAGGTCGGTGa              <  3:50876/135‑1 (MQ=255)
                                                                                                                    gcCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACGGCCAGGTCGGTGAg             >  3:50687/1‑135 (MQ=255)
                                                                                                                        aGCCAGGCCGGGGCGGGGGGCAGGTAGAAGCCCGCGGCGGGTGCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACGGCCAGGTCGGTGAGcgcg         <  6:59936/135‑1 (MQ=255)
                                                                                                                               cgggggCGGGGGGCCGGTAGAAGACCGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACGGCCAGGGCGGTGAGCGCGGGCGAgg  <  3:214711/133‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                 
GGGGTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACGGCCAGGTCGGTGAGCGCGGGCGAGG  >  KT2440/1720819‑1721080

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