Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 504,483 G→A 8.0% R136Q (CGG→CAG)  trpC → Indole‑3‑glycerol phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB122504,4830GA 8.0% 152.9 / 3.9 50R136Q (CGG→CAG) trpCIndole‑3‑glycerol phosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (22/24);  new base A (2/2);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

GGTGTTGACCGATGTCGATTATTTCCAGGGTGCCGATGAGTACTTGCAGCAGGCCCGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCTGAACTGGCCGCCACTGCCAAGGACGTCGGCCTCGACGTGCTGGTGGAAGTGCACGATGGCGATGAAC  >  ACB122/504358‑504611
                                                                                                                             |                                                                                                                                
ggTGTTGACCGATGTCGATTATTTCCCGGGTTCCGATGAGTACTTGCAGCAGGCCCGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTggg                                                                                                                         <  5:247667/135‑1 (MQ=255)
ggTGTTGACCGATGTCGATTATTTCCAGGGTGCAGATGAGTACTTGCAGCAGGCCCGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTggg                                                                                                                         <  8:143130/135‑1 (MQ=255)
   gtTGACCGATGTCGATTATTTCCAGGGTGCCGATGAGTACTTGCAGCAGGCCCGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGc                                                                                                                      <  5:124593/135‑1 (MQ=255)
   gtTGACCGATGTCGATTATTTCCAGGGTGCCGATGAGTACTTGCAGCAGGCCCGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGc                                                                                                                      <  3:166201/135‑1 (MQ=255)
   gtTGACCGATGTCGAGTATTTCCAGGGTGCCGATGAGTACTGGCAGCAGGCCCGCGCGGCGGTTTCGCGGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGc                                                                                                                      <  8:236866/135‑1 (MQ=255)
   ggtGACCGATGTCGATTATTTCCAGGGTGCCGATGAGTACTTGCAGCAGGCCCGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGTTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGc                                                                                                                      <  2:72273/133‑1 (MQ=255)
       aCCGATGTCGATTATTTCCAGGGTGCCGATGAGTACTTGCAGCAGGCCCGCGCGGCGGTTTCGCTGCCGGTGCTCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGATGCCCGGGCCCTGGGGGCAGAt                                                                                                                  >  7:54191/1‑135 (MQ=255)
       aCCGATGTCGATTATTTCCAGGGTGCCGATGAGTACTTGCAGCAGGCCCGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGCTGAt                                                                                                                  >  2:153209/1‑135 (MQ=255)
        ccGTTGTCGATTATTTCCAGGGTGCCGATGAGTACTTGCAGCAGGCCCGCGCGGCGGTTTCGCTGCCGGTGACCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGCAGAtt                                                                                                                 <  2:111622/135‑1 (MQ=255)
               cGATTATTTCCAGGGTGCCGATGAGTACTTGCAGCAGGCCCGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGCAGAtt                                                                                                                 <  5:275541/128‑1 (MQ=255)
               cGATTATTTCCAGGGTGCCGATGAGTACTTGCAGCAGGCCCGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGCAGAtt                                                                                                                 >  6:275541/1‑128 (MQ=255)
                  ttattTCCAGGGTGCCGATGAGTACTTGCAGCAGGCCCGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGCAGATTGCGTActgct                                                                                                       <  6:358732/135‑1 (MQ=255)
                        ccAGGGTGCCGATGAGTACTTGCAGCAGGCCCGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCg                                                                                                  <  6:177148/134‑1 (MQ=255)
                               gCCGATGAGTACTTGCAGCAGGCCCGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGgcg                                                                                          >  5:235363/1‑135 (MQ=255)
                                 cGATGAGTACTTGCAGCAGGCGCGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGTCAGATTGCGTACTGCTGATCGTGTCGGCGCt                                                                                        <  4:181462/135‑1 (MQ=255)
                                           ttGCAGCAGGCCCGCGCGGCGGTTTCGCTGCTGGTGATCCTCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGt                                                                               <  5:331169/134‑1 (MQ=255)
                                           ttGCAGCAGGCCCGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCAGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGt                                                                               >  5:359363/1‑134 (MQ=255)
                                           ttGCAGCAGGCCCGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCAGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGt                                                                               >  5:353807/1‑134 (MQ=255)
                                           ttGCAGCAGGCCCGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCAGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGt                                                                               <  6:353807/134‑1 (MQ=255)
                                           ttGCAGCAGGCCCGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCAGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGt                                                                               <  6:359363/134‑1 (MQ=255)
                                                   ggCCCGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTggg                                                                                                                         >  2:366311/1‑84 (MQ=255)
                                                   ggCCCGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTggg                                                                                                                         <  1:366311/84‑1 (MQ=255)
                                                    gCCCGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCt                                                                     <  7:334763/135‑1 (MQ=255)
                                                     cccGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGc                                                                                                                      <  7:150886/85‑1 (MQ=255)
                                                     cccGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCt                                                                     >  8:334763/1‑134 (MQ=255)
                                                     cccGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGGCGAAGCCCGGGCCCTGGGGGc                                                                                                                      >  8:150886/1‑85 (MQ=255)
                                                         cgcgGCGGTTTCGCTGCCGGTGATCCGAAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCTGAACt                                                                <  3:145103/135‑1 (MQ=255)
                                                            ggcggTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGCTCGTGTCGGCGCTGGATGCCGTGAAGATGGCTGAACTGgc                                                             >  6:336125/1‑135 (MQ=255)
                                                                           ggTGAGCCGCAAGGACTTCATGGTCGACCCTTACCAGAGCGGCGAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCTGAACTGGCCGCCACTGCCAAGGa                                              <  6:178270/135‑1 (MQ=255)
                                                                              gATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGAAGTGAAGATGGCTGAACTGGCCGCCACTGCCAAGGCCGt                                           >  5:256482/1‑135 (MQ=255)
                                                                                   gCAAGGACTTCATGGGCGACCCTTACCAGAGCGTCGAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCTGAACTGGCCGCCACTGCCAAGGACGTCGGcc                                      <  8:152247/135‑1 (MQ=255)
                                                                                    cAAGTAATTCATGGTCGACCCTTACCAGATCGGCGAAGACCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTTTCGGCGCTTGATGACGTTAAGATGGCTGAACTGGCCGCCACTGCCAAGGACGTCGGCCt                                     <  6:281709/135‑1 (MQ=255)
                                                                                     aaGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTGGGCGCGGGATGACGTGAAGATGGCTGAACTGGCCGCCACTGCCAAGGACGGCGGCCTc                                    <  7:141773/135‑1 (MQ=255)
                                                                                      aGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCTGAACTGGCCGCCACTGCCAAGGACGTCGGCCTCg                                   >  4:3525/1‑135 (MQ=255)
                                                                                               tGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCTGAACTGGCCGCCACTGCCAAGGACGTCGGCCTCGACGTGCtgg                          >  3:311381/1‑135 (MQ=255)
                                                                                                ggTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCTGAACTGGCCGCCACTGCCAAGGACGTCGGCCTCGACGTGCtggt                         >  4:23128/1‑135 (MQ=255)
                                                                                                        cTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTTTCGGCGCTGGAGGACGTGAAGATGGCTGAACTGGCCGCCACTGCCAAGGACGTCGGCCTCGACGTGCTGGTGGAAGTGc                 >  1:365730/1‑135 (MQ=255)
                                                                                                        cTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCTGAACTGGCCGCCACTGCCAAGGACGTCGGCCTCGACGTGCTGGTGGAAGTGc                 >  1:140443/1‑135 (MQ=255)
                                                                                                              aGATCGTCGAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCTGAACTGgccg                                                           <  6:19734/87‑1 (MQ=255)
                                                                                                              aGATCGTCGAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCTGAACTGgccg                                                           >  5:19734/1‑87 (MQ=255)
                                                                                                              aGATCGTCGAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCTGAACTGGCCGCCACTGCCAAGGACGTCGGCCTCGACGTGCTGGTGGAAGTGCACGATg           <  8:49930/135‑1 (MQ=255)
                                                                                                                atggtcgAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCTGAACTGGCCGCCACTGCCAAGGACGTCGGCCTCGACGTGCTGGTGGAAGTGCACGATGGc         <  5:231609/132‑1 (MQ=255)
                                                                                                                   gtcgAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCTGAACTGGCCGCCACTGCCAAGGACGTCGGCCTCGACGTGCTGGTGGAAGTGCACGATGGCGa       >  3:130344/1‑134 (MQ=255)
                                                                                                                   gtcgAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCTGAACTGGCCGCCACTGCCAAGGACGTCGGCCTCGACGTGCTGGTGGAAGTGCACGATGGCGa       >  2:147407/1‑134 (MQ=255)
                                                                                                                   gtcgAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCTGAACTGGCCGCCACTGCCAAGGACGTCGGCCTCGACGTGCTGGTGGAAGTGCACGATGGCGa       >  1:232212/1‑134 (MQ=255)
                                                                                                                   gtcgAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCTGAACTGGCCGCCACTGCCAAGGACGTCGGCCTCGACGTGCTGGTGGAAGTGCACGATGGCGAt      >  4:172786/1‑135 (MQ=255)
                                                                                                                   gtcgAAGCCCGGCCCCTGGGGGCAGATTGCGTACTGCTGATCGTGGCGGCGCTGGATGACGTGAAGATGGCTGAACTGGCCGCCACTGCCAAGGACGTCGGCCTCGCCGTGCGGGTGGAAGTGCCCGATGGCGAt      >  2:163109/1‑135 (MQ=255)
                                                                                                                       aaGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCTGAACTGgccg                                                           <  2:60000/78‑1 (MQ=255)
                                                                                                                       aaGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCTGAACTGgccg                                                           >  1:60000/1‑78 (MQ=255)
                                                                                                                       aaGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCTGAACTGGCCGCCACTGCCAAGGACGTCTGCCTCGACGTGCTGGTGGAAGTGCACGATGGCGATGAAc  >  8:264542/1‑135 (MQ=255)
                                                                                                                             |                                                                                                                                
GGTGTTGACCGATGTCGATTATTTCCAGGGTGCCGATGAGTACTTGCAGCAGGCCCGCGCGGCGGTTTCGCTGCCGGTGATCCGCAAGGACTTCATGGTCGACCCTTACCAGATCGTCGAAGCCCGGGCCCTGGGGGCAGATTGCGTACTGCTGATCGTGTCGGCGCTGGATGACGTGAAGATGGCTGAACTGGCCGCCACTGCCAAGGACGTCGGCCTCGACGTGCTGGTGGAAGTGCACGATGGCGATGAAC  >  ACB122/504358‑504611

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