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

GGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGG  >  ACB122/5301198‑5301455
                                                                                                                                  |                                                                                                                                    
ggAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGGGCGGCCCGGGCCCCGCCCGCcgcggata                                                                                                                               >  7:212152/1‑131 (MQ=255)
ggAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCCGGGGACGCCGGGCCCCCCCCCGCCGCGaata                                                                                                                               >  1:208820/1‑132 (MQ=255)
     tGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCa                                                                                                                        >  8:304522/1‑135 (MQ=255)
                      cGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGGGCGGGCCGGGCCCCTCCCCCCGCGCTTGGCGCCGGGGGGGGTTTTTTGcc                                                                                                        >  8:260091/1‑135 (MQ=255)
                         aaGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGGGCGGCCCGGCCCGCTCCCCCCGCCCGTCGCGCAGGTCGAGGTTTTTGGCcac                                                                                                     >  7:240771/1‑134 (MQ=255)
                                            aTGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGGGCCGCCGGGGCCCCGCCCCCCGCGCTTGGCGGCGGGGCCGTGTTTTTGCCCCAAAAACAAAACCCCCAccag                                                                                 >  4:311300/1‑133 (MQ=255)
                                                    gagGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGGGCCGGCGGGGCCCCGCCCCCCGCGGATAGCGGCGGGCGGGGTTTTTTGCCCCCAACACAAAACACCCACCAGGGGgcgcg                                                                         >  7:20955/1‑134 (MQ=255)
                                                    gagGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGCGCGGGCGGGGCCGCGCCCGCCGCGCTTGCCGCCCGGGGCGGTTTTTGGCCACAAAAACAAAAACCCAACCagcgtgc                                                                            >  2:258673/1‑132 (MQ=255)
                                                                       gTACATGAAACGGCCGAATGCCATCAGGCTGCGGGGGCGGCCGGGGCCCCGCCCCCCCCGCTTGGCGCAGGGGGAGGTTTTTTGCCACAAAGACAAAAACCCCACCAGCGGGCGCCCATAGGGGGTTTGGGCAAt                                                      >  2:237849/1‑135 (MQ=255)
                                                                         aCATGAAACGGCCGAATGCCATCAGGCTGCGGGCGCGGCCGTGCCCGCTCCCGCAGCGCTTCGCGCAGGTCGACTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCCCAGAAAGGGGTTTTGGGAAAtt                                                     >  8:317486/1‑133 (MQ=25)
                                                                           aTGAAACGGCCGAATGCCATCAGGCTGCGGGGGCGGGCCTTCCCCCGCCCCCCGCGCATACCGCAGGGCGAGGTTTTTTGCCACAAAAACAAAAACCCAAACAACGGGCGCCCACAGGGGTTTTTGGaaaataa                                                   >  6:199886/1‑130 (MQ=255)
                                                                                                      gggtggcgggcggGCCCGCGCCCGACCGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATgg                        <  8:330044/130‑1 (MQ=14)
                                                                                                          gccgggCCGGGCCCGGCCCCCGCCGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGa                     <  1:207496/131‑1 (MQ=14)
                                                                                                              ggCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAAccc                <  8:129861/135‑1 (MQ=255)
                                                                                                                  gTGCCCGCTCCCGCAGCGATTAGCGCAGGTGGAGTTTTTTTCCCACAAAGACAAAAACCCAACCAGCGTGCGGAGGTTGGGGGTTTGCGAAATGAATGTTGTAGATGACCTACTCTCACAAGGGGAAAccccccc            >  2:295903/1‑132 (MQ=255)
                                                                                                                    gcccgcccCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATgggg                      >  6:375220/1‑123 (MQ=255)
                                                                                                                    gCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAAAGAATCTTGACGATGCCCTACACGCACATGGGGCAACCCCCCACt          >  7:248739/1‑135 (MQ=17)
                                                                                                                    gCCCCCCCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATgggg                      <  5:375220/123‑1 (MQ=255)
                                                                                                                         gCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCa      >  4:150794/1‑135 (MQ=255)
                                                                                                                              cGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTCCTCTCACATGGGGAACCCCCACACTACCATCgg  >  7:254539/1‑134 (MQ=255)
                                                                                                                              cGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCgg  >  8:305414/1‑134 (MQ=255)
                                                                                                                                  |                                                                                                                                    
GGAGATGGGCGGTGTTGAGCATCGGAAGGTCTCTGGAGGCGCTGATGGCTGAGAGGGGAACGATATCACGCGTACATGAAACGGCCGAATGCCATCAGGCTGCGGGAGCGGCCGTGCCCGCTCCCGCAGCGATTCGCGCAGGTCGAGTTTTTTGCCCACAAAGACAAAACCCCTACCTGCGTGCGCAGATAGGGGTTTTGCGAAATGAATCTTGACGATGACCTACTCTCACATGGGGAAACCCCACACTACCATCGG  >  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: