Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,718,463 T→G 26.5% Q65P (CAG→CCG)  CJ019_01614 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,718,4630TG26.5% 25.8 / 2.7 45Q65P (CAG→CCG) CJ019_01614hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (11/19);  new base G (8/4);  total (22/23)
Fisher's exact test for biased strand distribution p-value = 9.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

TGATGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGG  >  ACB122/1718335‑1718592
                                                                                                                                |                                                                                                                                 
tgatgaACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGgcgcg                                                                                                                             >  7:26210/1‑135 (MQ=255)
tgatgaACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGgcgcg                                                                                                                             >  7:134574/1‑135 (MQ=255)
tgatcaCCAGCAACTTGCCAATCAGCTCAATGCGCGCACTTACCGACTAGCGGTGGCCAGGCTTGACCAGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGgcgc                                                                                                                              <  1:162683/134‑1 (MQ=255)
    gaACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTCCCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGCCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAg                                                                                                                         <  4:123130/135‑1 (MQ=255)
    gaACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAg                                                                                                                         <  8:134574/135‑1 (MQ=255)
          cAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCCggggccc                                                                                                                   >  1:192402/1‑131 (MQ=255)
          cAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCAc                                                                                                                   >  6:109253/1‑135 (MQ=255)
          cAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCAc                                                                                                                   >  7:278160/1‑135 (MQ=255)
          cAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCCCCAGGGTCAc                                                                                                                   >  4:151087/1‑135 (MQ=255)
               tGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCCCGGccg                                                                                                              >  6:120630/1‑133 (MQ=255)
               tGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCGGGTCCCGGCAg                                                                                                              >  7:130261/1‑135 (MQ=255)
                   ctgcAGCTCAAGGCGCGCACTACCCGCATCGCGGTTGCCAGCGGGGCCCTGCGCGGGCTCATCGGTGAGCTCGACCGGCAGCGGCACGGTCACTTTGACCGGGCAGGCCGGGCGCGCCAGGGTCACGGCAGGGGc                                                                                                          <  8:243556/132‑1 (MQ=255)
                       aGCTCAATGCGCGCACTTCCCGCATCGCGGGGGCCCGGCTGCACCTGCGCGGGGCCATAGGGGAGATGGACCGGCCGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTc                                                                                                       <  4:280485/134‑1 (MQ=255)
                               gcgcgcACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTGCCACAGCGa                                                                                              >  7:102508/1‑135 (MQ=255)
                                cgcgcACTTACGGCATAGCGGGGGCCAGGCTTGCCCGGCGCGGGGTCATCGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGac                                                                                             <  8:288973/135‑1 (MQ=255)
                                     aCTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCCCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATc                                                                                        <  3:151087/135‑1 (MQ=255)
                                                  cGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCCGGGTCACGGCAGGGGCGGCCCTCAGCGCCCCATCCTGCAGTTCCACg                                                                           >  2:182235/1‑135 (MQ=255)
                                                  cGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCAc                                                                            <  2:7453/134‑1 (MQ=255)
                                                  cGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCATTTGACCGGGCCGGCCAGGCGCGCCAGGGTCACGCGAGGGGCGTGCAGAAGCGCCACATGCGGCCGTTGCACg                                                                           >  4:57775/1‑135 (MQ=255)
                                                  cGGTGGCCAGGCTTGACCTGCGCGGGGCCATAGGGGAGATGGAACGGCATCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCAc                                                                            <  8:104358/134‑1 (MQ=255)
                                                          aGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCCGGGGGGCGGCCCGGGCGTCCCTCAGCGACGCCACCTGCCGTTACACGCAGAGGGt                                                                   >  3:29018/1‑135 (MQ=255)
                                                          aGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGACCCCCGGTCACCGCAGGGGCGGGCCACAGCGACAAATCCTGGAGTTGCCCGCTCCgggg                                                                   >  8:10154/1‑134 (MQ=255)
                                                              tggACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCGTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcg                                                               <  4:142999/133‑1 (MQ=255)
                                                                   cggcgcGGGGGCATCGGGCCGATGGAACGGCAGCGGCACGTGGCCTTTGCCCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcggcagg                                                          <  5:107165/133‑1 (MQ=255)
                                                                   cTGCGCGGGGGCATAGGGGAGATGGAACGGCAGCGGCACGTTGCCGTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTggcggcagg                                                          <  8:85112/135‑1 (MQ=255)
                                                                    tGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGGCGGGGGCGCCCGGGGCGCGGCACGGGCGGCCCGCAGCGCCCCCGCCTGCAGTTGCCCGCTCAGGGGggcggcaggc                                                         >  3:261905/1‑135 (MQ=255)
                                                                    tGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCCCGGCCGGGGGGGGCCGCCGCGCCCCCTCCGGCCGTTCCACGCTCCTGGTggcgggaggc                                                         >  2:195909/1‑135 (MQ=255)
                                                                    tGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGGGCGCCCGGGGCGCCGCAGGGGCGGGCCTCAGCCACCCCGCCAGCCGTGGCACCCGCAGGGGggcgggcggc                                                         >  1:2538/1‑130 (MQ=255)
                                                                         ggggTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCCGGGGGGTCCATGAGCGACACCTACTGCAGTTCCACGCGCAGGGGggccggcggcgggcc                                                    >  4:162190/1‑132 (MQ=255)
                                                                         ggggTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCCGGCGCGCCCGGGTGCAGGCAGGGGCGTCCCTCAGCGCCCCATCCTGCAGTTGACCCCTCAGGGTGgcggcacgcggcag                                                    >  8:254155/1‑135 (MQ=255)
                                                                         ggggTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCAGGCGCGCCAGGGTCCCGGCAGGGGGGTCCAACAGCGCCCCATCCAGCAGGTCCCCGCGCAGGGTggcggccggcggccg                                                    >  2:94250/1‑133 (MQ=255)
                                                                            gTCATTGTTGAGATGGAACGGCTGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGc                                                 <  6:302788/135‑1 (MQ=255)
                                                                              cATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCCGCAGCGACCCATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCCGGGCaa                                               >  7:148110/1‑135 (MQ=255)
                                                                                tAGGTGAGATGGAACGGCAGCGGCCCGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAAt                                              <  7:289484/134‑1 (MQ=255)
                                                                                   gTGAGAGGGACCGGCAGCGGCACGTTGCCTTGGACCGGGCCGGCCTGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATgcgc                                          <  5:203665/135‑1 (MQ=255)
                                                                                           gAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTa                                  <  8:186608/135‑1 (MQ=255)
                                                                                                gCAGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCCGGGTCCCGGCCGGGGCCGCCCTCCGCGACCCAACCTGCAGGTACCCGCACAGGGTGGGGGGCGGCGGGCGGGCCAAGCGCGGCAGGTagaaga                             >  4:139535/1‑135 (MQ=255)
                                                                                                  aGCGGCACGTTGCCTTTGACCGGGCCGGCCGGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACt                           <  1:182235/135‑1 (MQ=255)
                                                                                                           ttGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTcc                  <  2:192402/135‑1 (MQ=255)
                                                                                                            tGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCa                 <  8:278160/135‑1 (MQ=255)
                                                                                                                   gACCGGGCCGGCCTGGCGCGCCAGGGGCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTg            <  8:303414/133‑1 (MQ=255)
                                                                                                                         gccggccgGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGtt    <  5:20438/135‑1 (MQ=255)
                                                                                                                          ccggcgGGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCGACGCTCTGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTtg   <  1:286727/135‑1 (MQ=255)
                                                                                                                          ccggccTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTtg   <  5:106507/135‑1 (MQ=255)
                                                                                                                          ccggccTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTtg   <  8:61575/135‑1 (MQ=255)
                                                                                                                           cggccGGGCGCGCCCGGGGCCCGGCAGGGGCGGCCATCAGCGCCACATCCAGCCGTTACACGCTCCCGGTGGGGGCAGGCGGCCGGGCCATGCGCGGCCGGGAGAACACGTCACCGTCCAGGCTGGCCGGGGtgg  >  6:260209/1‑135 (MQ=255)
                                                                                                                                |                                                                                                                                 
TGATGAACAGCAACTTGCCATTCAGCTCAATGCGCGCACTTACCGCATAGCGGTGGCCAGGCTTGACCTGCGCGGGGTCATAGGTGAGATGGAACGGCAGCGGCACGTTGCCTTTGACCGGGCCGGCCTGGCGCGCCAGGGTCACGGCAGGGGCGTCCATCAGCGACACATCCTGCAGTTCCACGCTCAGGGTGGCGGCAGGCGGCAGGGCAATGCGCTGCAGGTAGAAGACTTCACCGTCCAGGCTGGCCTGGTTGG  >  ACB122/1718335‑1718592

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