Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 4,124,901 C→T 7.9% G234R (GGG→AGG)  CJ019_03883 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1224,124,9010CT 7.9% 145.9 / 3.7 52G234R (GGG→AGG) CJ019_03883hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (24/22);  new base T (2/2);  total (28/24)
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.86e-01

GAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCAAGCCCCGCACTCACCCCAAGCATTTCGCACGGCACCAGCACAATCCACAGAATA  >  ACB122/4124771‑4125034
                                                                                                                                  |                                                                                                                                     
gAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCCTCCAGCACAAACCCCAGCa                                                                                                                                  >  3:240648/1‑136 (MQ=255)
         aaGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCAACCAGCGCCGATGCAGCCCACGCGCCAGAGCAGCCAGCCCAACCCCCAGCACGCCAAt                                                                                                                           >  8:222761/1‑134 (MQ=255)
           ggCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCCGCACGCCAATcag                                                                                                                        >  8:247797/1‑135 (MQ=255)
           ggCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATcag                                                                                                                        >  4:202407/1‑135 (MQ=255)
           ggCCCAAGGCCGGT‑ATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGTCCTTTCAGCCATGAACCCAGCGCCGATGCAGCCCACCCGCCAGAGCATCCAGCCCACACCCCCGCACGCCAATcagc                                                                                                                       >  6:327029/1‑135 (MQ=38)
             cccAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATcagca                                                                                                                      >  1:224153/1‑135 (MQ=255)
                     cGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATc                                                                                                              >  1:361262/1‑135 (MQ=255)
                       gTGATGTACCTGGGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGc                                                                                                            <  2:224153/135‑1 (MQ=255)
                          aTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCa                                                                                                          >  1:189026/1‑134 (MQ=255)
                          aTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCAt                                                                                                         >  4:374541/1‑135 (MQ=255)
                                cTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCCTCGCCATCAGTTc                                                                                                   >  5:18195/1‑135 (MQ=255)
                                     ggAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCCGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCCTCAGTTCGGAAt                                                                                              >  5:229435/1‑135 (MQ=255)
                                     ggAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGCTGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATa                                                                                             >  8:297621/1‑136 (MQ=255)
                                           gccgcccccGCGAAGAAGCCCCATTGGCCTTTCGGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAg                                                                                       <  8:307106/130‑1 (MQ=255)
                                            gCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAg                                                                                       >  1:264292/1‑135 (MQ=255)
                                             cGAGCCCGCGAAGAAGCCCCATTGGCCTTTCCGCCATGCACCCAGCGCCGCTGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCCTCAGTTCGGAAGACGCCCGc                                                                                      >  2:276480/1‑135 (MQ=255)
                                               cgCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCgg                                                                                   <  8:256105/135‑1 (MQ=255)
                                                cccggcgAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCgg                                                                                   <  1:33988/131‑1 (MQ=255)
                                                cccccCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCgg                                                                                   <  5:344025/134‑1 (MQ=255)
                                                        agaagCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGt                                                                           >  1:11256/1‑135 (MQ=255)
                                                          aagCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTAt                                                                         <  2:124279/135‑1 (MQ=255)
                                                                  ttGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATg                                                                 <  6:375688/135‑1 (MQ=255)
                                                                  ttGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATg                                                                 <  2:189026/135‑1 (MQ=255)
                                                                            aGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCaa                                                       >  1:159410/1‑135 (MQ=255)
                                                                             gCCATGCACCCATCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACGCAAg                                                      <  2:361262/135‑1 (MQ=255)
                                                                             gCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCTCAGCACGCCAATcagca                                                                                                                      >  5:227960/1‑71 (MQ=255)
                                                                             gCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCTCAGCACGCCAATcagca                                                                                                                      <  6:227960/71‑1 (MQ=255)
                                                                             gCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCTCAGCACGCCAATcagca                                                                                                                      <  6:223217/71‑1 (MQ=255)
                                                                             gCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCTCAGAACGCCAATcagca                                                                                                                      >  5:223217/1‑71 (MQ=255)
                                                                             gCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAgcagca                                                                                                                   <  4:24736/74‑1 (MQ=255)
                                                                             gCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAgcagca                                                                                                                   >  3:24736/1‑74 (MQ=255)
                                                                             gCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCAAg                                                      >  5:146456/1‑135 (MQ=255)
                                                                             gCCATGCACCCAGCGCCGAGGCAGCCCTCGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCAAg                                                      <  1:323794/135‑1 (MQ=255)
                                                                                    aCCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCACAGCACGCCCATCAGCATCACCAACGCACGAAGGTACGAAAAGGCCAGACGGACACGCGTATAAAGAAAGAAATACGCCATCCCCGCa                                               >  6:196424/1‑135 (MQ=255)
                                                                                           gcCGATGCCGCCCACGCGCCACAGCATCCAGCTCAAACCACAGCACGCCACTCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCTCGGCTATCAAGAATGACATCCCCACGCCCCGCACTCcccc                                        >  8:141025/1‑135 (MQ=255)
                                                                                                tGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATTGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCAAGCCCCGCACTCACCCCAAGCa                                  <  3:202407/136‑1 (MQ=255)
                                                                                                     cccACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCAAGCCCCGCACTCACCCCAAGCATTTCg                             <  4:240648/136‑1 (MQ=255)
                                                                                                     cccACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCAAGCCCCGCACTCACCCCAAGCATTTCg                             <  7:297621/136‑1 (MQ=255)
                                                                                                     cccACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCAAGCCCCGCACTCACCCCAAGCATTTCg                             <  7:247797/136‑1 (MQ=255)
                                                                                                     cccACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTGCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCAAGCCCCGCACTCACCCCAAGCATTTc                              <  5:117163/135‑1 (MQ=255)
                                                                                                            gcCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCAAGCCCCGCACTCACCCCAAGCATTTCGCACGGc                       >  2:236373/1‑135 (MQ=255)
                                                                                                               agagCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCAAGCCCCGCACTCACCCCAAGCATTTCGCACGGCACCa                   <  3:374541/136‑1 (MQ=255)
                                                                                                                gagCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCAAGCCCCGCACTCACCCCAAGCATTTCGCACGGCACCa                   <  6:396155/135‑1 (MQ=255)
                                                                                                                      cctgctcaAACCCCAGCTCGCCAATCAGTAGCACCATGGCCATCTGTTCGTATTACGCCAGCCGGGCACGGGTAGCAAGAATGAAATACCCAAGCCCCGCACTCACCGCAAGCATTTCGAACGGCACCAGCACAAt            <  7:141025/130‑1 (MQ=255)
                                                                                                                         gCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCAAGCCCCGCACTCACCCCAAGCATTTCGCACGGCACCAGcac               <  3:197497/130‑1 (MQ=255)
                                                                                                                         gCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCAAGCCCCGCACTCACCCCAAGCATTTCGCACGGCACCAGcac               >  5:294392/1‑130 (MQ=255)
                                                                                                                         gCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCAAGCCCCGCACTCACCCCAAGCATTTCGCACGGCACCAGcac               >  4:197497/1‑130 (MQ=255)
                                                                                                                         gCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCAAGCCCCGCACTCACCCCAAGCATTTCGCACGGCACCAGcac               <  6:294392/130‑1 (MQ=255)
                                                                                                                         gCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCAAGCCCCGCACTCACCCCAAGCATTTCGCACGGCACCAGCACAAt            <  6:44573/133‑1 (MQ=255)
                                                                                                                         gCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCAAGCCCCGCACTCACCCCAAGCATTTCGCACGGCACCAGCACAAt            >  5:44573/1‑133 (MQ=255)
                                                                                                                         gCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCAAGCCCCGCACTCACCCCAAGCATTTCGCACGGCACCAGCACAATcc          >  6:334142/1‑135 (MQ=255)
                                                                                                                         gCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCAAGCCCCGCACTCACCCCAAGCATTTCGCACGGCACCAGCACAATCca         >  7:132862/1‑136 (MQ=255)
                                                                                                                                ccccAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCAAGCCCCGCACTCACCCCAAGCATTTCGCACGGCACCAGCACAATCCACAGAATa  <  8:29344/136‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                     
GAATACCGTAAGGCCCAAGGCCGGTGATGTACCTGTGGGAGCGGGCGAGCCCGCGAAGAAGCCCCATTGGCCTTTCAGCCATGCACCCAGCGCCGATGCAGCCCACGCGCCAGAGCATCCAGCACAAACCCCAGCACGCCAATCAGCAGCACCATCGCCATCAGTTCGGAATACGCCAGCCGGTCACGGGTATCAAGAATGAAATACCCAAGCCCCGCACTCACCCCAAGCATTTCGCACGGCACCAGCACAATCCACAGAATA  >  ACB122/4124771‑4125034

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