Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 2,696,522 C→T 8.5% Q208* (CAG→TAG)  CJ019_02551 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1222,696,5220CT 8.5% 148.0 / 3.8 47Q208* (CAG→TAG) CJ019_02551hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (17/26);  new base T (2/2);  total (19/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

GGCGTTATCAAAGCCATTTGCAATCGAAGCTTGAAGGCGAGGGGGGGGCCGTCGAGAGCACCCTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACAACGCCGTGTTGAAGTGTCACACCCTGCGGCTGTTGGGGAAGGAGGTCATTGGCGCACTGGTGATCGAAG  >  ACB122/2696388‑2696648
                                                                                                                                      |                                                                                                                              
ggCGTTATCAAAGCCATTTGCAATCGAAGCTTGAAGGCGAGGGGGGGGCCGTCGAGCGCACCCTGGAGGAGGCGTTGCGGGCCAACCTGGCGCTTGCGGCCATCGCCCCCAGGCTCAAGGATGAAATCCATGAGc                                                                                                                                >  2:335556/1‑135 (MQ=255)
        cAAAGCCATTTGCAATCGAAGCTTGAAGGCGAGGGGGGGGGCGTCGAGAGCACCCTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTg                                                                                                                         <  3:150132/134‑1 (MQ=255)
        aaaaGCCATTTGCAATCTAAGCTGGCAGGCGAGGGGGGGGCCGTCGAGAGCACCCTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTg                                                                                                                         <  1:71468/133‑1 (MQ=255)
            ggtaTTTGCAATCGACGCTTGATGGCGAGGGGGGGGCCTTCGAGAGCAGCCTGGAGGAGGCGTTGAGGGCGAACCGGGCGCTTGCTGCAATCGCCTCTAGGATCAAGGATGATATCGATGAGCAGACCTGTCAGTg                                                                                                                   <  8:262919/133‑1 (MQ=255)
            gCCATTTGCAATCGAAGCTTGAAGGCGAGGGGGGGGCCGTCGAGAGCACCCTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTg                                                                                                                   <  8:329610/136‑1 (MQ=255)
            gCCATTTGCAATCGAAGCTTGAAGGCGAGGGGGGGGCCGTCGAGAGCACCCTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTg                                                                                                                   <  8:324000/136‑1 (MQ=255)
              cATTTGCAATCGAAGCTTGAAGGCGAGGGGGGGGCCGTCGAGAGCACCCTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGCCCTGTCAGTgg                                                                                                                  >  6:127922/1‑135 (MQ=255)
                   gCAATCGAAGCTTGAAGGCGAGGGGGGGGCCGTCGAGAGCACCCTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCa                                                                                                              <  3:301308/134‑1 (MQ=255)
                    cAATCGAAGCTTGAAGGCGAGGGGGGGGCCGTCGAGAGCACCCTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCaa                                                                                                             >  1:37597/1‑134 (MQ=255)
                     aaTCGAAGCTTGAAGGCGAGGGGGGGGCCGTCGAGAGCACCCTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAAc                                                                                                            <  2:254374/134‑1 (MQ=255)
                          aaGCTTGAAGGCGAGGGGGGGGCCGTCGAGAGCACCCTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGc                                                                                                       <  3:195169/134‑1 (MQ=255)
                              ttGAAGGAGAGGGGGTGGCCGGCGAGAGCACCCTTGAGGAGGCGTTGAGGGCGAACCTGGCGCTTGGTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTgg                                                                                                  <  3:34075/135‑1 (MQ=255)
                                  cggCGAGGGGGGGGCCGTCGAGAGCACCCTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGcac                                                                                              <  3:135781/134‑1 (MQ=255)
                                   ggCGAGGGGGGGGCCGTCGAGAGCACACTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGcaca                                                                                             <  3:56581/135‑1 (MQ=255)
                                          ggggggCCGTCGAGAGCACCCTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTTCTGCAATCGCCGCCAGGATCAAGGATGAAAACGATGAGTATACCTGTCAGTGGATAAACCAGGCGGTGGGCACAGCGTCt                                                                                       >  6:148795/1‑134 (MQ=255)
                                          ggggggCCGTCGAGAGCACCCTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTGCTGCCATCGCCGCAAGGATCAAGGATGAAATCGATGAGTAGACCGGTCAGTGGCTCAACCAGGCGGTGGGCACCGCGTCtt                                                                                      >  6:155719/1‑135 (MQ=255)
                                          ggggggCCGTCGAGAGCACCCTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTGCTGCCATCGCCGCAAGGATCAAGGATGAAATAGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTcat                                                                                      >  3:153622/1‑133 (MQ=255)
                                          ggggggCCGTCGAGAGCACCCTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCtt                                                                                      >  6:204001/1‑135 (MQ=255)
                                             gggCCGTCGAGAGCACCCTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTcc                                                                                   <  3:347056/135‑1 (MQ=255)
                                               gCCGTCGAGAGCACCCTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTg                                                                                 <  5:204001/135‑1 (MQ=255)
                                                 cGTCGAGAGCACCCTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGcc                                                                               >  4:195372/1‑135 (MQ=255)
                                                        agCACCCTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACa                                                                       <  7:350715/136‑1 (MQ=255)
                                                         gCACCCTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACa                                                                       <  3:327547/135‑1 (MQ=255)
                                                         gCACCCTGGAGGAGGCGGTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGGCTTTCCTGCCAGCCGACa                                                                       <  4:21542/135‑1 (MQ=255)
                                                                ggaggaGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGCGCAGACCTGTCAGTGGATAAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACAACGCCgt                                                                >  3:52333/1‑135 (MQ=255)
                                                                ggaggaGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACAACGCCgtg                                                               >  7:314770/1‑136 (MQ=255)
                                                                ggaggaGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATAGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCTGTGGGCACAGCGTCTTTCCTGCCAGCCGACAACGCCgt                                                                >  6:378924/1‑135 (MQ=255)
                                                                       gCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCAGCCAGCCGACAACGCCGTGTTGAAg                                                         >  1:23363/1‑135 (MQ=255)
                                                                         gTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACAACGCCGTGTTGAAgtg                                                       <  4:353169/135‑1 (MQ=255)
                                                                         gTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTATTTCCTGCCAGCCGACAACGCCGTGTTGAAgtg                                                       >  1:341656/1‑135 (MQ=255)
                                                                              gggCCAACCTGGCGCTTGCTGCAAACGCCGCAAGGATCAAGGAAGAAATCGATGAGCAGCCCGGTCAGGGGATAAACCAGCCCGGGGGCACAGCGGCGTTCAGGCCCGCCGCAAACGCCGGGATGAAGTGTcaca                                                  >  6:79212/1‑135 (MQ=255)
                                                                                 ccAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACAACGCCGTGTTGAAGTGTCACAccc                                               <  4:59203/135‑1 (MQ=255)
                                                                                  cAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACAACGCCGTGTTGAAGTGTCACACCCTg                                             <  8:293839/136‑1 (MQ=255)
                                                                                  cAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACAACGCCGTGTTGAAGTGTCACACCCTg                                             <  8:314770/136‑1 (MQ=255)
                                                                                     ccTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACAACGCCGTGTTGAAGTGTCACACCCTGCg                                           >  6:139843/1‑135 (MQ=255)
                                                                                       tGGCGCTTGCTGCAAACCCCGCCAGGATCAAGGATGAAATCGATGAGCACACCTGTCAGTGGATCAACCACGCGGTGGGCACAGCCGCTTTCCTGCCCGCCGCCAACGCCGTGGTGAAGTGTAACACCCAGCGGc                                         >  2:358507/1‑135 (MQ=255)
                                                                                                ctgcAAGCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACAACGCCGTGTTGAAGTGTCACACCCTGCGGCTGTTgggg                                 <  1:335556/134‑1 (MQ=255)
                                                                                                 tgcAATCGCCGCAAGGATCAAGGATTAAATCGATGAGCAGACCTGTCAGTGGATCAACCACGCGGTGGGCCCAGCGTCTTTCCTGCCAGCCGACAACGCCGTGTTGAAGTGTCACACCCTGCGGCTGTTGGGGaa                               >  6:335946/1‑135 (MQ=255)
                                                                                                  gcAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACAACGCCGTGTTGAAGTGTCACACCCTg                                             <  2:288594/120‑1 (MQ=255)
                                                                                                  gcAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACAACGCCGTGTTGAAGTGTCACACCCTg                                             >  1:288594/1‑120 (MQ=255)
                                                                                                    aaTCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACAACGCCGTGTTGAAGTGTCACACCCTGCGGCTGTTGGGGAagga                            <  2:37597/135‑1 (MQ=255)
                                                                                                       cgccgcAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACAACGCCGTGTTGAAGTGTCACACCCTGCGGCTGTTGGGGAAGGAGGt                         <  1:9592/135‑1 (MQ=255)
                                                                                                                 atcaaggatGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACAACGCCGTGTTGAAGTGTCACACCCTGCGGCTGTTGGGGAAGGAGGTCATTGGCGCa               <  2:31957/135‑1 (MQ=255)
                                                                                                                  tcaaggatGAAATCGATGAGTAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACAACGCCGTGTTGAAGTGTCACACCCTGCGGCTGTTGGGGAAGGAGGTCATTGGCGCAc              <  5:155719/135‑1 (MQ=255)
                                                                                                                  tcaaggatGAAATCGATGAGTAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACAACGCCGTGTTGAAGTGTCACACCCTGCGGCTGTTGGGGAAGGAGGTCATTGGCGCAc              <  5:148795/135‑1 (MQ=255)
                                                                                                                            aaTCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACAACGCCGTGTTGAAGTGTCACACCCTGCGGCTGTTGGGGAAGGAGGTCATTGGCGCACTGGTGATc      >  1:280045/1‑133 (MQ=255)
                                                                                                                            aaTCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACAACGCCGTGTTGAAGTGTCACACCCTGCGGCTGTTGGGGAAGGAGGTCATTGGCGCACTGGTGATCg     <  3:195372/134‑1 (MQ=255)
                                                                                                                             aTCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACAACGCCGTGTTGAAGTGTCACACCCTGCGGCTGTTGGGGAAGGAGGTCATTGGCGCTCTGGTGATCGAAg  >  7:74262/1‑136 (MQ=255)
                                                                                                                                      |                                                                                                                              
GGCGTTATCAAAGCCATTTGCAATCGAAGCTTGAAGGCGAGGGGGGGGCCGTCGAGAGCACCCTGGAGGAGGCGTTGAGGGCCAACCTGGCGCTTGCTGCAATCGCCGCAAGGATCAAGGATGAAATCGATGAGCAGACCTGTCAGTGGATCAACCAGGCGGTGGGCACAGCGTCTTTCCTGCCAGCCGACAACGCCGTGTTGAAGTGTCACACCCTGCGGCTGTTGGGGAAGGAGGTCATTGGCGCACTGGTGATCGAAG  >  ACB122/2696388‑2696648

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