Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,128,518 C→T 10.2% T130T (ACG→ACA)  CJ019_04823 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,128,5180CT10.2% 113.2 / 3.4 39T130T (ACG→ACA) CJ019_04823hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (16/19);  new base T (2/2);  total (18/21)
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

CGCCAGCAGTGCATCGAAGCAGGCCGGTTGCTCCAGCGTGTTCAATGCCCCGACCACCGCCTCGATCGTCGACAACGCACCCGCCTCCGGCGCTTTGCGCAACCGATAGCGTGAAGCGGCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTGCGCGCCTTGCGCCAGGTGCCGTCGGGCACGATCAGCAGCAAAGGTTTGTCGCCCGCCTCGCCATAAGCTGTCAGCACCTGGGCGT  >  ACB122/5128386‑5128651
                                                                                                                                    |                                                                                                                                     
cgcCAGCAGTGCATCGAAGCAGGCCGGTTGCTCCAGCGTGTTCAATGCCCCGACCACCGCCTCGATCGTCGACAACGCACCCGCCTCCGGCGCTTTGCGCAACCGATAGCGTGAAGCGGCCACCTGCCCCAGCGt                                                                                                                                     <  7:138735/135‑1 (MQ=255)
   cagcagTGCATCGAAGCAGGCCGGTTGCTCCAGCGTGTTCAATGCCCCGCCCACCGCCTCGATCGTCGACAACGCACCCGCCTCCGGCGCTTTGCGCAACCGATAGCGTGAAGCGGCCACCTGCCCCAGCGTCAc                                                                                                                                  >  3:210126/1‑135 (MQ=255)
          gCATCGAAGCAGGCCGGTGGCTCCCGCGTGTTCAATGCCCCGACCACCGCCTCGATCGTCGACAACGCACCCGCCTCCGGCGCTTTGCGCAACCGATAGCGTGAAGCGGCCACCTGCCCCAGCGTCACCCGTGGc                                                                                                                           <  5:191627/135‑1 (MQ=255)
           cATCGAAGCAGGCCGGTTGCTCCAGCGTGTTCAATGCCCCGCCCACCGCCTCGATCGTCGACAACGCACCCGCCTCCGGCGCTTTGCGCAACCGATAGCGTGAAGCGGCCACCTGCCCCAGCGTCACCCGTGGCa                                                                                                                          >  8:172892/1‑135 (MQ=255)
            aTCGAAGCAGGCCGGTTGCTCCAGCGTGTTCAATGCCCCGACCACCGCCTCGATCGTCGACAACGCACCCGCCTCCGGCGCTTTGCGCAACCGATAGCGTGAAGCGGCCACCTGCCCCAGCGTCACCCGTGGCAg                                                                                                                         >  1:158247/1‑135 (MQ=255)
             tCGAAGCAGGCCGGTTGCTCCAGCGTGTTCAATGCCCCGACCACCGCCTCGATCGTCGACAACGCACCCGCCTCCGGCGCTTTGCGCAACCGATAGCGTGAAGCGGCCACCTGCCCCAGCGTCACCCGTGGCAGc                                                                                                                        >  3:258731/1‑135 (MQ=255)
                 aGCAGGCCGGTTGCTCCAGCGTGTTCAATGCCCCGACCACCGCCTCGATCGTCGACAACGCCCCCGCCTCCGGCGCTTTGCGCAACCGATAGCGTGAAGCGGCCACCAGCCCCAGCGTCACCCGTGGCAGccccc                                                                                                                    >  4:42893/1‑135 (MQ=255)
                     ggCCGGTTGCTCCAGCGTGTTCAATGCCCCGACCACCGCCTCGATCGTCGACAACGCACCCGCCTCCGGCGCTTTGCGCAACCGATAGCGTGAAGCGGCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAg                                                                                                                <  5:112358/135‑1 (MQ=255)
                      gCCGGTTGCTCCAGCGTGTTCAATGCCCCGACCACCGCCTCGATCGTCGACAACGCACCCGCCTCCGGCGCTTTGCGCAACCGATAGCGTGAAGCGGCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAg                                                                                                                >  8:210554/1‑134 (MQ=255)
                                           aaTGCCCCGACCACCGCCTCGATCGTCGACAACGCACCCGCCTCCGGCGCTTTGCGCAACCGATAGCGTGAAGCGGCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACagcagc                                                                                          >  4:212328/1‑135 (MQ=255)
                                                cccGACCACCGCCTCGATCGTCGACAACGCACCCGCCTCCGGCGCTTTGCGCAACCGATAGCGTGAAGCGGCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTgcg                                                                                     >  7:105907/1‑135 (MQ=255)
                                                   gACCACCGCCTCGATCGTCGACAACGCACCCGCCTCCGGCGCTTTGCGCAACCGATAGCGTGAAGCGGCCACCTGCCCCAGCGGCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTgcgcgc                                                                                  <  5:207063/135‑1 (MQ=255)
                                                       accGCCTCGATCGTCGACAACGCACCCGCCTCCGGCGCTTTGCGCAACCGATAGCGTGAAGCGGCCACCTGCCCCCGCGTCACCCGTGGCAGCCCCCCCCGCAACGGGTTCAGGTACAGCAGCTTGCGCGCCTTg                                                                              >  4:165117/1‑135 (MQ=255)
                                                          gCCTCGATCGTCGACAACGCACCCGCCTCCGGCGCTTTGCGCAACCGATAGCGTGAAGCGGCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTGCGCGCCTTgcgc                                                                           <  8:134152/135‑1 (MQ=255)
                                                            ctcgctcgTCGACACCGCACCCGCCTCCGGCGCTTTGCGCACCCGAGAGCGTGACGCGGCCACCTGCCCCCGTGTCCCCCGTGGCAGCCCCCCCCGCCCCGGGTTCAGGTCCCGCAGCTGGCGCGCCTTGCGCCa                                                                         >  2:231961/1‑135 (MQ=255)
                                                            cTCGATCGTCGACAACGCACCCGCCTCCGGCGCTTTGCGCAACCGATAGCGTGAAGCGGCCACCTGCCCCAGTGTCACCCGTGGCAGCCCCCCCAGCACCGGGTTCAGGTACAGCAGCTTGCGCGCCTTGCGCCa                                                                         >  2:231962/1‑135 (MQ=255)
                                                                                   ccTCCGGCGCTTTGCGCAACCGATAGCGTGAAGCGGCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTAGAGCAGCTTGCGCGCCTTGCGCCAGGTGCCGTCGGGCACGATcagca                                                  <  7:196566/135‑1 (MQ=255)
                                                                                   ccTCCGGCGCTTTGCGCAACCGATAGCGTGAAGCGGCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTGCGCGCCTTGCGCCAGGTGCCGTCGGGCACGATcagca                                                  <  4:210126/135‑1 (MQ=255)
                                                                                                 cgcAACCGATAGCGTGAAGCGGCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTGCGCGCCTTGCGCCAGGTGCCGTCGGGCACGATCAGCAGCAAAGGTTTGTcg                                    >  6:193504/1‑135 (MQ=255)
                                                                                                   cAACCGATAGCGTGAAGCGGCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTGCGCGCCTTGCGCCAGGTGCCGTCGGGCACGATCAGCAGCAAAGGTTTGTcgcc                                  <  5:193504/135‑1 (MQ=255)
                                                                                                    aaCCGATAGCGTGAAGCGGCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTGCGCGCCTTGCGCCAGGTGCCGTCGGGCACGATCAGCAGCAAAGGTTTGTcgccc                                 <  3:212328/135‑1 (MQ=255)
                                                                                                    aaCCGATAGCGTGAAGCGGCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTGCGCGCCTTGCGCCAGGTGCCGTCGGGCACGATCAGCAGCAAAGGTTTGTcgccc                                 <  8:105907/135‑1 (MQ=255)
                                                                                                         aTAGCGTGAAGCGGCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTGCGCGCCTTgcgc                                                                           <  6:28203/88‑1 (MQ=255)
                                                                                                         aTAGCGTGAAGCGGCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTGCGCGCCTTgcgc                                                                           >  5:28203/1‑88 (MQ=255)
                                                                                                              gTGAAGCGGCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTGCGCGCCTTGCGCCAGGTGCCGTCGGGCACGATCAGCAGCAAAGGTTTGTCGCCCGCCTCGCCAt                       <  2:15301/135‑1 (MQ=255)
                                                                                                                  aGCGGCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAGCTACGGGTTCAGGTACAGCAGATTGCGCGCCTTGCGCCAGGTTCCGTCGGGCACGATCAGCAGCATAGGTTTGTCGCCCGCCTCTCCATAAGc                   <  4:234641/135‑1 (MQ=255)
                                                                                                                  aGCGGCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTGCGCGCCTTGCGCCAGGTGCCGTCGGGCACGATCAGCAGCAAAGGTTTGTCGCCCGCCTCGCCATAAGc                   <  2:208840/135‑1 (MQ=255)
                                                                                                                   gCGGCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTGCGCGCCTTGCGCCAGGTGCCGTCGGGCACGATCAGCAGCAAAGGTTTGTCGCCCGCCTCg                           >  7:302456/1‑126 (MQ=255)
                                                                                                                   gCGGCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTGCGCGCCTTGCGCCAGGTGCCGTCGGGCACGATCAGCAGCAAAGGTTTGTCGCCCGCCTCg                           <  8:302456/126‑1 (MQ=255)
                                                                                                                   gCGGCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTGCGCGCCTTGCGCCAGGTGCCGTCGGGCACGATCAGCAGCAAAGGTTTGTCGCCCGCCTCGCCATAAGCt                  >  4:95860/1‑135 (MQ=255)
                                                                                                                     ggCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTGCGCGCCTTGCGCCAGGTGCCGTCGGGCACGATCAGCAGCAAAGGTTTGTCGCCCGCCTCGCCATAAGCTg                 >  8:164394/1‑134 (MQ=255)
                                                                                                                         aCCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTAGAGCAGCTTGCGCGCCTTGCGCCAGGTGCCGTCGGGCACGATCAGCAGCAAAGGTTTGTCGCCCGCCTCGCCATAAGCTGTCAGc            <  8:135869/135‑1 (MQ=255)
                                                                                                                            tGCCCCAGTGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTGCGCGCCTTGCGCCAGGTGCCGTCGGGCACGATCAGCAGCAAAGGTTTGTCGCCCGCCTCGCCATAAGCTGTCAGCAcc         <  1:231962/135‑1 (MQ=255)
                                                                                                                            tGCCCCAGTGGCACCCGTGGCAGCCCCCCCAGCCACGGGTTCAGGTACAGCAGCTTGCGCGCCTTGCGCCAGGGGCCGTCGGGCACGATCAGCAGCAAAGGTTTGTCGCCCGCCTCGCCATAAGCTGTCAGCAcc         <  1:231961/135‑1 (MQ=255)
                                                                                                                              ccccAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTGCGCGCCTTGCGCCAGGTGCCGTCGGGCACGATCAGCAGCAAAGGTTTGTCGCCCGCCTCGCCATAAGCTGTCAGCACCTg       <  5:22616/135‑1 (MQ=255)
                                                                                                                                 cAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCCGGTTCAGCAGCTTGCGCGCCTTGCGCCAGGTGCCGTCGGGCACGATCAGCAGCAAAGGTTTGTCGCCCGCCTCGCCATAAGCTGTCAGCCCCTGGGc    >  3:267129/1‑135 (MQ=255)
                                                                                                                                 cAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTgcgc                                                                                    <  2:52668/55‑1 (MQ=255)
                                                                                                                                 cAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTgcgc                                                                                    >  1:52668/1‑55 (MQ=255)
                                                                                                                                   gCGTCACCCGGGGCAGCCCCCCCCGCAACGGGTTCAGGTACAGCAGCTTGCGCGCCTTGCGCCAGGTGCCGTCGGGCACGATCAGCAGCAAAGGTTTGTCGCCCGCCTCGCCATAAGCTGTCAGCACCTGGGCGt  <  4:295273/135‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                     
CGCCAGCAGTGCATCGAAGCAGGCCGGTTGCTCCAGCGTGTTCAATGCCCCGACCACCGCCTCGATCGTCGACAACGCACCCGCCTCCGGCGCTTTGCGCAACCGATAGCGTGAAGCGGCCACCTGCCCCAGCGTCACCCGTGGCAGCCCCCCCAGCAACGGGTTCAGGTACAGCAGCTTGCGCGCCTTGCGCCAGGTGCCGTCGGGCACGATCAGCAGCAAAGGTTTGTCGCCCGCCTCGCCATAAGCTGTCAGCACCTGGGCGT  >  ACB122/5128386‑5128651

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