Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 2,729,080 T→A 9.1% C980S (TGC→AGC)  CJ019_02576 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1432,729,0800TA 9.1% 136.1 / 2.4 43C980S (TGC→AGC) CJ019_02576hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (27/12);  new base A (2/2);  total (29/14)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GACGCGCAGGCCAGGTTATCCATGCTGGAGCAGCAATATGGGGTGCTCAAAAGGCACTTGACCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGGCGCTTCATGCTGACCCAGACCATCAGCCGCTTGCGCCATCTGCCGGCCCTGCCCGCGCAAGTGGGCTTTCCCCACG  >  ACB143/2728945‑2729209
                                                                                                                                       |                                                                                                                                 
gACGCGCAGGCCAGGTTATCCATGCTGGAGCAGCAATATGGGGTGCTCAAAAGGCACTTGACCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAAct                                                                                                                                   >  6:8655/1‑136 (MQ=255)
    cgcAGGCCAGGTTATCCATGCTGGAGCAGCAATATGGGGTGCTCAAAAGGCACTTGACCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAActgctg                                                                                                                               <  3:299859/136‑1 (MQ=255)
                          ggAGCAGCAATATGGGGTGCTCAAAAGGCACTTGACCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGACCTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTg                                                                                                         >  7:264696/1‑136 (MQ=255)
                           gAGCAGCAATATGGGGTGCTCAAAAGGCACTTGACCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTg                                                                                                         >  4:32918/1‑135 (MQ=255)
                               agcaATATGGGGTGCTCAAAAGGCACTTGACCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGATGTCTTGTGCCGGAACTGCAAgg                                                                                                    >  3:31997/1‑136 (MQ=255)
                               agcaATATGGGGTGCTCAAAAGGCACTTGACCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTCGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAgg                                                                                                    >  8:74929/1‑136 (MQ=255)
                                    tatGGGGTGCTCAAAAGGCACTTGACCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAg                                                                                               <  5:8655/136‑1 (MQ=255)
                                              tCAAAAGGCACTTGACCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGa                                                                                     >  7:103831/1‑136 (MQ=255)
                                              tCAAAAGGCACTTGACCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGa                                                                                     >  1:97271/1‑136 (MQ=255)
                                              tCAAAAGGCACTTGACCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGa                                                                                     >  2:202916/1‑136 (MQ=255)
                                              tCAAAAGGCACTTGACCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGa                                                                                     >  2:31185/1‑136 (MQ=255)
                                              tCAAAAGGCACTTGACCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGa                                                                                     >  5:49864/1‑136 (MQ=255)
                                              tCAAAAGGCACTTGACCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGATGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAATCGGCACTGGa                                                                                     >  6:240599/1‑136 (MQ=255)
                                              tCAAAAGGCACTTGACCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACAGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGa                                                                                     >  5:222909/1‑136 (MQ=255)
                                              tCAAAAGGCACTTGACCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACAGCTGGCGCTGTCTTGTGCCGGAAATACAAGGGCAAGCGGCACTGGa                                                                                     >  5:222913/1‑136 (MQ=255)
                                               cAAAAGGCACTTGACCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTTTGCCGGAACTGCAAGGGCAAGCGGCACTGGAt                                                                                    >  5:240546/1‑136 (MQ=255)
                                                aaaaGGCACTTGACCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATg                                                                                   >  1:41893/1‑136 (MQ=255)
                                                 aaaGGCACTTGACCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGa                                                                                  >  3:89104/1‑136 (MQ=255)
                                                      cACTTGACCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGgc                                                                             <  7:35075/136‑1 (MQ=255)
                                                          tGACCCAGTGGGCACGCAGTGCCTTGCTGTCA‑CAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGGCGCtt                                                                         >  7:240458/1‑135 (MQ=255)
                                                            aCCCAGTTGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACAGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGGCGCTTc                                                                        <  6:222913/135‑1 (MQ=255)
                                                            aCCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGGCGCTTCa                                                                       <  3:32918/136‑1 (MQ=255)
                                                            aCCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACAGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGGCGCTTc                                                                        <  6:222909/135‑1 (MQ=255)
                                                               cAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGGCGCTTCATGc                                                                    >  5:37985/1‑136 (MQ=255)
                                                                 gTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGGCGCTTCATGCt                                                                   >  3:257583/1‑135 (MQ=255)
                                                                      cACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACTGGCGCTTCATGCTg                                                                  <  7:153621/131‑1 (MQ=255)
                                                                      cACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACTGGCGCTTCATGCTg                                                                  >  8:153621/1‑131 (MQ=255)
                                                                         gCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGGCGCTTCATGCTGACCCAGAc                                                          >  1:185181/1‑136 (MQ=255)
                                                                                  tGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGc                                                                                           <  5:263254/94‑1 (MQ=255)
                                                                                  tGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGc                                                                                           >  6:263254/1‑94 (MQ=255)
                                                                                       tCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGGCGCTTCATGCTGACCCAGACCATCAGCCGCTTgc                                            <  3:154727/136‑1 (MQ=255)
                                                                                          ccAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGGCGCTTCATGCTGACCCAGACCATCAGCCGCTTgcgc                                          >  7:252204/1‑135 (MQ=255)
                                                                                          ccAGAACTGGCCGCCCGGCGGGAGTTCAGGAAGGGGCTGAGCAACTGCTGGCGCTGGCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGGCGCTTCATGCTGACCCAGACCATCAGCCGCTTGCGcc                                         <  4:210706/136‑1 (MQ=255)
                                                                                            aGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGCCGGGCGCTTCATGCTGACCCAGACCATCAGCCGCTTGCGCCAt                                       >  1:89114/1‑136 (MQ=255)
                                                                                            aGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGGCGCTTCATGCTGACCCAGACCATCAGCCGCTTGCGCCAt                                       >  1:21112/1‑136 (MQ=255)
                                                                                                      gccCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGa                                                                                  <  7:54577/83‑1 (MQ=255)
                                                                                                      gccCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGa                                                                                  >  8:54577/1‑83 (MQ=255)
                                                                                                      gccCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGGCGCTTCATGCTGACCCAGCCCATCAGCCGCTTGCGCCATCTGCCGGccc                             >  1:305735/1‑136 (MQ=255)
                                                                                                      gccCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGGCGCTTCATGCTGACCCAGACCATCAGCCGCTTGCGCCATCTGCCGGccc                             >  4:243312/1‑136 (MQ=255)
                                                                                                      gccCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGGCGCTTCATGCTGACCCAGACCATCAGCCCCTTGCGCCATCTGCCGGccc                             >  4:174063/1‑136 (MQ=255)
                                                                                                       cccGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGGCGCTTCATGCTGACCCAGACCATCAGCCGCTTGCGCCATCTGCCGGCCCt                            <  8:240458/136‑1 (MQ=255)
                                                                                                                aGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGGCGCTTCATGCTGACCCAGACCATCAGCCGCTTGCGCCATCTGCCGGCCCTGCCCgcgc                    <  6:283763/135‑1 (MQ=255)
                                                                                                                     aGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGGCGCTTCATGCTGACCCAGACCATCAGCCGCTTGCGCCATCTGCCGGCCCTGCCCGCGCAAGTgg              >  4:224837/1‑136 (MQ=255)
                                                                                                                      ggAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGGCGCTTCATGCTGACCCAGACCATCAGCCGCTTGCGCCATCTGCCGGCCCTGCCCGCGCAAGTggg             <  5:103536/136‑1 (MQ=255)
                                                                                                                                 aTCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGGCGCTTCATGCTGACCCAGACCATCAGCCGCTTGCGCCATCTGCCGGCCCTGCCCGCGCAAGTGGGCTTTCCCCACg  <  6:240546/136‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                 
GACGCGCAGGCCAGGTTATCCATGCTGGAGCAGCAATATGGGGTGCTCAAAAGGCACTTGACCCAGTGGGCACGCAGTGCCTTGCTGTCACCAGAACTGGCCGCCCGGCGTGAGTTCAGGAAGGGGCTGATCAACTGCTGGCGCTGTCTTGTGCCGGAACTGCAAGGGCAAGCGGCACTGGATGACGGGCGCTTCATGCTGACCCAGACCATCAGCCGCTTGCGCCATCTGCCGGCCCTGCCCGCGCAAGTGGGCTTTCCCCACG  >  ACB143/2728945‑2729209

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