Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,047,512 A→C 41.2% intergenic (‑104/+38) CJ019_00985 ← / ← CJ019_00986 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,047,5120AC41.2% 6.5 / 2.6 28intergenic (‑104/+38)CJ019_00985/CJ019_00986hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (14/2);  new base C (9/3);  total (23/5)
Fisher's exact test for biased strand distribution p-value = 6.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCGGCACAA‑GGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCA  >  ACB122/1047385‑1047638
                                                                                                                                 |                                                                                                                               
ccAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑CCCCCGCCCAA‑GGcc                                                                                                                         >  3:225249/1‑135 (MQ=255)
     cATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCGGCACAA‑GGCCGGGcc                                                                                                                    >  7:265462/1‑134 (MQ=255)
                           gCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑CCCCGGCACAA‑GGGCGGGCGTTTTCCGTGAAGGGCTTccgca                                                                                              >  4:73901/1‑135 (MQ=255)
                             ctcctGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCGGCACAACGGCCGGGC‑TTTGCGGGGAAGCGCAAcagca                                                                                              >  7:303219/1‑133 (MQ=255)
                             ctcctGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCGGCACAA‑GGTCGGGCGTTTGCGGTGAAGCGCTACAGCAgg                                                                                            >  8:288510/1‑135 (MQ=255)
                             ctcctGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCCGCCCAA‑GGCCGGGCGGTTGCCGTGAAGCGCTTCACCAgg                                                                                            >  8:380700/1‑135 (MQ=255)
                                  gTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCCGCACAA‑GGCCGGGCGTTTTGGGGGGAGCCCGACAGCACGATCa                                                                                        >  5:196244/1‑134 (MQ=255)
                                  gTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCCGCACAA‑GGCCGGGCGGTTGCGGTGAAGCGCTACAGCAGGAtaa                                                                                        >  6:149152/1‑132 (MQ=255)
                                     aGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCGGCACAA‑CGCCGGGCGTTTGCGGTGAAGCGCTACAGCCCGATCAGAAc                                                                                    >  3:349887/1‑135 (MQ=255)
                                      ggTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCCGCACAA‑GGCCGGGCGTTTGCGGTGCAGCGCGACACCACGATCAGAAca                                                                                   >  1:107592/1‑134 (MQ=255)
                                       gTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCCGCCCAA‑GGCCCGGCCTTTTCGGTGAAGCGCTTCAGCAGGAAAAGAAATc                                                                                  >  8:87772/1‑135 (MQ=255)
                                           gATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCCGCACAA‑GGCCGGGCGTTTGCGGTGAAGCGCTACAGCCGGCTCAGAACTAGTGa                                                                              >  6:190308/1‑135 (MQ=255)
                                                    ccTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCCGCCCAA‑GGCCGGGCGGTTTGGGGGGAGCCCGTCAGCACGCTCAGAACTCGTGTTACGcccgg                                                                     >  7:335158/1‑131 (MQ=255)
                                                          ggATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCGGCCCCA‑GGCCGGGCGGTTGCGGGGTAGCGCTACCGCAGGCTCAGAACTCGTGATCCGCAGTGTCCggg                                                               >  3:306666/1‑135 (MQ=255)
                                                                gCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCGGCACAA‑GGGCGGGCGTTTTGGGTGCAGCGCGTCAGCAGGGTCAGCACCCGGGATCCCCAGTGTGCGGGttgacg                                                         >  4:271284/1‑131 (MQ=255)
                                                                gCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCCGCCCAA‑GGCCGGGCGTTTGCGGGGGATCGCGACAGCAGGATCCGAACTCGGGCTCCGCCGGGTGCGGGGTGAgg                                                         >  8:57101/1‑135 (MQ=255)
                                                                gCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCCGCACAA‑GGGCCGGCGTTTTGGGTGGAGCGCTACAGCAGGGTCAGAACTCGGGGTCCGCAGTGTGCGGGTTCAgg                                                         >  1:234998/1‑135 (MQ=255)
                                                                  aTCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCCGCCCAA‑GGCCCGGCGTTTGCGGTGAAGCGCTACACCAGGATAAGAACTCGTGATCCCCAGTGTCCGGGGTCAGGtt                                                       >  3:254461/1‑135 (MQ=255)
                                                                        gTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCGGCCCAA‑GGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGGGATGCGCAGTGTCCGGGTTCAGGTTGGCAAt                                                 >  8:222516/1‑135 (MQ=255)
                                                                        gTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCGGCACAA‑GGCCGGGCGTTTGGGGTGTAGCGCTACAGGAGGGTCAGAACTCGGGCTGCGCCGTGTCCGGGGTCAGGTTGGCAAt                                                 >  8:216989/1‑135 (MQ=255)
                                                                              gCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCGGCCCAA‑GGCCGGGCGGTTGGGGGGAAGCGCTACAGCAGGATGATCACTCGTGATCCGCAGTGTGCGGGGTCAGGTTGGCAATGccccg                                           >  4:270126/1‑133 (MQ=255)
                                                                                            aGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCGGCACAA‑GGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGGTCAGGTTGGGCATGCCCAGCTTTCCCGCAGccg                             >  5:169092/1‑134 (MQ=255)
                                                                                                     cccgcccgC‑AAACGAAAAC‑GCCCGGCACAA‑GGCGGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTtcgatc                     <  4:304812/134‑1 (MQ=255)
                                                                                                      ccGGCATC‑AAACCCACACGGCCCGGCCCCT‑GGCGGGGGGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATGCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTtcgatcg                    <  3:284917/135‑1 (MQ=255)
                                                                                                      ccGGCACCAAAACCACAAC‑GCCCGGCCCAC‑GGCGGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTtcgatcg                    <  7:311735/135‑1 (MQ=255)
                                                                                                        ggCAGC‑AAACGAAAAC‑GCCCGGCACAA‑GGTCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGAcc                 <  7:288510/135‑1 (MQ=255)
                                                                                                             c‑aAACGAAAAC‑GCCCGGCACAA‑GGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCGGGTCGATCGACCCGGTc            >  7:302956/1‑135 (MQ=255)
                                                                                                                         cccccGCCCCAA‑GGCCGGGCGGTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGGTCGATCGACCCGGTCTGCTTGACCa  <  7:380700/133‑1 (MQ=255)
                                                                                                                                 |                                                                                                                               
CCAGACATCAAGAATCTGCAGCACTGTGCCTCCTGTTAGGTTTGATCGAAGGCCTACAGGATGAGCATCTGCGTGCACGCCAGGCAGCGGGAAGCGTTTTACCCGGCAGC‑AAACGAAAAC‑GCCCGGCACAA‑GGCCGGGCGTTTGCGGTGAAGCGCTACAGCAGGATCAGAACTCGTGATCCGCAGTGTCCGGGTTCAGGTTGGCAATGCCCAGCTTGCCCGCAGCCTGTTCGATCGACCCGGTCTGCTTGACCA  >  ACB122/1047385‑1047638

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