Predicted mutation
evidence seq id position mutation freq annotation gene description
RA KT2440 1,720,951 T→G 39.0% 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,9510TG39.0% 23.0 / 9.1 39D90A (GAC→GCC) plsBGlycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (15/8);  new base G (4/12);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 2.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGGGTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACGGCCAGGTCGGTGAGCGCGGGCGAGGGC  >  KT2440/1720817‑1721082
                                                                                                                                      |                                                                                                                                   
ccGGGGTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGt                                                                                                                                     <  3:248371/135‑1 (MQ=255)
   gggTCTGGCCCCAGTACTCGCTGAGAGGATTGTTCTGCGCTTCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGGGGGCGGGGCGCCTCTTTTGTCGTGGCGGCCGAGCCAGTCCGGGTCgg                                                                                                                                  <  8:167455/135‑1 (MQ=255)
   gggTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCgg                                                                                                                                  <  1:106000/135‑1 (MQ=255)
   gggTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCgg                                                                                                                                  <  8:264758/135‑1 (MQ=255)
   gggTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCgg                                                                                                                                  <  8:276579/135‑1 (MQ=255)
   gggTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCgg                                                                                                                                  <  8:332644/135‑1 (MQ=255)
   gggTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGGCgg                                                                                                                                  <  7:18031/135‑1 (MQ=255)
        tGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTGCTGGCGGCCGACCCCGGGCgggggggggggg                                                                                                                             >  1:168025/1‑133 (MQ=255)
                 aaCACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTGTAGTCCTGGCCGCCGCGGCCGGCCAGGTCGGGGGTAAgggtaga                                                                                                                     >  1:314344/1‑129 (MQ=255)
                  acacGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGGCGGGGGTCAGGTtgaaga                                                                                                                   >  5:94031/1‑135 (MQ=255)
                  acacGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCGGGCGGCCGCGCCAGTCCGGGTCGGGGGTCAGGTagaaga                                                                                                                   >  6:220998/1‑135 (MQ=255)
                       gtGACAGGAATGATCTGCGCATCTTCTGCGGCAGGCGGGCTGAGGGCAGCCACCCGGCGCTCGAGGGGGGGGGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGGGGGGGGTCTGGTAGAAGAAcgcg                                                                                                              <  6:94031/134‑1 (MQ=255)
                       cTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAAcgcg                                                                                                              >  8:276834/1‑135 (MQ=255)
                       cTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGGGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAAcgcg                                                                                                              <  8:270085/135‑1 (MQ=255)
                         gACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCCAGCCAGTCCGGGTCGGGGGGCAGGTATAACAAcgccgc                                                                                                            >  5:284735/1‑135 (MQ=255)
                         gACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCGTTTGTGCGGGCGGCCGCGCCCGGGCGCGTCGGGGGGCCGGTTGGAGACCgcggc                                                                                                            >  5:239163/1‑135 (MQ=255)
                          aCAGGAATGATCTGCGCATCTGCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGGGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGGCAGGTAGAAGAACgcggcg                                                                                                           <  8:41129/135‑1 (MQ=255)
                              gAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGtt                                                                                                       >  8:64972/1‑135 (MQ=255)
                                         gcATCTTCTTCGGCATGCTGGCTGACGGCAGCCCCCAGGCGCTCGAGCGTGGGCGGGGCGCCGCGTTTGTCCGGGCGGGCGCGCCAGTCCGGGGCGGGGGGGAGGTAGAAGCAACCGGCGGGTTGCTGCAGTggg                                                                                            >  2:135240/1‑135 (MQ=255)
                                          cATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTTTCCGTGCGGGCGCGGCCGGCCGGGGCCGGGGTGGGGGAGAAGACCGCGGCGGGTTCCTGCAGTGGGc                                                                                           >  1:156535/1‑135 (MQ=255)
                                           aTCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGCGCCAGGGCGGGTCGGGGGGCAGGTAGAAGAACGCGGCGGGGTCCGGCAGTGGGcc                                                                                          >  5:24759/1‑135 (MQ=255)
                                           aTCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGcc                                                                                          >  1:198134/1‑135 (MQ=255)
                                           aTCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGcc                                                                                          >  6:122557/1‑135 (MQ=255)
                                           aTCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCGTTTTTCCTGGCGGCCGAGCCAGTCCGGGGCGGGGGTGAGGGTGAAGAACGCGGCGGGGTCCTGCCGTggggc                                                                                          >  5:20473/1‑133 (MQ=255)
                                                 ttcGGCCGGCTGGCGGCCGGCAGCCACCCGGCGCTCGAGCGGGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGGCCGGGGCGGGGGGCAGGTAGACGACCGCGGCGGGTTCCTGCAGTGGGCCGACg                                                                                      <  4:71721/133‑1 (MQ=255)
                                                 ttcGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGGCCGGGCGGCCGAGCCAGTCCGGGGCGGGGGGCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACg                                                                                      <  8:11821/133‑1 (MQ=255)
                                                 ttcGGCATGCTGGCGGACGGCACCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACg                                                                                      <  4:309944/133‑1 (MQ=255)
                                                    ggCTTGCGGGCTGCCGGCGGCCACCCGGCGCTGGCGCGGGGGCGGGGGGCCGCTTTTGTCTGGGCGGCCGAGCCAGTCCGGGGCGGGGGTCGGTTAGAAGCCCGCGGCGGGTTCCTGCAGTGGGCCGACGGCTAc                                                                                 <  1:217952/135‑1 (MQ=255)
                                                     gCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCGTTTGTGCTGGGGGCCGAGCCAGTCCGGGTCGGGGGGGAGGTAGCATCACGCGGCGGGTTCCTGGAGTGGGGCGCCGGCGAcc                                                                                >  2:143418/1‑135 (MQ=255)
                                                     gCAGGCGGGCTCACGCCCCCCACCCGGCGCGGGCGCGGGGGCGGGGCGCCGCTTTTGCCCTGGCGGCCGAGCCAGTCGGGGGCGGGGGTCCGTTACAACAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTAcc                                                                                <  5:122557/135‑1 (MQ=255)
                                                            ggctgACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCGTTTGTCCTGGCGGGCGGGCCAGTGCGGGTGGGGGGGGAGGTTGAAGAACGCGGCGGGTTGCTGCAGTGGGCCGGCGGCTACCGgcagca                                                                         >  1:303014/1‑135 (MQ=255)
                                                            ggctgACGGCAGCCACAAGGCTCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCAGGGTCGGGGGTCAGGTATAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGGTACCGgccgca                                                                         >  6:18947/1‑135 (MQ=255)
                                                                                     gAGCGTGGGCGGGGCGCCGCTTTTGTCCTTGCGGCCGAGCCAGTCAGGGTCGGGGGTCAGGTAGAATAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGt                                                >  5:5329/1‑135 (MQ=255)
                                                                                         ggggggggggcgCGGTTTTGGCCCGGGGGGCCGCGCCAGGCGGGGGGGGGGGTCTGGTAGAAGACCGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCAc                                            <  5:220998/133‑1 (MQ=255)
                                                                                         gggggcggggcgCCGCTTTTGTCCTGGCGGCCGAGCCAGGCCGGGGCGGGGGGCAGGTAGAAGACCGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCAc                                            <  7:193594/133‑1 (MQ=255)
                                                                                                    gccgcTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGCCGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCGTGGTGCACTCGTGGTCGAg                                 >  1:110654/1‑135 (MQ=255)
                                                                                                     ccgcTTTTGTCCAGGCGGCCGCGCCAGTCCGGGTCGGGGGTCAGGTATAAGAAAGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCCGCAGCACCCGGCGCGGCCGCCCCCCCATGGGGCCCTCCGGGTcccgc                                >  6:254157/1‑131 (MQ=255)
                                                                                                        cTTTGGGCCGGGCGCCCGCGCCAGGCGGGGGGGGGGGGCTGGTACAACCCCGCGGCGGGTGCCGGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTGGGTGCACTCGTGGTCGAGCACg                             <  4:195008/135‑1 (MQ=255)
                                                                                                                cgggcggcCGAGCCAGGCCGGGGGGGGGGGCAGGTAGAAGACCGCGGCGGGGTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACGGCCAGGTc                     <  2:198134/133‑1 (MQ=255)
                                                                                                                     ggcCGCGCCTGGCGGGGGCGGGGGCCTGGTCGCACCCCGCGGCGGGTGCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGGCGAGCACGGCCAGGTCGGTGa                <  1:225273/135‑1 (MQ=255)
                                                                                                                                ccgggggCGGGGGCCAGGTACACCCCCGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACGGCCAGGTCGGTGAGCGCGGGCGAg     <  3:221288/132‑1 (MQ=255)
                                                                                                                                 ccGGGGCGGGGGGCAGGTCGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACGGCCAGGTCGGTGAGCGCGAGCGAgg    <  2:218722/135‑1 (MQ=255)
                                                                                                                                   ggggCGGGGGTCCGGGAGAAACCCGCGGCGGGTGCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTGGGTGCACTCGGGGGCGAGCACGGCCAGGTCGGTGAGCGCGGGCGAGGGc  <  3:151770/135‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                   
CCGGGGTCTGGCCCCAGAACACGCTGACAGGAATGATCTGCGCATCTTCTTCGGCATGCTGGCTGACGGCAGCCACCAGGCGCTCGAGCGTGGGCGGGGCGCCGCTTTTGTCCTGGCGGCCGAGCCAGTCCGGGTCGGGGGTCAGGTAGAAGAACGCGGCGGGTTCCTGCAGTGGGCCGACGGCTACCGGCAGCACCGGGCGCGGCAGCCCCGCCTTGGTGCACTCGTGGTCGAGCACGGCCAGGTCGGTGAGCGCGGGCGAGGGC  >  KT2440/1720817‑1721082

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: ‑