Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,895,319 C→T 9.5% R118C (CGC→TGC)  CJ019_03685 → putative protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,895,3190CT 9.5% 119.9 / 4.0 42R118C (CGC→TGC) CJ019_03685putative protein
Reads supporting (aligned to +/- strand):  ref base C (25/13);  new base T (2/2);  total (27/15)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-01

CAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGAGGACTGCGACAAGGTCTACGCCCAGGCCCTGGCCGCCGGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGTCCACCCACAAAGAGGACCTGACCCCCGACGAGATTCGCGCCCGTGCCGCGAAGATGTTCGGCGGTAGCCAAGCGA  >  ACB122/3895193‑3895439
                                                                                                                              |                                                                                                                        
cagAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGAGGACTGCGACAAGGTCTACGCCCAGGCCCTGGCCGCCGGGGGCATCCAGTTACATCCCTTGCCCGACCAGTTTTACGGTGACCGCAGCGGc                                                                                                                  >  1:103051/1‑135 (MQ=255)
           cGGTGCGGGTGTTGGCCTGCACTTGTATGTTGAGGACTGCGACAAGGTCTACGCCCAGGCCCTGGCCGCCGGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGCCCGCAGCGGCCCCCTg                                                                                                            >  2:148759/1‑130 (MQ=255)
                       tGGCCTGCACTTGTATGTTGAGGACTGCGACAAGGTCTACGCCCAGGCCCTGGCCGCCGGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCaa                                                                                           <  5:229225/135‑1 (MQ=255)
                                 ttGTATGTTGAGGACTGCGACAAGGTCTACGCCCAGGCCCTGGCCGCCGGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCCCCCTGAAAGATCCGTTTGGCAACATCTGGTTc                                                                                 >  5:144590/1‑135 (MQ=255)
                                 ttGTATGTTGAGGACTGCGACAAGGTCTACGCCCAGGCCCTGGCCGCCGGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTc                                                                                 >  8:86093/1‑135 (MQ=255)
                                       gTTGAGGACTGCGACAAGGTCTACGCCCAGGCCCTGGCCGCCGGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGTcc                                                                           >  1:158563/1‑135 (MQ=255)
                                         tGAGTACTGCGACAAGGTCTACGCCCAGGCCCTGGCCGCCGGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGTcca                                                                          <  5:165809/134‑1 (MQ=255)
                                            ggACTGCGACAAGGTCTACGCCCAGGCCCTGGCCGCCGGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGTccaccc                                                                       >  8:87025/1‑134 (MQ=255)
                                             gACTGCGACAAGGGCTACGCCCAGGCCCTGGCCGCCGGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGTccacccac                                                                     <  8:207969/135‑1 (MQ=255)
                                              aCTGCGACAAGGTCTACGCCCAGGCCCTGGCCGCCGGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGttt                                                                                                >  8:20076/1‑107 (MQ=255)
                                              aCTGCGACAAGGTCTACGCCCAGGCCCTGGCCGCCGGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGttt                                                                                                <  7:20076/107‑1 (MQ=255)
                                                   gACAAGGTCTACGCCCAGGCCCTGGCCGCCGGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCgt                                                                               <  3:297527/119‑1 (MQ=255)
                                                   gACAAGGTCTACGCCCAGGCCCTGGCCGCCGGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCgt                                                                               >  4:297527/1‑119 (MQ=255)
                                                           cTACGCCCAGGCCCTGGCCGCCGGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGCCCGCAGCGGCACCCTGAAAGATCCGATTGGCAACATCTGGTTCGTCTCAACCCACAAAGAGGACCTGAc                                                       >  3:268866/1‑135 (MQ=255)
                                                               gCCCAGGCCCTGGCCGCCGGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGTCCACCCACAAAGAGGACCTGAccccc                                                   >  1:55518/1‑135 (MQ=255)
                                                               gCCCAGGCCCTGGCCGCCGGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGCTGCCCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGACCACCCACAAAGAGGACCTGAcgcc                                                    >  7:60586/1‑134 (MQ=255)
                                                                cccAGGCCCTGGCCGCCGGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGTCCACCCACAAAGAGGACCTGACCCCcg                                                  <  4:161311/135‑1 (MQ=255)
                                                                 ccAGGCCCTGGCCGCCGGGGGCATCCAGTTACATCCCTTGCCCGACCAGTTTTACGGTGCCCGCAGCGGCACCCTGAAAGATCCGTTTGGCACCATCTGGTTCGTGTCCACCCACAAAGAGGACCTGACCCCcga                                                 >  8:281020/1‑135 (MQ=255)
                                                                      cccTGGCCGCCGGGGGCATCCAGTTACATCCCTTGCCCGCCCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGAGCCGTTTGGCGACATCTGGTTCGTGTCCACCCACCAAGAGGACCTGCCCCCCGAAgaga                                            >  5:167846/1‑135 (MQ=255)
                                                                       ccTGGCCGCCGGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGTCCACCCACAAAGAGGACCTGACCCCCGACGAGAt                                           >  8:202147/1‑135 (MQ=255)
                                                                       ccTGGCCGCCGGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGTCCACCCACAAAGAGGACCTGACCCCCGACGAGAt                                           >  2:21989/1‑135 (MQ=255)
                                                                       ccTGGCCGCCGGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGTCCACACACAAAGAGGACCTGAACCCCGACGAGAt                                           >  3:201550/1‑135 (MQ=255)
                                                                             cgccgGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAATATCCGTTTGGCAACATCTGGTTCGGGTCCACCCACAAAGAGGACCTGACCCCCGACGAGATTCgcgc                                     <  5:151835/135‑1 (MQ=255)
                                                                             cgccgGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGGCCACCCACAAAGAGGACCTGACCCCCGACGAGATTCgcgc                                     <  8:252534/135‑1 (MQ=255)
                                                                              gccgGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGTCCACCCACCAAGAGGACCTGCCCCCCGACGAGATTCGCGcc                                    >  7:248895/1‑135 (MQ=255)
                                                                              gccgGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGTCCACCCACAAAGAGGCCCTGACCCCCGACGAGATTCGCGcc                                    >  4:165883/1‑135 (MQ=255)
                                                                                   gggCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGTCCACCCACAAAGAGGACCTGACCCCCGACGAGATTCGCGCCCGTGc                               >  3:119222/1‑135 (MQ=255)
                                                                                       aTCCAGTTACATCCCCTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGTCCACCCACAAAGAGGACCTGACCCCCGACGAGATTCGCGCCCGTGCcgcg                           <  8:208363/135‑1 (MQ=255)
                                                                                            gTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGGCCACCCACAAAGAGGACCTGACCCCCGACGAGATTCGCGCCCGTGCCGCGAAGAt                      <  7:86093/135‑1 (MQ=255)
                                                                                                          ccgaccAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGTCCACCCACAAAGAGGACCTGACCCCCGACGAGATTCGCGCCCGTGCCGCGAAGATGTTCGGCGGTAGc         >  5:242423/1‑134 (MQ=255)
                                                                                                          ccgaccAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGTCCACCCACAAAGAGGACCTGACCCCCGACGAGATTCGCGCCCGTGCCGCGAAGATGTTCGGCGGTAGc         >  8:49037/1‑134 (MQ=255)
                                                                                                              ccAGTTTTCCGGTGACCGCAGCGGCACCCTGAACGAGCCGTTTGGCAACATCTGGTGCGTGCCCACCCACAACGAGGACCTGACCCCCGCCGAGATGCGCGCCCGTGCCGCGAAGATGTTCGGCGGTAGCCAAGc    <  3:165883/135‑1 (MQ=255)
                                                                                                              ccAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGTCCACCCACAAAGAGGACCTGACCCCCGACGAGATTCGCGCCCGTGCCGCGAAGATGTTCGGCGGTAGCCAAGc    <  7:202147/135‑1 (MQ=255)
                                                                                                              ccAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCCTTTGGCAACAACTGGTTCGTGTCCACCCCCAAAGACGCCCTGACCCCCCACGAGATTCGCGCCCGTGCCGCGAAGATCTTCGGCGGTAGCCAAGc    >  2:187774/1‑135 (MQ=255)
                                                                                                                aGTTTTACGGTGACTGCAGCGGCACCCTGAAAGATc                                                                                                     <  2:132545/36‑1 (MQ=255)
                                                                                                                aGTTTTACGGTGACTGCAGCGGCACCCTGAAAGATc                                                                                                     <  2:128680/36‑1 (MQ=255)
                                                                                                                aGTTTTACGGTGACTGCAGCGGCACCCTGAAAGATc                                                                                                     >  1:132545/1‑36 (MQ=255)
                                                                                                                aGTTTTACGGTGACTGCAGCGGCACCCTGAAAGATc                                                                                                     >  1:128680/1‑36 (MQ=255)
                                                                                                                aGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCgt                                                                               <  3:73009/58‑1 (MQ=255)
                                                                                                                aGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCgt                                                                               >  4:73009/1‑58 (MQ=255)
                                                                                                                aGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGTCCACCCACAAAGAGGACCTGACCCCCGACGAGATTCGCGCCCGTGCCGCGAAGATGTTCGGCGGTAGCCAAGcga  >  7:161401/1‑135 (MQ=255)
                                                                                                                aGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGTCCACCCACAAAGAGGACCTGACCCCCGACGAGATTCGCGCCCGTGCCGCGAAGATGTTCGGCGGTAGCCAAGcga  >  8:293203/1‑135 (MQ=255)
                                                                                                                              |                                                                                                                        
CAGAAGCTCAACGGTGCGGGTGTTGGCCTGCACTTGTATGTTGAGGACTGCGACAAGGTCTACGCCCAGGCCCTGGCCGCCGGGGGCATCCAGTTACATCCCTTGACCGACCAGTTTTACGGTGACCGCAGCGGCACCCTGAAAGATCCGTTTGGCAACATCTGGTTCGTGTCCACCCACAAAGAGGACCTGACCCCCGACGAGATTCGCGCCCGTGCCGCGAAGATGTTCGGCGGTAGCCAAGCGA  >  ACB122/3895193‑3895439

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: