Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 2,094,808 A→C 34.9% L186R (CTG→CGG)  CJ019_01995 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1432,094,8080AC34.9% 20.1 / 5.6 44L186R (CTG→CGG) CJ019_01995hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (14/9);  new base C (9/7);  total (23/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAA  >  ACB143/2094692‑2094938
                                                                                                                    |                                                                                                                                  
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)
                                                                                                                    |                                                                                                                                  
GGTCCAGTACCTGGGCATTGGCCCGGGCCAACAGTTTCTCAAGATACGTCAGCTGCTCTTCATACACCTTCGGCTCCGGCTGCTTGCGCAAATACTGCACGCCGCGCTCGAAGGCCAGGCGCGCCTGGCCGGGCTGGCCTTTTTGCAGGGCCTGTTGGCCGAGGTTGTTGAAGAACTCTATATGCAGCAGCACCAGGATATGGCGGATTTCCTTGACCCAATGCTTGGCCTCGACAGCGGGCAAGAA  >  ACB143/2094692‑2094938

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