Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 3,214,279 G→A 20.5% I474I (ATC→ATT)  CJ019_03017 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1433,214,2790GA20.5% 83.7 / 15.9 39I474I (ATC→ATT) CJ019_03017hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (10/21);  new base A (4/4);  total (14/25)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.13e-02

CGGTGGCGTCGGCACCTTGGGTGGTGGCCTGGGTGTCCTGTTGCAGCCTCTCGGTGTCTTGCTGGATCGAACGGGAAAGGCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAATGCTTGACCTGGCCTTTTCGCCGTTCAGGTTCATCGCGTGCAAGACACTGCTGGCGCTGTTTTGCAG  >  ACB143/3214152‑3214407
                                                                                                                               |                                                                                                                                
cggcggCGTCGGCACCTTGGGTGGTGGCCTGGGTGTCCTGTTGCAGCCTCTCGGTGTCTTGCTGGATCGAACGGGAAAGGCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGc                                                                                                                           <  1:26309/135‑1 (MQ=255)
 ggTGGCGTCGGCACCTTGGGTGGTGGCCTGGGTGTCCTGTTGCAGCCTCTCGGTGTCTTGCTGGATCGAACGGGAAAGGCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACGCATCAATCTGGCGa                                                                                                                         >  5:71307/1‑136 (MQ=255)
    ggCGTCGGCACCTTGGGTGGTGGCCTGGGTGTCCTGTTGCAGCCTCTCGGTGTCTTGCTGGATCGAACGGGAAAGGCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACg                                                                                                                      <  5:235135/136‑1 (MQ=255)
     gCGTCGGCACCTTGGGTGGTGGCCTGGGTGTCCTGTTGCAGCCTCTCGGTGTCTTGCTGGATCGAACGGGAAAGGCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACgg                                                                                                                     <  3:173273/136‑1 (MQ=255)
           gCACCTTGGGTGGTGGCCTGGGTGTCCTGTTGCAGCCTCTCGGTGTCTTGCTGGATCGAACGGGAAAGGCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCgg                                                                                                               >  6:93788/1‑136 (MQ=255)
                  gggTGGTGGCCTGGGTGTCCTGTTGCAGCCTCTCGGTGTCTTGCTGGATCGAACGGGAAAGGCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATAGATCTGGCGAACGGCCCCGGTAATCGc                                                                                                        >  3:189825/1‑136 (MQ=255)
                            cTGGGTGTCCTGTTGCAGCCTCTCGGTGTCTTGCTGGATCGAACGGGAAAGGCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTc                                                                                              <  1:83061/136‑1 (MQ=255)
                            cTGGGTGTCCTGTTGCAGCCTCTCGGTGTCTTGCTGGATCGAACGGGAAAGGCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTc                                                                                              >  3:130615/1‑136 (MQ=255)
                                           cAGCCTCTCGGTGTCTTGCTGGATCTAACGGGAAAGGCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTc                                                                                <  1:270173/135‑1 (MQ=255)
                                                tctcGGTGTCTTGCTGGATCGAACGGGAAAGGCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCAATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAgcg                                                                          >  2:219780/1‑136 (MQ=255)
                                                  tcGGTGTCTTGCTGGATCGAACGGGAAAGGCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAgcg                                                                          <  6:221672/134‑1 (MQ=255)
                                                  tcGGTGTCTTGCTGGATCGAACGGGAAAGGCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAgcg                                                                          <  4:25178/134‑1 (MQ=255)
                                                  tcGGTGTCTTGCTGGATCGAACGGGAAAGGCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCAATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAgcg                                                                          <  7:135079/134‑1 (MQ=255)
                                                     gtgtCTTGCTGGATCGAACGGGAAAGCCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATAGCCTCCAGGGTCTGCGTGGCATGATCAGAGCGCTGa                                                                      >  5:269966/1‑135 (MQ=255)
                                                      tgtCTTGCTGGATCGAACGGGAAAGGCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCAATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAAt                                                                    <  5:35065/136‑1 (MQ=255)
                                                            gCTGGATCGAACGGGAAAGGCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAATGCtt                                                                <  8:22610/134‑1 (MQ=255)
                                                             cTGGATCGAACGGGAAAGGCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCAATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAATGCTTGAc                                                             >  7:129559/1‑136 (MQ=255)
                                                                      aCGGGAAAGGCCGATCTGCTCGTTTGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGATCGCTGAATGCTTGACCTGGCCttt                                                    <  3:8356/136‑1 (MQ=255)
                                                                        gggAAAGGCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATATGGCGAACGGCCCCGGTAATAGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCGGAATGCTTTACCTGGCCTTTTc                                                  >  6:183714/1‑136 (MQ=255)
                                                                        gggAAAGGCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTGTAACTCATCGATCTGGCGAACGGCCCCGGTTATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAATGCTTGACCTGGCCTTTTc                                                  >  6:180453/1‑136 (MQ=255)
                                                                           aaaGGCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAATGCTTGACCTGGCCTTTTCGcc                                               <  5:87265/136‑1 (MQ=255)
                                                                                      tGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAATGCTTGACCTGGCCTTTTCGCCGTTCAGGTTCa                                    >  2:222624/1‑136 (MQ=255)
                                                                                       gCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAATGCTTGACCTGGCCTTTTCGCCGTTCAGGTTCAt                                   >  3:123509/1‑136 (MQ=255)
                                                                                        cTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAATGCTTGACCTGGCCTTTTCGCCGTTCAGGt                                       <  6:35107/131‑1 (MQ=255)
                                                                                        cTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAATGCTTGACCTGGCCTTTTCGCCGTTCAGGt                                       <  8:82590/131‑1 (MQ=255)
                                                                                        cTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAATGCTTGACCTGGCCTTTTCGCCGTTCAGGTTc                                     <  5:250478/133‑1 (MQ=255)
                                                                                        cTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCAATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAATGCTTGACCTGGCCTTTTCGCCGTTCAGGTTc                                     <  7:95715/133‑1 (MQ=255)
                                                                                        cTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCAATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAATGCTTGACCTGGCCTTTTCGCCGTTCAGGTTc                                     <  1:219780/133‑1 (MQ=255)
                                                                                           gTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAATGCTTGACCTGGCCTTTTCGCCGTTCAGGTTCATcgc                                <  5:190161/135‑1 (MQ=255)
                                                                                            ttGGTGAAACGGGCGATGCCGGCGTTTAACTCATCAATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTCGCATGTTCAGAGCGCTGAATGCTTGACCTGGCCTTTTCGCCGTTCAGGTTCATCGCGt                              >  2:164464/1‑136 (MQ=255)
                                                                                                tGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAATGCTTGACCTGGCCTTTTCGCCGTTCAGGTTCATCGCGTGCaa                          <  7:81578/136‑1 (MQ=255)
                                                                                                     cGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGCTCTGCGTGGCATGTTCAGAGCGCTGAATGCTTGACCTGGCCTTTTCGCAGTTCAGGTTCATCGCGTGCAAGacac                     >  7:95462/1‑136 (MQ=255)
                                                                                                        gCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAATGCtt                                                                >  4:167289/1‑90 (MQ=255)
                                                                                                        gCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAATGCtt                                                                <  3:167289/90‑1 (MQ=255)
                                                                                                          gATGGCGGCGTTTAACTCATCGATCTGGCGAACGGACCCGGTAATAGCCTCGAGTTTCTGCTTGTCATGTTCATAGCGCTGAATGCTTGACCTGTCCTTTTCGCCGTTCTGGTTCATCGCGTGCAAGACActgctt                <  1:180207/136‑2 (MQ=255)
                                                                                                                   ggttAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAATGCTTGACCTGGCCTTTTCGCCGTTCAGGTTCATCGCGTGCAAGACACTGCTGGCGCTGttt       <  8:166470/134‑1 (MQ=255)
                                                                                                                   gTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAATGCTTGACCTGGCCTTTTCGCCGTTCAGGTTCATCGCGTGCAAGACACTGCTGGCGCTGttt       <  4:130615/136‑1 (MQ=255)
                                                                                                                        aCTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAATGCTTGACCTGGCCTTTTCGCCGTTCAGGTTCATCGCGTGCAAGACACTGCTGGCGCTGTTTTGCAg  <  1:135805/136‑1 (MQ=255)
                                                                                                                        aCTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAATGCTTGACCTGGCCTTTTCGCCGTTCAGGTTCATCGCGTGCAAGACACTGCTGGCGCTGTTTTGCAg  <  3:69034/136‑1 (MQ=255)
                                                                                                                               |                                                                                                                                
CGGTGGCGTCGGCACCTTGGGTGGTGGCCTGGGTGTCCTGTTGCAGCCTCTCGGTGTCTTGCTGGATCGAACGGGAAAGGCCGATCTGCTCGTTGGTGAAACGGGCGATGCCGGCGTTTAACTCATCGATCTGGCGAACGGCCCCGGTAATCGCCTCCAGGGTCTGCGTGGCATGTTCAGAGCGCTGAATGCTTGACCTGGCCTTTTCGCCGTTCAGGTTCATCGCGTGCAAGACACTGCTGGCGCTGTTTTGCAG  >  ACB143/3214152‑3214407

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