Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 2,094,817 T→G 30.5% Q183P (CAG→CCG)  CJ019_01995 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1432,094,8170TG30.5% 26.7 / 5.8 47Q183P (CAG→CCG) CJ019_01995hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (16/14);  new base G (6/7);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCATG  >  ACB143/2094692‑2094947
                                                                                                                             |                                                                                                                                  
ggTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAACGCCCCGCGCGCCGGGCCCGGcgg                                                                                                                          >  2:245335/1‑134 (MQ=255)
     aGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGccgt                                                                                                                     >  2:291734/1‑134 (MQ=255)
            gggCTTGGGCCCGGCCAAACAGTTTCTCAATAAACGTCAGCTGATCTACAAACACCTTCGGCGCCGGCTGCTGGCTCAAATACTCCACCCCGCGCTCTCCGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCa                                                                                                               <  6:375700/135‑1 (MQ=255)
               cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCCCGCCGCGCTCGAAGGCCCGGCGCGCCCGGGCGGGGGGGGCGGTTTGGAGGGc                                                                                                           >  1:351410/1‑136 (MQ=255)
               cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCGGGCCCGGGGGGCCGGTTTGCAGGGc                                                                                                           >  1:65433/1‑136 (MQ=255)
               cATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCACGCCCCCCTGGCCGGGCGGGGCGTTTTTCAaggg                                                                                                           >  5:103342/1‑135 (MQ=255)
                  gggCCCGGGCCAACTGTTTCTCAAGATCCGTCCGCGGCTCTTCAACCACCTTCGGGTCCGGCTGCTTGCCCAAATACTGCCCGCCGCGCGCGACGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGc                                                                                                           <  6:141488/132‑1 (MQ=255)
                     cccGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCCCAAATACTGCACGCCGCGCTCGAAGGCCCGGCCCCCCTGGCCTGGCAGGCCGGTTTGGAAGGCCTgagg                                                                                                     >  5:149473/1‑133 (MQ=255)
                            cAACAGTTTCTCAAGATACGTCAGCTGCTCTTAATAAACCTTCGGCTCCGGCTGCTTTCGCAAATACTGCACGCCGCGCTCGAAGGCCCCCCCCGCCATGCCGGGCGGGGCATTTTTGAGGGCCCGTGGGCCGAgg                                                                                              >  8:10922/1‑136 (MQ=255)
                                 ttttCTCAAGATCCGTCAGCTGCTCTTCATACACCTTCGGCTCGGGCTGCTTGCGCAAATACTCCCCGCCGCGCGCGAAGCCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGgttgtt                                                                                         <  2:109054/135‑1 (MQ=255)
                                 gTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGTCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGgttgtt                                                                                         <  1:388081/136‑1 (MQ=255)
                                  tttCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCGGGGTGGGCGTTTTTCAGGGCCGGTTGGCCGAGGttgttg                                                                                        >  6:297441/1‑136 (MQ=255)
                                   ttCTCAAGAGCCTGCAGCTGCTATAAAAACCCCTGCGGCGCGGGGTGCTGCCCAAAATCCAGCACGCCGCGCTCGACGGCCCGGGGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGttgttg                                                                                        <  4:185765/135‑1 (MQ=255)
                                     cccAAGATACGTCAGCTGCTGTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTgaa                                                                                      <  2:64771/133‑1 (MQ=255)
                                       cAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTTCAGGGCCTGTTGGCCGAGGTTGTTgaaga                                                                                    >  2:271883/1‑135 (MQ=255)
                                       cAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCCCCCAGCCGCGCGGGCCGTTTTGCACGGCCTGGTGGCCGAGGGTGTTtaata                                                                                    >  6:128632/1‑133 (MQ=255)
                                             aCGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCGGCGCCCCATGACGGGCTGGGCGTTTTGCAGGGCCTGTTTGCCGAGGTTGTTGTAGAAcacga                                                                              >  4:403687/1‑131 (MQ=255)
                                                tCAGCTGCTCTTAATACACCTTCGGCTCCGGGTGCTTGCGCAAATACTGCACGCCGCGCTCGACGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATg                                                                          <  6:306626/136‑1 (MQ=255)
                                                ggcgcggctCTTAATACACCTGGGGCGCGGGGTCTGTGCCAAAATCCTCCCCGCCGCGCCCTACGGCCTGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATg                                                                          <  4:275288/130‑1 (MQ=255)
                                                      gctCTTCATACACGTGGGGCTCCGGCTGCTGGCGAAAATTCTCCCCGCGGCGCGCGAAGCCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagcag                                                                    <  1:225445/136‑1 (MQ=255)
                                                       ctctTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCCCGCGCGCCGGGCCGGGCTGGCCTTTTTTCAGGGCCTGTTGGCCGAGGTTGTTGAAGACCACAATATgcagcac                                                                    >  5:256645/1‑134 (MQ=255)
                                                       ctctTCAAACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCCCGCCGCGCTCGACGGCCTGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATgcagcag                                                                    <  1:97816/135‑1 (MQ=255)
                                                            cATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCGGGCCGGGCTGGCCGTTTTGCAGGGCCGGTTGGCCGCGGTTGTTGTAGAACTCTATATGCAGCAGCCCCAg                                                              >  1:45221/1‑136 (MQ=255)
                                                                    ttCGGCTCGGGGTGATTGCCCAAATACTGCACGCCGCGCTCGATGGCCTGGCGGGCTTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGc                                                      <  8:258229/136‑1 (MQ=255)
                                                                      cGGCTCCGTCTGCTTTCGCATATCCTCCACGCGGCGCTCGATGGCGAGGCGCGGCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCg                                                     <  7:245962/135‑1 (MQ=255)
                                                                      cGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCGTTTTGCAGGGCCGGTTGGCCGAGGTTGTTGAAGAGCTCTATAAGCAGCAGCCACAGGATATGGCgg                                                    >  4:114514/1‑136 (MQ=255)
                                                                       ggCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCGGTTGGGCGAGGTTGTTTAAGAACTCTATATTCAGCAGCCCCAGGATATTGCGGa                                                   >  2:41329/1‑136 (MQ=255)
                                                                        gCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGGCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGGAGAACTATATATGCAGCAGCACCAGGATATGGGGGAt                                                  >  7:352092/1‑136 (MQ=255)
                                                                         cTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGGCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCCGCACCAGGATATGGGGGAtt                                                 >  7:320041/1‑136 (MQ=255)
                                                                         cTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCGGGCCGGGCTGGGCTTTTTGCAGGGCCTGTTGGCCGACGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATTGCGGAtt                                                 >  5:114367/1‑136 (MQ=255)
                                                                          gccGGCTGCTTGCCCAAATACTCCCCGCCGCGCTCGACGGCCCGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGAttt                                                <  7:11193/135‑1 (MQ=255)
                                                                          gccGGCTGCTGTCGCAAATAGTGCACGCCGCGCTCGACGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGAttt                                                <  6:306381/135‑1 (MQ=255)
                                                                             ggCTGCTGGCGCAAATCCGGCCCCCCGCGCGCGCCGGCCTGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCt                                             <  6:285499/136‑1 (MQ=255)
                                                                                      cgcAAATACTGCCCGCCGCGCTCGAAGGCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCaa                                     <  8:251050/135‑1 (MQ=255)
                                                                                                  ccgccgcGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGa                        <  2:320485/135‑1 (MQ=255)
                                                                                                    gccgcGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCGTTTTGCAGGGGCTGTTGGGCGAGGTTGTTGTAGAACTCTATATTGAGCAGCACCAGGATATTGCGGATTTGCTTGACCCAATGCTTGGGCTCGa                        >  5:868/1‑134 (MQ=255)
                                                                                                         gcTCGAAGGCCCGGCGCGCCGGGCCGGGCTGGGCGTTTTTCAGGGCCGGTTGGGCGAGGTTGTTGAAGAACTCTATATGCAGCCGCACCAGGATATGGCGGATTTCCTTGAACCAATGCTTGGGCTCGACAGCggg                 >  5:212546/1‑136 (MQ=255)
                                                                                                         gcTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCAGGTGGGCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTGCTTGAACCACTGCTTGGCCTCGACAGCggg                 >  3:92297/1‑136 (MQ=255)
                                                                                                         gcTCGAAGGCCAGGCGCGCCTGGCCGGGCGGGCCTTTTTGGCGGGCCGGGTGGGCGAGGTTGTTGAAGAACTCTATATGCCGCAGCACCAGGATATTGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCggg                 >  8:313355/1‑136 (MQ=255)
                                                                                                         gcTCGAAGGCCACGCCCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGGTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGGGGATTTCCTTGACCAAATGCTTGGCCTCGACAGCggg                 >  1:371520/1‑136 (MQ=255)
                                                                                                             gacggCGAGGCGGGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGATAGCGGGCa               <  1:231612/131‑1 (MQ=255)
                                                                                                             gacggCCTGGCGCGCCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag             <  4:369842/133‑1 (MQ=255)
                                                                                                             cagggCCCGGCGCCGCGGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaag             <  2:157582/133‑1 (MQ=255)
                                                                                                               cggCCCGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCaagaa           <  5:204796/135‑1 (MQ=255)
                                                                                                                        gcgcCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCATg  >  5:340712/1‑136 (MQ=255)
                                                                                                                        gcgcCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCCAGGTTGTTGAATAACTCTATATGCAGCAGCACCAGCATAAGGCGGATTTCCTTTAACCAATGCTTGGCCAAGACAGCCGGCAAGAAACCCTCAtc  >  2:67081/1‑135 (MQ=255)
                                                                                                                        gcgcCTGGCCGCGCTGGCCTATTTGCAGGGCCTGTTGTCCCAGGTTGTTGAAAAACTCTATATGCAGCAGCAACAGGATATTGCGGATTTCCTTGACCCCATTCTTGGCCTCGACAGCGGGCAACAAGACcacaag  >  6:270242/1‑134 (MQ=255)
                                                                                                                             |                                                                                                                                  
GGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAAGCCCTCATG  >  ACB143/2094692‑2094947

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: