Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,491,319 A→C 48.4% intergenic (‑681/+40) CJ019_03306 ← / ← CJ019_03307 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,491,3190AC48.4% ‑1.2 / 9.5 20intergenic (‑681/+40)CJ019_03306/CJ019_03307hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (5/4);  new base C (8/2);  total (13/7)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.49e-01

CGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATG  >  ACB122/3491212‑3491433
                                                                                                           |                                                                                                                  
gggtgggGCCGCCCGGTTGCAGGAGGGGGTTACCCCGCCAATAGGCCGCACCGGCCGGCAACTATCCCTGGCAATCCCCCGACACACCCCCCACCGCAGCCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCg                                                                                         <  8:214970/134‑1 (MQ=255)
 ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGGCGGCGGCCCCCGGACCCCGGCgccgcccc                                                                                        >  1:170049/1‑133 (MQ=18)
 ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACAACGGCGCAGCCgg                                                                                        >  6:114566/1‑135 (MQ=255)
 ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCCGGCCGCCGGCCCGCCCACACCCGCGCCGCCgg                                                                                        >  6:375645/1‑135 (MQ=21)
     gtgtCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGtt                                                                                      >  8:350478/1‑133 (MQ=255)
      gggcGCCAGGTTGTCGGCGGTGGTTTCCCCGCAAATCGGCCGCCCCGGCCGGCCACTAGCCCTGGCCATACCCCGCCCACCCCCCCACCGCCTGCGCGGGCCGCCGGCCCGCGAACACCGGCGCAGCCGGTTTCa                                                                                   <  8:145957/132‑1 (MQ=255)
              ggTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCAATATCCGTGGGCATTCACCGTCCCAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTc                                                                           <  8:286533/134‑1 (MQ=255)
                  gTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCGGCCCCCGCACACCGGCGCCCCCCGGTTGATCCCACACAAcc                                                                       >  3:331970/1‑134 (MQ=17)
                      gAGGTGGTTTCCCCGCTATTTGGCCGTACCGGCGGGCCAATACCCCTGGGAATACCCGGTCCCCCCCGCCCCCGCATGCGCGGGCCGCCGGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGAcc                                                                    <  2:205214/134‑1 (MQ=255)
                        ggtgttttACCCGCTAATTGGCGGGACCTCCCGCCAAATACCCCTGTGCATTCCCCGCCCCACCCCCCCCCGCATCCGCGGGCAGCCTGCCCGCCAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATc                                                                 <  6:281146/130‑1 (MQ=255)
                             tttACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCGGCCCCCGAACACCGGCCCAGCCCGTTTGCTCCAACTCAACGCCCCTCGACCt                                                            >  6:63060/1‑135 (MQ=21)
                              ttACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCTGCCCGCGAACACCCGCGCCCCCGGGGTCATCCCACTCCACGCCCCGCGCCCTg                                                           >  4:26927/1‑135 (MQ=18)
                                         tAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCGGACCGCGCACACCCGCCCCCCCGCTTGTATCCCACAAAACGACCCACCCCCCGTACCCAAAccc                                                >  2:286262/1‑133 (MQ=12)
                                                                          ttCCCCGTCACCGCCGCCCCGGCATCCGGGGGCTGCGTGCCGGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGt               <  4:267041/135‑1 (MQ=255)
                                                                             aCCGTCACAGCCGCCACCGCATTCGCGGGGCGCCTGCCCGCGAACACCGGCCCCGCCGGGTTGATACAACTAAACCCCCAGCGACCTGTACGCAAACGTCTCTGCCCACGAATTCCTTGCAGGGCCAATTTTCGt            >  8:2830/1‑135 (MQ=255)
                                                                             aCCGTCACAGCCGCCACCGCATTCGCGGGCCGCCTGCCCGCGAACACCGGCGCAGCCGGGTTCATACAGCTCAACGACCATCGCCCTGTACGCCAACGTCTCTGCCATCGAATGCCTTGCATGGCCAATTGTCGt            >  2:60430/1‑135 (MQ=255)
                                                                              ccGTCACAGCCGCCACCGCATTCGCGGGGAGCCTGCCCGCGCACACCGGCGCCGCCGGGTTGCTCCAACACAACGCCCATCGACCAGTACCCAAAAGGCACTGCCCTCGATTGGCGTGCATGGCCGATTGTCGtg           >  8:42965/1‑134 (MQ=255)
                                                                               cGTCACAGCCGCCACCGCATTCGCGGGCCGCCTGCCCGCGAACAACCGCGCCGCCCGTTTAATACAACAAAAACAACATCGCCCAGTTCCGACAAGGCACAGACAACGATTTCCTTTCATGGCCAATTGTCggtt          >  2:147412/1‑132 (MQ=255)
                                                                                  cacccccgccACCGATTGCGCGGGCAGCGGGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATg       <  6:288724/130‑1 (MQ=255)
                                                                                        cgccACCGCATTCGCGGGCCGCCGGCCCGCGCACACCGGCGCAGCCGGGTTGATCCAACTCAACGCCCATCGACCGGTACGCAAACGTCTCTGGCATCGATTGGCTTGGATTGGCAATTGTGGTTGATGTTATg  >  4:174434/1‑134 (MQ=255)
                                                                                                           |                                                                                                                  
CGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATG  >  ACB122/3491212‑3491433

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: