Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 3,969,743 A→C 51.1% A251A (GCT→GCG CJ019_03750 ← putative protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1433,969,7430AC51.1% ‑3.4 / 7.6 25A251A (GCT→GCGCJ019_03750putative protein
Reads supporting (aligned to +/- strand):  ref base A (3/8);  new base C (9/4);  total (12/13)
Fisher's exact test for biased strand distribution p-value = 9.95e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

GCGGCTGTGGCGCCGACGGCGGGCTGGCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTC  >  ACB143/3969610‑3969874
                                                                                                                                     |                                                                                                                                     
gcggcTGTGGCGCCGACGGCGGGCTGGCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCgcgcgccgc                                                                                                                                     >  4:152103/1‑136 (MQ=255)
   gtgggggCGCCGACGGCGGGCTGGCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCgcgcgcagcgcg                                                                                                                                  <  6:266649/131‑1 (MQ=255)
                      gCTGGCCGGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCGCGACATCCACTTCTTCGATAGCTCCAGCCCCCCGCCAGGCACCGGGGGCGCGGGCAGCGCGCAGCCACACCAGGTCGgcg                                                                                                               <  8:298099/136‑1 (MQ=255)
                         ggcggggggCGTCATACGGCGGCTAGGGCGCCGGGCTAATTCGCCCACCGCCGCTACACCCACTTCTTCGATAGCGGCCCCCCCACGCCCGCCACCGTGGGGGCGCGCAGCGCGCAGCCACACCAGGTCGgcgcgc                                                                                                            <  6:269596/131‑1 (MQ=255)
                         ggcggggggCGCCTTCCGGCGGCCTGGGCCCCGGGCGAATTCGCCCCCCGCCCCCACACCCACTTCTTCGATCGCGCCCCCCCCACGCCCGGCACCGGGGGCGCCGGCAGCCCGCAGCCACACCAGGTCGgcgcgc                                                                                                            <  6:291211/131‑1 (MQ=255)
                         ggCCGGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGAACGCTGCACCCCCACGCCAGCCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGgcgcgc                                                                                                            <  7:371473/136‑1 (MQ=255)
                         ggCCGGGGGCGCCATCCGCCGGCGCGGGCCCCGGGCTAATTCGGCCACCCCCGCCACCCCCACTGCTACGATCGCGCCCCCCCCCCCCCCGGCCCCGGGGGCGCGCGCCGCGCGCAGCCACACCAGGTCGgcgcgc                                                                                                            <  4:42012/136‑1 (MQ=255)
                                  cGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCCGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTc                                                                                                    >  8:99805/1‑135 (MQ=255)
                                  cGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCACCCCCCCCCCGGCGGCGCGCAGGCCGTcc                                                                                                   >  6:319885/1‑136 (MQ=255)
                                       tACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCCGCGCGCACCCCCACCAGGTCGGCGCGCAGGCCCTCCAcccc                                                                                              >  3:265517/1‑136 (MQ=255)
                                       tACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCCCCCGGCCCCCCCCCCCCGGCCGCGCGCAGGCCGGCCAcccc                                                                                              >  7:323164/1‑136 (MQ=255)
                                       tACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCCCCCCGCAGCCCCACCCGGTCGGCGCGCCGGCCCTCCAcccc                                                                                              >  1:236784/1‑136 (MQ=255)
                                       tACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGTGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGcc                                                                                              >  3:199739/1‑136 (MQ=255)
                                        aCGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCCGCCCCACCCGGTCGGCGCGCCGGCCGTCCACGCcg                                                                                             >  5:338228/1‑136 (MQ=255)
                                              gcgAGGGCGCAGGGCTAATGCGCCCCCCGCCGCCACACCCACTTCTTCGACCGCGCCCCCCCCACCCCCGCCACAGGGGGCGCGCCCCGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGgcggcg                                                                                       <  7:347678/136‑1 (MQ=255)
                                                gttGGCGCTGGGCTAATGCGCCCCCCCCCGCCACACCCACTTCTCCCATCGCTCCCCCCCCCCCCCAGCCACTGGGGGCGCGCGCCGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTa                                                                                     <  4:357951/134‑1 (MQ=255)
                                                 aTGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCCGCCCGCCGCCCCCCCCGGGCCGGGGGCCGGCCGGGCACCCCCgcggcgggt                                                                                    >  7:12365/1‑133 (MQ=255)
                                                   ggCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCCGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTagca                                                                                  <  7:243437/136‑1 (MQ=255)
                                                   ggCCCAGGGCTAATTCGCCCCCCCCCTCGACACCCACTTCTTCGATCGCTGCCCCCCCCCCCCCGCCACTGGGGGCGCGCGCCGCGCGAAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTagca                                                                                  <  8:336922/136‑1 (MQ=255)
                                                    gcgcAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCCGCGCGCCGCCCCCCCCGGGCGGCGCGCCGGCCGTCCCCCCCCGCGGCGGaggcg                                                                                 >  2:234690/1‑133 (MQ=255)
                                                        aGGGCGAATTCGGCCACCGCCGCGACATCCACTTCTTCGATCGCTGCAGCCCCCCGCCAGGCACAGGGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCt                                                                             <  4:326489/136‑1 (MQ=255)
                                                                    cccACCGCCTCGACATCCATTTTTTTGTTAGCGGCCGCCCCCGGCCAGGCACTGTGGGCGGGCGCTGCGCGAAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAgg                                                                <  2:19124/135‑1 (MQ=255)
                                                                            cTCGACATCCACTTCTTCGATAGCGGCAGCCCCACGCCCGGCACCGGGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGcc                                                        <  8:89831/136‑1 (MQ=255)
                                                                                 cATCCACTGCTTCTATAGCGGCACCCCCCCGCCAGGCACTGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTc                                                   <  8:200520/136‑1 (MQ=255)
                                                                                                         ccccACGCCAGGCACAGTGGGCGCGCGCCGCGCGCCGCCCCACCAGGTCGGCGCGCCGGCCGTCCACGCCGGCGGCGTTGCCGCGCGCGGTAATCCAGGCCAGCGCCTGGTCATCCAGGGCAATGTAGGGCAAGGc                          >  2:29913/1‑136 (MQ=255)
                                                                                                         ccccACGCCAGGCACAGTGGGCGCGCGCCCCCCGCACCCCCCCCCGGTCGGCGCCCAGGGCCTACCCCCCGCCGGCGGTGCCGCGCGCGGGTAAACAAGCCACCCGCTGGTCATCACAGGGCATGTCGGCCCAggg                          >  8:263864/1‑135 (MQ=255)
                                                                                                                ccAGGCACAGTGGGCGCGCGCCGCGCGCAGCCCCACCCGGTCCGCGCGCACGCCGTCCACGCCGGCGGCGGAGCAGCGCGCGGTAAAACCGGCCAACGCCGGGTGATCCAGGGCAATGTAGGCCAATGCCGGgcgc                   >  2:328560/1‑136 (MQ=255)
                                                                                                                                cgcgccgcgcgcagCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCt   <  8:13147/136‑1 (MQ=255)
                                                                                                                                 gcgcagcgcgcagcCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTc  <  7:99805/136‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                     
GCGGCTGTGGCGCCGACGGCGGGCTGGCCTGGGGCGTCATACGGCGGCGATGGCGCAGGGCGAATTCGGCCACCGCCTCGACATCCACTTCTTCGATAGCTGCAGCCCCACGCCAGGCACAGTGGGCGCGCGCAGCGCGCAGCCACACCAGGTCGGCGCGCAGGCCGTCCACGCCGGCGGCGTAGCAGCGCTCGGTAATCCAGGCAAGCGCCTGGTCATCCAGGGCAATGTCGGCCAAGGCCTGGCGCGCCGCGTGGCAGCGCTC  >  ACB143/3969610‑3969874

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: