Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,301,325 A→C 49.1% intergenic (‑106/+3467) recC ← / ← PP_RS08750 RecBCD enzyme subunit RecC/tRNA‑Val

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,301,3250AC49.1% 1.9 / 3.4 20intergenic (‑106/+3467)recC/PP_RS08750RecBCD enzyme subunit RecC/tRNA‑Val
Reads supporting (aligned to +/- strand):  ref base A (8/2);  new base C (8/2);  total (16/4)
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.00e+00

GGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGC‑GGCCGTGCCC‑GCTCCCG‑CAGCGATTCGC‑G‑CAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGG  >  ACB122/5301198‑5301455
                                                                                                                                  |                                                                                                                                    
ggAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGGGC‑GGCCCGGGCC‑CCGCCCG‑Ccgcggata                                                                                                                               >  7:212152/1‑131 (MQ=255)
ggAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCCGGGGA‑CGCCGGGCCC‑CCCCCCG‑CCGCGaata                                                                                                                               >  1:208820/1‑132 (MQ=255)
     tGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGC‑GGCCGTGCCC‑GCTCCCG‑CAGCGATTCGC‑G‑Ca                                                                                                                        >  8:304522/1‑135 (MQ=255)
                      cGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGGGCGGGCCGGGCCC‑‑CTCCCC‑CCGCGCTTGGC‑GCCGGGGGGGGTTTTTTGcc                                                                                                        >  8:260091/1‑135 (MQ=255)
                         aaGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGGGC‑GGCCCGGCCC‑GCTCCCC‑CCGCCCGTCGC‑G‑CAGGTCGAGGTTTTTGGCcac                                                                                                     >  7:240771/1‑134 (MQ=255)
                                            aTGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGGGC‑CGCCGGGGCC‑CCGCCCC‑CCGCGCTTGGCGG‑CGGGGCCGTGTTTTTGCCC‑CAAAAACAAAACCCCCAccag                                                                                 >  4:311300/1‑133 (MQ=255)
                                                    gagGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGGGC‑CGGCGGGGCC‑CCGCCCC‑CCGCGGATAGCGG‑CGGGCGGGGTTTTTTGCCC‑CCAACACAAAACACCCACCAGGGGgcgcg                                                                         >  7:20955/1‑134 (MQ=255)
                                                    gagGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGCGC‑GGGCGGGGCC‑GCGCCCG‑CCGCGCTTGCC‑G‑CCCGGGGCGGTTTTTGGCCACAAAAACAAAAACCCAACCagcgtgc                                                                            >  2:258673/1‑132 (MQ=255)
                                                                       gTACATGAAACGGCCGAATGCCATCAGGCTGCGGGGGC‑GGCCGGGGCC‑CCGCCCC‑CCCCGCTTGGC‑G‑CAGGGGGAGGTTTTTTGCCACAAAGACAAAAACCCCACCAGCGGGCGCCCATAGGGGGTTTGGGCAAt                                                      >  2:237849/1‑135 (MQ=255)
                                                                         aCATGAAACGGCCGAATGCCATCAGGCTGCGGGCGC‑GGCCGTGCCC‑GCTCCCG‑CAGCGCTTCGC‑G‑CAGGTCGACTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCCCAGAAAGGGGTTTTGGGAAAtt                                                     >  8:317486/1‑133 (MQ=25)
                                                                           aTGAAACGGCCGAATGCCATCAGGCTGCGGGGGC‑GGGCCTTCCC‑CCGCCCC‑CCGCGCATACC‑G‑CAGGGCGAGGTTTTTTGCCACAAAAACAAAAACCCAAACAACGGGCGCCCACAGGGGTTTTTGGaaaataa                                                   >  6:199886/1‑130 (MQ=255)
                                                                                                      gggtggc‑gggcggGCCC‑GCGCCCGACCGCGATTCGC‑G‑CAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATgg                        <  8:330044/130‑1 (MQ=14)
                                                                                                          gccgggCCGGGCCCGGCCCCCG‑CCGCGATTCGC‑G‑CAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGa                     <  1:207496/131‑1 (MQ=14)
                                                                                                              ggCCGTGCCC‑GCTCCCG‑CAGCGATTCGC‑G‑CAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAAccc                <  8:129861/135‑1 (MQ=255)
                                                                                                                  gTGCCC‑GCTCCCG‑CAGCGATTAGC‑G‑CAGGTGGAGTTTTTTTCCCACAAAGACAAAAACCCAACCAGCGTGCGGAGGTTGGGGGTTTGCGAAATGAATGTTGTAGATGACCTACTCTCACAAGGGGAAAccccccc            >  2:295903/1‑132 (MQ=255)
                                                                                                                    gccc‑gcccCCG‑CAGCGATTCGC‑G‑CAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATgggg                      >  6:375220/1‑123 (MQ=255)
                                                                                                                    gCCC‑GCTCCCG‑CAGCGATTCGC‑G‑CAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAAAGAATCTTGACGATGCCCTACACGCACATGGGGCAACCCCCCACt          >  7:248739/1‑135 (MQ=17)
                                                                                                                    gCCC‑CCCCCCG‑CAGCGATTCGC‑G‑CAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATgggg                      <  5:375220/123‑1 (MQ=255)
                                                                                                                         gCTCCCG‑CAGCGATTCGC‑G‑CAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCa      >  4:150794/1‑135 (MQ=255)
                                                                                                                              cG‑CAGCGATTCGC‑G‑CAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTCCTCTCACATGGGGAACCCCCACACTACCATCgg  >  7:254539/1‑134 (MQ=255)
                                                                                                                              cG‑CAGCGATTCGC‑G‑CAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCgg  >  8:305414/1‑134 (MQ=255)
                                                                                                                                  |                                                                                                                                    
GGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGC‑GGCCGTGCCC‑GCTCCCG‑CAGCGATTCGC‑G‑CAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGG  >  ACB122/5301198‑5301455

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