Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,616,780 G→A 8.5% G110G (GGC→GGT)  hutI_1 ← Imidazolonepropionase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,616,7800GA 8.5% 137.5 / 3.8 47G110G (GGC→GGT) hutI_1Imidazolonepropionase
Reads supporting (aligned to +/- strand):  ref base G (28/15);  new base A (2/2);  total (30/17)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.77e-01

ACGCGGGTGGGGAACCAGTCGTTCAGGCTCTGTCGATTGGTGATGACTCGGCTGCTTGCAGCGATACGGGGGCCTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGT  >  ACB122/3616649‑3616913
                                                                                                                                   |                                                                                                                                     
aCGCGGGTGGGGAACCAGTCGTTCATGGTGAGTCGATTGGTGATGACTCGGCTGCTTGCAGCGATACCGGGGCCCTAGAACAGTCCGGCGTTGACCCTATTTCGCCCGGCAATGGTGAGTTGGCCCCAGTCGCCg                                                                                                                                    >  6:91294/1‑135 (MQ=255)
aCGCGGGTGGGGAACCAGTCGTTCAGGCTCTGTCGATTGGTGATGACTCGGCTGCTTGCAGCGATACGGGGGCCTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCg                                                                                                                                    <  3:261630/135‑1 (MQ=255)
  gcgGGTGGGGAACCAGTCGTTCAGGCTCTGTCGATTGGTGATGCCTCGGCTGCTTGCAGCGATACGGGGGCCTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCggg                                                                                                                                  <  1:98545/135‑1 (MQ=255)
              ccAGTCGTTCAGGCTCTGTCGATTGGTGATGACTCGGCTGCTTGCAGCGATACGGGGGCCTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGa                                                                                                                      <  3:193532/135‑1 (MQ=255)
                   cGTTCAGGCTCTGTCGATTTGTGTTTACTCGGCTGCTTGAAGCGATACGGGGGCCTTCGAACAGTCCGGCGTTTCACGTATTTCGCACCGCCATTGTTAGTTTGCCACAGTCGCCGGGGGAGACCATGGTCGTCa                                                                                                                 >  6:27359/1‑135 (MQ=255)
                                   aTTGGTGATGACTCGGCTGCTTGCAGCGATACGGGGGCCTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAAc                                                                                                 >  7:202394/1‑135 (MQ=255)
                                              gtCGGCTGCTTGCAGCGATACGGGGGCCTTCGAACAGTCCGGAGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGTCAGCCAGAACGTGCTGAGCga                                                                                      <  1:29959/134‑1 (MQ=255)
                                              cTCGGCTGCTTGCAGCGATACGGGGGCCTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGGGCga                                                                                      <  3:156717/135‑1 (MQ=255)
                                                cGGCTGCTTGCAGCGATACGGGGGCCGTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCgaag                                                                                    <  7:247317/135‑1 (MQ=255)
                                                  gctgctTGCAGCGATACGGGCTCCTTCGAACAGTCCGGCGATGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCgaag                                                                                    <  7:240878/133‑1 (MQ=255)
                                                    tgctTGCAGCGATACGGGGGCCTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGCACGTGCTGAGCGCAGAACa                                                                                >  5:161091/1‑135 (MQ=255)
                                                    tgctTGCAGCGATACGGGGGCCTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACa                                                                                >  8:261585/1‑135 (MQ=255)
                                                    tgctTGCAGCGATACGGGGGCCTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACa                                                                                >  7:36684/1‑135 (MQ=255)
                                                         gCAGCGATACGGGGGCCTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCg                                                                                                                                    <  3:83803/78‑1 (MQ=255)
                                                         gCAGCGATACGGGGGCCTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCg                                                                                                                                    >  4:83803/1‑78 (MQ=255)
                                                           aGCGATACGGGGGCCTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGCACGTGCTGAGCGAAGAACAGCCCCCg                                                                         >  8:271556/1‑135 (MQ=255)
                                                           aGCGATACGGGGGCATTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCCCGCCGGCAGCCAGAACGTGATGCGCGAAGAACAGCCCCCg                                                                         >  5:217129/1‑135 (MQ=255)
                                                                          ttCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCg                                                                                                                                    <  3:244681/61‑1 (MQ=255)
                                                                          ttCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCg                                                                                                                                    >  4:244681/1‑61 (MQ=255)
                                                                             gAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCg                                                                                                                                    <  8:187714/58‑1 (MQ=255)
                                                                             gAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCg                                                                                                                                    >  7:187714/1‑58 (MQ=255)
                                                                                      ggCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGt                                                 >  7:117362/1‑132 (MQ=255)
                                                                                                gTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGCACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAgatgtcgatg                                    >  7:294077/1‑135 (MQ=255)
                                                                                                gTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGaaa                                                             >  5:272469/1‑110 (MQ=255)
                                                                                                gTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGaaa                                                             <  6:272469/110‑1 (MQ=255)
                                                                                                gTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAgatgtcgatg                                    >  5:296453/1‑135 (MQ=255)
                                                                                                gTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGAGCTGAGCGAAGAACAGCCCCCGCAGCGCGCGAAACTCGTCGCTGGTCCAgatgtcgagg                                    >  5:10157/1‑133 (MQ=255)
                                                                                                gTAATGCGCACGGCAATGCTGAGTTCGCCACTTTCGAATGTGAAGACCATGGACGAACCGCCTGCAGCCAGAACGTGCTAAGCTAAGAACAGACACCGCAGGGCGCGAAACTCCTAGCTGGTAAAgatctcca                                      >  8:268381/1‑131 (MQ=255)
                                                                                                    tGCGCACGGCAATGCTGAGTGGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCt                                 <  3:42348/134‑1 (MQ=255)
                                                                                                     gcgcACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTc                               >  5:197742/1‑135 (MQ=255)
                                                                                                      cgcACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCt                                 <  7:292952/132‑1 (MQ=255)
                                                                                                        cACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTc                            <  1:241657/135‑1 (MQ=255)
                                                                                                         aCGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCg                           >  6:158430/1‑135 (MQ=255)
                                                                                                           ggCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGCGCGAAGCACCGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTACAGCTGTCGATGTCTTCCTCGct                         >  3:116359/1‑135 (MQ=255)
                                                                                                           ggCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGct                         >  3:71622/1‑135 (MQ=255)
                                                                                                           ggCAATGCTGAGTTGGCCACTGTCACCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGccccc                                                                          >  5:215414/1‑86 (MQ=255)
                                                                                                           ggCAATGCTGAGTTGGCCACTGTCACCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGccccc                                                                          >  5:211469/1‑86 (MQ=255)
                                                                                                           ggCAATGCTGAGTTGGCCACTGTCACCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGccccc                                                                          <  6:215414/86‑1 (MQ=255)
                                                                                                           ggCAATGCTGAGTTGGCCACTGTCACCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGccccc                                                                          <  6:211469/86‑1 (MQ=255)
                                                                                                              aaTGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTg                      >  2:177557/1‑135 (MQ=255)
                                                                                                               aTGCTGTGTTGGCCAGTGTCGCCGGGGCCGACCATGGACGACACGCCGGCAGCCAGAACGTGCTGAGCGATGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTgg                     <  5:274747/135‑1 (MQ=255)
                                                                                                                        tGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGcc             >  5:167446/1‑134 (MQ=255)
                                                                                                                        tGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATCGcc             >  2:205169/1‑134 (MQ=255)
                                                                                                                          gCCACTGTCGCCGGGGCAGACCATGGCCGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTGGCTCTGGGCATTGCCaaa          >  6:95546/1‑135 (MQ=255)
                                                                                                                          gCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCaaa          >  6:87333/1‑135 (MQ=255)
                                                                                                                                  cGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAggt  >  3:136763/1‑135 (MQ=255)
                                                                                                                                  cGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAggt  >  1:53122/1‑135 (MQ=255)
                                                                                                                                   |                                                                                                                                     
ACGCGGGTGGGGAACCAGTCGTTCAGGCTCTGTCGATTGGTGATGACTCGGCTGCTTGCAGCGATACGGGGGCCTTCGAACAGTCCGGCGTTGACCGTATTGCGCACGGCAATGCTGAGTTGGCCACTGTCGCCGGGGCAGACCATGGACGTCACGCCGGCAGCCAGAACGTGCTGAGCGAAGAACAGCCCCCGCAGGGCGCGAAACTCGTCGCTGGTCCAGATGTCGATGTCTTCCTCGCTCTGGGCATTGCCAAACGCCAGGT  >  ACB122/3616649‑3616913

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