Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 5,171,928 A→G 6.9% Y15H (TAC→CAC)  CJ019_04861 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1225,171,9280AG 6.9% 182.6 / 2.0 58Y15H (TAC→CAC) CJ019_04861hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (20/33);  new base G (2/2);  total (23/35)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACCGTCGCGTTGCCAGCGTTCGACGGCAATGAAGAAGGCATCATGCATGGCCTCCTCGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAGGCGGTCCAGGCTTTCTCCCCAGCCCTGGTGGAAGCCCATTTCTTCATGGGCCCGGC  >  ACB122/5171805‑5172061
                                                                                                                           |                                                                                                                                     
tACCGTCGCGTTGCCAGCGTTCGACGGCAATGAAGAAGGCATCATGCATGGCCTCCTCGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGa                                                                                                                            <  2:395864/135‑1 (MQ=255)
 aCCGTCGCGTTGCCAGCGTGCGACGGCAATGAAGAAGGCAGCATGCATGGCCTCCGCGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGa                                                                                                                            <  6:71585/134‑1 (MQ=255)
    gTCGCGTTGCGAGCGGTCGACGGCAATGAAGACGGCATCATGCATGGCGGCCTCGGCCAGGTCGAAACCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTGGACTTc                                                                                                                        <  6:119714/135‑1 (MQ=255)
     tCGCGTTGCCAGCGTTCGACGGCAATGAAGAAGGCATCATGCATGGCCTCCTCGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCg                                                                                                                       <  5:350854/135‑1 (MQ=255)
     tCGCGTTGCCAGCGTTCGACGGCAATGAAGAAGGCATCATGCATGGCCTCCTCGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCg                                                                                                                       <  4:159558/135‑1 (MQ=255)
     tCGCGTTGCCAGCGTTCGACGGCAATGAAGAAGGCATCATGCATGGCCTCCTCGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCg                                                                                                                       <  6:329417/135‑1 (MQ=255)
     tCGCGTTGCCAGCGTTCGACGGAAATGAAGAAGGCATCATGCATGGCCTCCTCGGCCAGGTCTAAATCACCCAGCAAGCGGATCTGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCg                                                                                                                       <  1:102201/135‑1 (MQ=255)
          ttGCCAGCGTTCGACGGCAATGAAGAAGGCATCATGCATGGCCTCCTCGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGgcacg                                                                                                                  <  4:279976/135‑1 (MQ=255)
                gCGTTCGCCGGCATTGAAGAAGGCATCATGCATGGCCTCCTCGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGGATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTg                                                                                                            <  7:417131/135‑1 (MQ=255)
                gCGTTCGACGGCAATTAAGAAGGCATCATGCATGGCCTCCTCGGCCAGGTCGAAATCACCCAGCAAGCGGAGCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTg                                                                                                            <  6:449170/135‑1 (MQ=255)
                gCGTTCGACGGCAATGAAGAAGGCATCATGCATGGCCTCCGCGGCCAGGTCGAAATCACCCAGCAAGCGGAGCAGTGGCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTg                                                                                                            <  2:170698/135‑1 (MQ=255)
                       aCGGCAATGAAGAAGGCATCATGCATGGCCTCCTCGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGt                                                                                                     <  6:106899/135‑1 (MQ=255)
                                 agaagGCATCATGCATGGCCTCCTCGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATc                                                                                            <  2:316600/134‑1 (MQ=255)
                                            tgcatgGCCTCCTCGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATc                                                                                 <  4:114985/134‑1 (MQ=255)
                                             gcatgGCCTCCTCTGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTCCACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCg                                                                               >  5:296018/1‑135 (MQ=255)
                                             gcatgGCCTCCTCGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTGCACAGCCTCGACTTCGGCACGCACCTGGGCCa                                                                                                       <  3:444144/111‑1 (MQ=255)
                                             gcatgGCCTCCTCGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTGCACAGCCTCGACTTCGGCACGCACCTGGGCCa                                                                                                       <  3:451649/111‑1 (MQ=255)
                                             gcatgGCCTCCTCGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTGCACAGCCTCGACTTCGGCACGCACCTGGGCCa                                                                                                       >  4:444144/1‑111 (MQ=255)
                                             gcatgGCCTCCTCGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTGCACAGCCTCGACTTCGGCACGCACCTGGGCCa                                                                                                       >  4:451649/1‑111 (MQ=255)
                                             gcatgGCCTCCTCGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGa                                                                                                                            <  7:453850/90‑1 (MQ=255)
                                             gcatgGCCTCCTCGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGa                                                                                                                            >  8:453850/1‑90 (MQ=255)
                                             gcatgGCCTCCTCGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCg                                                                               >  6:449973/1‑135 (MQ=255)
                                             gcatgGCCTCCTCGGCCAGGTCGAAATCACCCAGCAAGAGTATCAGTGTCGCCAGTATGCGCAGCGACTCCTGACGGTACACAGCCTCTACTTCCGCACGCACCATGGCCACTGAGGCCCACAGGAGGGCAGccc                                                                               >  1:204622/1‑134 (MQ=255)
                                             gcatgGCCTCCTCGGCAAGGTCGAAATCAACCAGAAAGCGGATCAGTGTCGCCAGTATTCGCCGCGACTCCTGCCGGTACACAGCATCGACTTCGGCACGCAAATGGGCCACTGCGGCCATCAGTCGGGCATACg                                                                               >  3:37748/1‑135 (MQ=255)
                                                 tgCCTCCTCGGCCAGGTAGAAATCACCCAGCAAGCTGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCATTGCGGCCATCAGTCTGGCATCAGCTgg                                                                           <  8:87961/134‑1 (MQ=255)
                                                   cctcctCGGCCAGGTCGAAATCGCCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTggg                                                                          >  7:254741/1‑134 (MQ=255)
                                                   cctcctCGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGt                                                                         <  3:209066/135‑1 (MQ=255)
                                                    ctcctcGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCCGTGCGGCCATCAGTCGGGCATCCGCTGGGTc                                                                        >  4:31715/1‑135 (MQ=255)
                                                    ctcctcGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGACGGTACACAGCCTCGACTTCGGCACGACCCTGGGCCACTGAGGACATAACTCCGGCCTCCGCTGGGTc                                                                        >  7:177433/1‑135 (MQ=255)
                                                        tcGGCCAGGTCGAAATCACCCAGCCAGCGGATCAGTGTCGCCAGTATGCGCTGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTcac                                                                      <  5:143467/133‑1 (MQ=255)
                                                          ggCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGCCTCCTGCCGGTACACAGCCTCGACTACGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCaccacc                                                                  >  3:270559/1‑135 (MQ=255)
                                                          ggCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTcac                                                                      <  5:103120/131‑1 (MQ=255)
                                                          ggCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTcac                                                                      >  6:103120/1‑131 (MQ=255)
                                                                gTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTccg                                                              >  6:167579/1‑132 (MQ=255)
                                                                   gAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCCCGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAgg                                                         >  1:320754/1‑135 (MQ=255)
                                                                    aaaTCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCCCGCACCTGGGCCCGTGCGGCCCTCAGTCGGGCCTCCGCTGGGTCCCCCCCTCCCCCAGGc                                                        >  5:330841/1‑135 (MQ=255)
                                                                    aaaTCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCCCCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCCCCTCCCCCCGGc                                                        >  5:336885/1‑135 (MQ=255)
                                                                      aTCGCCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAGGCgg                                                      <  8:254741/135‑1 (MQ=255)
                                                                             aGCAAGCGGATCAGTGTCGCAAGTCTGCTCCGCGACTCCTGCAGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAGGCGGTCCAGGc                                               <  2:204622/135‑1 (MQ=255)
                                                                                  ggggATAAGTGCCGCCAGTAGGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAGGCGGTCCAGGCTTtct                                          <  7:396135/133‑1 (MQ=255)
                                                                                  gCGGATCAGTGTCGCCATTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAGGCGGTCCAGGCTTtct                                          <  6:237399/135‑1 (MQ=255)
                                                                                  gCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCCCCTCCACCAGGCGGTCCAGGCGTtct                                          >  8:460805/1‑135 (MQ=255)
                                                                                  gCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAGGCGGTCCAGGCTTtct                                          <  3:445355/135‑1 (MQ=255)
                                                                                  gCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAGGCGGTCCAGGCTTtct                                          <  6:208354/135‑1 (MQ=255)
                                                                                  gCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAGGCGGTCCAGGCTTtct                                          <  3:418632/135‑1 (MQ=255)
                                                                                  gCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAGGCGGTCCAGGCTTtct                                          <  5:462590/135‑1 (MQ=255)
                                                                                  gCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAGGCGGTCCAGGCTTtct                                          <  8:361414/135‑1 (MQ=255)
                                                                                  gCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAGGCGGTCCAGGCTTtct                                          <  3:113472/135‑1 (MQ=255)
                                                                                       tCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAGGCGGTCCAGGCTTTCTCCCCa                                     <  4:270559/135‑1 (MQ=255)
                                                                                        cAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAGGCGGTCCAGGCTTTCTCCCCAg                                    <  7:10899/135‑1 (MQ=255)
                                                                                                   gTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAGGCGGTCCAGGCTTTCTCCCCAGCCCTGGTGGaa                         >  5:137601/1‑135 (MQ=255)
                                                                                                   gTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAGGCGGTCCAGGCTTTCTCCCCAGCCCTGGTGGaa                         >  7:316995/1‑135 (MQ=255)
                                                                                                            gcgACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAGGCGGTCCAGGCTTTCTCCCCAGCCCTGGTGGAAGCCCATttc                >  7:307599/1‑135 (MQ=255)
                                                                                                              gACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAGGCGGTCCAGGCTTTCTCCCCCGCCCTGGTGGAAGCCCATttct               >  2:48648/1‑134 (MQ=255)
                                                                                                                cTCCTGCCGGTACACAGCCTCGACTTCGGCACGCCCCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAGGCGGTCCCGGCTTTCTCCCCAGCCCTGGTGGAAGCCCATTTCTTCa            >  7:77274/1‑135 (MQ=255)
                                                                                                                    tGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAGGCGGTCCAGGCTTTCTCCCCAGCCCTGGTGGAAGCCCATTTCTTCATgg         <  7:239952/134‑1 (MQ=255)
                                                                                                                        ggTACACAGCCTCGACTTCGGCACGCCCCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCCCCTCCCCCAGGCGGTCCCGGCTTTCTCCCCAGCCCTGGTGGAAGCCCAGTTCTTCATGGGCCCg    >  4:194215/1‑135 (MQ=255)
                                                                                                                          tACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAGGCGGTCCAGGCTTTCTCCCCAGCCCTGGTGGAAGCCCATTTCTTCATGGGCCCGGc  <  3:207541/135‑1 (MQ=255)
                                                                                                                           |                                                                                                                                     
TACCGTCGCGTTGCCAGCGTTCGACGGCAATGAAGAAGGCATCATGCATGGCCTCCTCGGCCAGGTCGAAATCACCCAGCAAGCGGATCAGTGTCGCCAGTATGCGCCGCGACTCCTGCCGGTACACAGCCTCGACTTCGGCACGCACCTGGGCCAGTGCGGCCATCAGTCGGGCATCCGCTGGGTCACCACCTCCACCAGGCGGTCCAGGCTTTCTCCCCAGCCCTGGTGGAAGCCCATTTCTTCATGGGCCCGGC  >  ACB122/5171805‑5172061

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