Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,674,731 A→C 35.9% D65A (GAC→GCC)  recJ → Single‑stranded‑DNA‑specific exonuclease RecJ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,674,7310AC35.9% 19.1 / 5.8 37D65A (GAC→GCC) recJSingle‑stranded‑DNA‑specific exonuclease RecJ
Reads supporting (aligned to +/- strand):  ref base A (15/8);  new base C (7/7);  total (22/15)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATGG  >  ACB122/1674602‑1674866
                                                                                                                                 |                                                                                                                                       
tGACCCGCCCTTACGCCGCGCGCGGCGTGCAGCGCGAGGCCGAGCTGGACAAGCGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGGCGAAGCCCTCGCcccgc                                                                                                                                    >  3:378787/1‑135 (MQ=255)
                     gcggcgTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGGCGAAGCCCTCGACCAGCGCCCGCGCATCCTCCtcgtcg                                                                                                               >  1:144351/1‑135 (MQ=255)
                     gcggcgTGCAGAGCGAGGCCGAGCTGGACAAGAGCCAGGCACGCCTGCAGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCAGGTAGCAGCCATCGCCCCGCGCCCGCGCCTCCTCCtcggccg                                                                                                              >  7:395518/1‑132 (MQ=255)
                                     ggCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGAc                                                                                              >  7:68608/1‑136 (MQ=255)
                                        cGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATTCCGACgg                                                                                            >  6:360788/1‑135 (MQ=255)
                                        cGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGCCgg                                                                                            >  3:306701/1‑135 (MQ=255)
                                        cGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACgg                                                                                            >  1:252558/1‑135 (MQ=255)
                                        cGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGt                                                                                           >  8:304573/1‑136 (MQ=255)
                                         gAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGCCCCGCGCCAGAGCCTCCTCATCGTCGGCGACTTCGATTCCGACGGt                                                                                           >  2:98578/1‑135 (MQ=255)
                                         gAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGt                                                                                           >  2:117486/1‑135 (MQ=255)
                                                  acgggCCGGGCACGCCTGCTGCCCTACCCGCAGCTCAAGGGCATCGCGGCCGGTGGGGACCTGCGGGTGGAAGCCAGCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCa                                                                                 <  8:278796/132‑1 (MQ=255)
                                                  aaGCGCCTGGCACGCCTGCGGCCTTACCAGAAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCTACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGc                                                                                   <  2:328103/134‑1 (MQ=255)
                                                  aaGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGCCCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGCCGGTGCCACCGcc                                                                                  >  3:238943/1‑135 (MQ=255)
                                                   ccagCCTGGCCCGCCGGCTGCCCTCCCACCAGCTCACGGGCACGGCGGCGGGTGTGGACCGGCGGGTCGAAGCCAGCCACAAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAg                                                                                <  8:138720/134‑1 (MQ=255)
                                                         tgccACGCCTGCTGCCCTACCAGCAGCTCAAGGGCAGCGAGGCCGGGGGGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCg                                                                           <  4:227565/132‑1 (MQ=255)
                                                          ggCACGCCTGCTGCCCTACCAGCAGCTCATGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGCCCTGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGt                                                                          <  4:99101/135‑1 (MQ=255)
                                                                    ctgcCCTACCAGCCGCTCAAGGGCAGCGAGGCCGGGGTGGACCTGCTGGGCTAAGCCAGCGCCCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTgg                                                               <  8:241297/136‑1 (MQ=255)
                                                                    ctgcCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCCTCGCCCCGCGCCCGCGCCTCCGCATCGTCGGCGCCGTCGATGCCGCCGGGGCCCCCGCCAGCCCCGGCGGGGTgggg                                                                >  5:39071/1‑132 (MQ=255)
                                                                       cccTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCTCCAGCGCATCCTCATAGTCGGCGACTTCGATGCCGACGGTGCCACCGCCCGCCCCGTCGGGGTGCTGGgc                                                             >  4:334454/1‑135 (MQ=255)
                                                                              agcagcTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGCCCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGgcctgcgcct                                                     >  7:132790/1‑136 (MQ=255)
                                                                                   cgcAAGGGCTGCGTGGCCGGGGGGGACCTGCTGGTCGACCCCATCCCCCAGCGCCAGCGCATCCTCATCGGCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCctgctg                                                 <  2:211295/133‑1 (MQ=255)
                                                                                    tCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCCGCGCATCCTCATCGTCGGCGACTTCGATGCCGCCGGTGCCACCGCCCGCACCGGCGGGGTGCTGGGCCGGCGCCTGCTggg                                               >  7:363314/1‑136 (MQ=255)
                                                                                    tCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTTCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGACTGCTggg                                               >  7:278488/1‑136 (MQ=255)
                                                                                          gCATCGCGGCCGGGGGGGACCTGCTGGTCTAAGCCAGCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGc                                         <  8:278488/136‑1 (MQ=255)
                                                                                          gCAGCGTGGCGGGTGGGCACGTGCGGGGCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGGCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCgg                                          <  2:313332/135‑1 (MQ=255)
                                                                                           cATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCCGCGCCAGCGCATCCTCATCGTCGGCGACGTCGATGCCGACGGTGCCACCGCCCGCACCGTCGGGGTGCTGGGCCTGCGCCGGCTGGGTGCggg                                         >  6:26998/1‑134 (MQ=255)
                                                                                                aGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCCTg                                    >  1:313332/1‑135 (MQ=255)
                                                                                                              ctgctgGTCGAAGCCATCGACCCGCGCCAGCGCATCCTCATCGTCGGCGCCTTCGATGCCGACGGTGCCACCGCCACCACCGTCGGGGTGCGGGGCCTGCGCCTGCGGGGTGCGGCCCCTGTCGAGTACCTGGTg                      >  6:129325/1‑135 (MQ=255)
                                                                                                                ggggggcGAAGCCGTGGCCCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGccc                   <  8:395518/130‑1 (MQ=255)
                                                                                                                gctgGGCGACGCCATCGCCCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGcc                    <  5:377180/135‑1 (MQ=255)
                                                                                                                   ggTCGGAGCCCTCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGCTGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGGGCGGCCCAGGGCGATTACCTGGTGccgca                 >  2:69209/1‑132 (MQ=255)
                                                                                                                     tCGATGCCTTCGCCCCGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAAcc               <  3:208877/135‑1 (MQ=255)
                                                                                                                     gcGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGGCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAAcc               <  2:144351/134‑1 (MQ=255)
                                                                                                                           ccATCGCCCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTTGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGCTTACCGGGTGCCCCACCGCATCg         >  2:394562/1‑135 (MQ=255)
                                                                                                                            cgtCGCCCAGCGCCAGCGCAGCCGCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGGGCGGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGa        <  3:49659/133‑1 (MQ=255)
                                                                                                                            cATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGa        <  4:306701/135‑1 (MQ=255)
                                                                                                                            cATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGa        >  2:6590/1‑135 (MQ=255)
                                                                                                                                 aCCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCCACCGCTTCGAGGATg   >  1:146014/1‑135 (MQ=255)
                                                                                                                                 aCCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATgg  <  7:304573/136‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                       
TGACCCGCCTTTACGCCGCGCGCGGCGTGCAGAGCGAGGCCGAGCTGGACAAGAGCCTGGCACGCCTGCTGCCCTACCAGCAGCTCAAGGGCATCGAGGCCGGTGTGGACCTGCTGGTCGAAGCCATCGACCAGCGCCAGCGCATCCTCATCGTCGGCGACTTCGATGCCGACGGTGCCACCGCCAGCACCGTCGGGGTGCTGGGCCTGCGCCTGCTGGGTGCGGCCCATGTCGATTACCTGGTGCCCAACCGCTTCGAGTATGG  >  ACB122/1674602‑1674866

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