Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 3,935,250 C→T 9.5% E137K (GAA→AAA)  CJ019_03720 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1433,935,2500CT 9.5% 132.5 / 2.9 42E137K (GAA→AAA) CJ019_03720hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (19/19);  new base T (2/2);  total (21/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

GCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTAGAACTTGGCGCCCTTGGCACTGC  >  ACB143/3935115‑3935384
                                                                                                                                       |                                                                                                                                      
gCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTc                                                                                                                                        <  3:196927/136‑1 (MQ=255)
                 caccaGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGt                                                                                                                       <  5:137360/136‑1 (MQ=255)
                 caccaGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGt                                                                                                                       <  4:55554/136‑1 (MQ=255)
                      gcagTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCgggg                                                                                                                  <  8:46263/136‑1 (MQ=255)
                             gATCAGGGTGAAGCCTTTGCTGCGTTTCATGTCTGTCACCATTGCCCGTAAGGGATATTGTTAAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTTGCCCTCGCGCGGGTTGTCGGGGctgctg                                                                                                            >  8:303857/1‑135 (MQ=255)
                             gATCAGGGTGAAGCCTTTGATGAGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGACGAGGGCCCTGGAATACACGTCGAACACGTCATTCCCCTAGCGCGGGTTGTCGGGGCTGCAGt                                                                                                           >  8:303858/1‑136 (MQ=255)
                                    gggAAGCCGTTGCTGCGTTTCATGGCTGTCACCATTGCCCGGAAGGGATCTTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGGCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGgc                                                                                                    <  4:166326/134‑1 (MQ=255)
                                           cgttGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAgg                                                                                             <  7:266632/134‑1 (MQ=255)
                                            tttGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGc                                                                                            >  1:52031/1‑136 (MQ=255)
                                                ctgcGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGcccc                                                                                         <  4:61879/135‑1 (MQ=255)
                                                ctgcGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGccc                                                                                          >  5:372670/1‑134 (MQ=255)
                                                   cGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAg                                                                                       <  3:389564/134‑1 (MQ=255)
                                                                   ccATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGGCGTCATcc                                                                     >  6:216387/1‑136 (MQ=255)
                                                                          ccGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTgccg                                                              >  1:95728/1‑136 (MQ=255)
                                                                          ccGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGAGCGCAGGCCCCAGCCGCCTTCGGCTTCCTCCTTTgccg                                                              >  7:97358/1‑136 (MQ=255)
                                                                               aGGGATTTTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGGCGAACACGTCTTTGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGccac                                                          <  7:303858/135‑1 (MQ=255)
                                                                               aGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTTGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGccac                                                          <  7:303857/135‑1 (MQ=255)
                                                                                 ggATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTTCCGCCAACAg                                                       >  8:377923/1‑136 (MQ=255)
                                                                                        gTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTc                                                 <  1:159129/135‑1 (MQ=255)
                                                                                        gTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTc                                                 <  5:382235/135‑1 (MQ=255)
                                                                                        gTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTc                                                 <  6:357272/135‑1 (MQ=255)
                                                                                              gCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCggg                                          >  3:175505/1‑136 (MQ=255)
                                                                                              gCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCggg                                          >  1:378835/1‑136 (MQ=255)
                                                                                               ccTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCCGCGGGTCGTGggggg                                         >  5:220960/1‑135 (MQ=255)
                                                                                                    gCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCCCCAGCGGGTCGTGGGGGGTTcgc                                    >  6:213645/1‑136 (MQ=255)
                                                                                                     ccGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATAcgcc                                   <  1:346730/136‑1 (MQ=255)
                                                                                                     ccGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATAcgcc                                   <  5:341066/136‑1 (MQ=255)
                                                                                                     ccGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATAcgcc                                   <  2:373930/136‑1 (MQ=255)
                                                                                                      cgcgGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCCGCGGGTCGTGCGGGATAcgccg                                  >  5:98472/1‑136 (MQ=255)
                                                                                                      cgcgGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGTTAcgccg                                  >  4:398795/1‑136 (MQ=255)
                                                                                                      cgcgGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGGTAcgccg                                  >  2:366905/1‑136 (MQ=255)
                                                                                                      cgcgGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATAcgccg                                  >  3:298763/1‑136 (MQ=255)
                                                                                                      cgcgGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCGTCGTCTTCATCCTTTGCCGCCCCCCGCGGGTCGTGGGGGATAcgccc                                  >  7:376704/1‑135 (MQ=255)
                                                                                                      cgcgGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGGTGTCGGGGCTGCTGTCGTATGCGCGCAGGCCCCCGCCGCCGTCGTCTTGATCCTTAGCCGCCCCCAGCGGGTCGTGCGGGAtcaggag                                  >  4:86587/1‑130 (MQ=255)
                                                                                                       gcgGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATAcgccgc                                 <  2:393958/136‑1 (MQ=255)
                                                                                                                      cacGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTAGAACtt                   <  4:16028/135‑1 (MQ=255)
                                                                                                                      cacGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTAGAACtt                   <  8:134534/135‑1 (MQ=255)
                                                                                                                      cacGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTAGAACtt                   <  1:164071/135‑1 (MQ=255)
                                                                                                                      cacGTCGAACACGTCTTCGCCCTCGCGCGGGGTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCTCATGAAGTAGAACtt                   <  5:105450/135‑1 (MQ=255)
                                                                                                                            gAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTAGAACTTGGCGccc            >  8:99560/1‑136 (MQ=255)
                                                                                                                              acacGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCCCCAGCGGGTCGTGCGGGATACGCCGCAGGAAGGATCACTTGGCGCCCtt          >  1:276195/1‑136 (MQ=255)
                                                                                                                              acacGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTAGAACTTGGCGCCCtt          <  7:1666/136‑1 (MQ=255)
                                                                                                                                      tCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATACTTTGCCGCCACCAGCGGGTCGTGGGGGGTTCGCCGCAGGGAGTAGCACTTTGCGCCCTTGGCACTgc  >  3:303438/1‑136 (MQ=255)
                                                                                                                                       |                                                                                                                                      
GCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGGCGCGCAGGCCCCAGCCGCCTTCGTCTTCATCCTTTGCCGCCACCAGCGGGTCGTGCGGGATACGCCGCAGGAAGTAGAACTTGGCGCCCTTGGCACTGC  >  ACB143/3935115‑3935384

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