Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,719,504 G→A 9.5% P572P (CCC→CCT)  plsB ← Glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,719,5040GA 9.5% 122.6 / 4.1 42P572P (CCC→CCT) plsBGlycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (20/18);  new base A (2/2);  total (22/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-01

GTTGTTGCGGTAGTAGGTCATCAGCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGA  >  ACB122/1719372‑1719634
                                                                                                                                    |                                                                                                                                  
gttgttGCGGTAGTAGGTCATCAGCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGc                                                                                                                                 <  5:120565/136‑1 (MQ=255)
          tagtagGTCATCAGCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGgtgt                                                                                                                       <  5:13327/136‑1 (MQ=255)
            gtagGTCATCAGCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTgc                                                                                                                     <  4:8914/136‑1 (MQ=255)
            gtagGTCATCAGCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGGCCTGATGCTCGGCCACCAGGTTCAAGCGGCGCCCCTGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGTGGGGTGTgc                                                                                                                     >  8:78241/1‑136 (MQ=255)
               ggTCATCAGCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTgcggc                                                                                                                  >  6:72600/1‑136 (MQ=255)
                 tcatcaGCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGa                                                                                                                <  8:127479/136‑1 (MQ=255)
                       gCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGGTCGGCGAGTAggg                                                                                                          >  5:125332/1‑136 (MQ=255)
                             cgTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTg                                                                                                    <  5:72600/136‑1 (MQ=255)
                              gTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTgg                                                                                                   >  7:38422/1‑136 (MQ=255)
                                  ccTACATCCAGGTTGAGCATCGGGCCCCGGGCGGCCTGCTGCTCGGCCACCAGGTTAATGCGCCGCACTTGTTCGATCCGTGCCTGGCCGTCGCCGTCGGGCAGGGGGGTGTGCGGCGAGAAGGGCACTTGGCgca                                                                                               <  4:7487/136‑1 (MQ=255)
                                       aTCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGcagcag                                                                                           <  6:178633/135‑1 (MQ=255)
                                       aTCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCagcagc                                                                                          <  4:195949/136‑1 (MQ=255)
                                          cAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAgcg                                                                                         >  8:114194/1‑134 (MQ=255)
                                          cAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGcc                                                                                       >  2:75806/1‑136 (MQ=255)
                                          cAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGcc                                                                                       >  4:106782/1‑136 (MQ=255)
                                          cAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGGGTGCGGCGAGTAGGGCACTTGGCGCAGCAgcg                                                                                         <  7:114194/134‑1 (MQ=255)
                                               agagGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTa                                                                                  >  8:31971/1‑136 (MQ=255)
                                               agagGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTa                                                                                  >  7:19686/1‑136 (MQ=255)
                                                         gCCCAGTGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTTCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGt                                                                        <  3:38337/136‑1 (MQ=255)
                                                           ccAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGt                                                                             <  7:146601/129‑1 (MQ=255)
                                                           ccAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGt                                                                             >  8:146601/1‑129 (MQ=255)
                                                                    ccTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGt                                                                        >  1:73911/1‑125 (MQ=255)
                                                                    ccTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGt                                                                        <  2:73911/125‑1 (MQ=255)
                                                                           cTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGgcg                                                         <  1:75806/133‑1 (MQ=255)
                                                                                ccaccaGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTc                                                 >  3:79451/1‑136 (MQ=255)
                                                                                 caccaGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCgca                                                                                               <  8:114649/89‑1 (MQ=255)
                                                                                 caccaGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCgca                                                                                               >  7:114649/1‑89 (MQ=255)
                                                                                   ccaGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAg                                              >  5:123672/1‑136 (MQ=255)
                                                                                                gcgcACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGGGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCtg                                 <  8:37694/136‑1 (MQ=255)
                                                                                                 cgcACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCtgg                                >  7:185137/1‑136 (MQ=255)
                                                                                                  gcACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCtggt                               <  8:38422/136‑1 (MQ=255)
                                                                                                    acaTGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCtggtgg                             >  7:167771/1‑136 (MQ=255)
                                                                                                              tCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGc                   >  5:69666/1‑136 (MQ=255)
                                                                                                              tCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGc                   >  3:199571/1‑136 (MQ=255)
                                                                                                              tCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGAGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGc                   >  1:191192/1‑136 (MQ=255)
                                                                                                                   gCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAggg              <  8:19686/136‑1 (MQ=255)
                                                                                                                     cTGGCCGTCGCCCTCAGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCgcagca                                                                                            >  1:123691/1‑56 (MQ=255)
                                                                                                                     cTGGCCGTCGCCCTCAGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCgcagca                                                                                            <  2:126877/56‑1 (MQ=255)
                                                                                                                     cTGGCCGTCGCCCTCAGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCgcagca                                                                                            <  2:123691/56‑1 (MQ=255)
                                                                                                                     cTGGCCGTCGCCCTCAGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCgcagca                                                                                            >  1:126877/1‑56 (MQ=255)
                                                                                                                              gCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTg   <  8:185137/136‑1 (MQ=255)
                                                                                                                               cccTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCCGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTGCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGTCCACCAGGTTGa  >  3:61003/1‑136 (MQ=255)
                                                                                                                                    |                                                                                                                                  
GTTGTTGCGGTAGTAGGTCATCAGCACCGCGTTGCCTTCATCCAGGTAGAGGATGCGGCCCAGGGCGTCCTTCTGCTCGGCCACCAGGTTCATGCTGCGCACATGTTCGATCAGTGCCTGGCCGTCGCCCTCGGGCAGGGTGGTGTGCGGCGAGTAGGGCACTTGGCGCAGCAGCGCCAGGTACAGGTCGAGTACGCGGGTCAGGGCGCGCTCGTCCAGGGCAAGGCGGCTGGTGGACAGCAGTGCCAGGGCCACCAGGTTGA  >  ACB122/1719372‑1719634

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