Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,491,323 T→G 51.2% intergenic (‑685/+36) 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,3230TG51.2% ‑2.3 / 9.0 25intergenic (‑685/+36)CJ019_03306/CJ019_03307hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (4/7);  new base G (10/2);  total (16/9)
Fisher's exact test for biased strand distribution p-value = 3.61e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

CTCCTACATCGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACACCC  >  ACB122/3491199‑3491453
                                                                                                                            |                                                                                                                                  
cTCCTACATCGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGcaca                                                                                                                         >  8:87821/1‑136 (MQ=255)
    tACATCGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGGCGCCGGCCCCCGCACACCg                                                                                                                      >  2:339986/1‑135 (MQ=21)
              ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGGCGCCGGGCCGCGCACACCCGCGCCGCCggg                                                                                                           >  7:52024/1‑135 (MQ=18)
              ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGGCGCCGGCCCCCGCACACCCCCGCCGCCgg                                                                                                            >  5:327310/1‑135 (MQ=18)
              ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGGCCCCCGCCCCCCGACACCCGCgccgcccc                                                                                                            >  1:149583/1‑133 (MQ=18)
              ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCGGCCCGCGGACACCCGCGCCGCCgg                                                                                                            >  3:42641/1‑135 (MQ=21)
              ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCCCAGCCgg                                                                                                            >  4:124180/1‑135 (MQ=255)
              ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCCCGGGCCGCCGGCCCGGGAACACCCGCGCCGCCgg                                                                                                            >  6:321132/1‑135 (MQ=21)
              ggtggtGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCCCGGGCCGCCGGCCCCCCCAACCCGGCACAGcccgg                                                                                                           >  7:175024/1‑133 (MQ=18)
                   gggcGCCTCGTTGTCGGCGCGGGTTTCCCCGCAAATTGCGCGGCCCTGCCGGCCACTATCCCTGGCCATGCCCCGACCCCCCCCCCCCCGCATGCGCGGGCAGCCGCCCCGCCAACACCGGCGCAGCCGGTTTCa                                                                                                       <  4:53660/132‑1 (MQ=255)
                                           ttACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCGGCCCGCGCACACCCGCGCAGCCCGGTTCATCCAACTCAACCCCCAACGACCt                                                                                >  2:216770/1‑134 (MQ=21)
                                                      tAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGGCGCCGGCCCCCGCACACCCCCCCCGCCGCGTTGCTCCCACTCAACGCCCCGCGCCCGGTGCGCCAACgc                                                                    >  1:263829/1‑134 (MQ=9)
                                                      tAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCCGCCGGCCCGCGCACACCGGCGCCGCCGGGTTGATACCACTCCACCCCCCACGCCCCGTACGCAAACGt                                                                    >  5:106149/1‑135 (MQ=18)
                                                         gCCGGACCTGCCGGCCACTATCCCTGGTCATGCCCCGACACAGCGCCCACCGCATTGGCGGGCAGCCTGCCGGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGtctc                                                                 <  7:251489/135‑1 (MQ=255)
                                                            ggACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCCGGCCGCCGGCCCGCGCACACCCGCGCCGCCGGGTTCATCCAACACAACGACCCTCCCCCGGTACGCCAACGTCTCTGcg                                                             >  7:216343/1‑135 (MQ=21)
                                                                      gCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCCGGGCGCCTGCCCCCGAACACCGGCGCAGCCGGTGTACTACCACACCACGCCCATCGCCCGGTACGCCAACGTCTCGGCCATCGATTGc                                                    >  6:188142/1‑135 (MQ=25)
                                                                       ccACTATCCCTGGGCATACCCCGTCACAGCCGCCACCGCATCCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCt                                                  <  8:52024/136‑1 (MQ=255)
                                                                        cACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCCCGGGCAGCCTCCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCtt                                                 <  7:329438/136‑1 (MQ=255)
                                                                        acagaccccTGGTCTTCCCCCGGCCCCGCCGCCCCCGCATTCGGGGGCCGCCTCCCGGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCt                                                  <  3:31446/129‑1 (MQ=255)
                                                                                       ttCACCGTCACCCCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGt                                   <  3:238239/135‑1 (MQ=255)
                                                                                          cccGCCACCGCCGCCCCCGCATGCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCg                                 <  4:42641/133‑1 (MQ=255)
                                                                                          aCCGTCACAGCCGCCACCGCATTCGCGGGCCGCCAGCCCGCGAACACCCGCGCAGCCGGGTTGATCCCACACAACCCCCCACCACCGGTACGCCAACCTCTCTGCCCTCGATTTCCATGGATGGCCAAtagctgt                                >  4:31446/1‑129 (MQ=255)
                                                                                                 caGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACCTCTCTGCCATCGATTGCCTTGCATGGCCCAATGTGGTTTATCTa                         >  6:106883/1‑135 (MQ=255)
                                                                                                        ccccGCATTCGGGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCg                 <  7:33061/134‑1 (MQ=255)
                                                                                                                       cAGCCGGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACAccc  <  8:98465/136‑1 (MQ=255)
                                                                                                                            |                                                                                                                                  
CTCCTACATCGACCGGTGGTGTCGCCACGTTGTAGGAGCTGGTTTACCCGCGAATAGGCCGGACCTGCCGGCCACTATCCCTGGTCATTCACCGTCACAGCCGCCACCGCATTCGCGGGCAGCCTGCCCGCGAACACCGGCGCAGCCGGTTTCATCCAACTCAACGACCATCGACCTGTACGCAAACGTCTCTGCCATCGATTGCCTTGCATGGCCAATTGTCGTTGATGTAATGATCCGGCTTGAAGCACACCC  >  ACB122/3491199‑3491453

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: