Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 4,119,178 G→A 7.1% A94A (GCC→GCT)  CJ019_03877 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1434,119,1780GA 7.1% 175.4 / 3.6 57A94A (GCC→GCT) CJ019_03877hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (13/39);  new base A (2/2);  total (15/42)
Fisher's exact test for biased strand distribution p-value = 2.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-02

ACGCACGCTTGGCCGCCGCACGGTCGGGGTGGCCGATGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATT  >  ACB143/4119045‑4119302
                                                                                                                                     |                                                                                                                            
acgcacgcTTGGCCGCCGCACGGTCGGGGTGGCCGATGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGc                                                                                                                            <  8:245514/136‑1 (MQ=255)
  gcacgcTTGGCCGCCGCACGGTCGGGGTGGCCGATGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCAt                                                                                                                          >  2:285741/1‑136 (MQ=255)
   cacgcTTGGCCGCCGCACGGTCGGGGTGGCCGATGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATc                                                                                                                         <  1:1734/136‑1 (MQ=255)
   cacgcTTGGCCGCCGCACGGTCGGGGTGGCCGATGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATc                                                                                                                         <  4:93824/136‑1 (MQ=255)
        ttGGCCGCCGCACGGTCGGGGTGGCCGATGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCcac                                                                                                                     >  4:169045/1‑135 (MQ=255)
           gccgccGCACGGTCGGGGTGGCCGATGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGc                                                                                                                 <  1:248941/136‑1 (MQ=255)
           gccgccGCACGGTCGGGGTGGCCGATGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGc                                                                                                                 <  6:385336/136‑1 (MQ=255)
           gccgccGCACGGTCGGGGTGGCCGATGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGc                                                                                                                 <  6:199792/136‑1 (MQ=255)
           gccgccGCACGGTCGGGGTGGCCGATGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGc                                                                                                                 <  2:10189/136‑1 (MQ=255)
           gccgccGCACGGTCGGGGTGGCCGATGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGc                                                                                                                 <  2:116536/136‑1 (MQ=255)
           gccgccGCACGGTCGGGGTGGCCGATGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGc                                                                                                                 <  3:122544/136‑1 (MQ=255)
           gccgccGCACGGTCGGGGTGGCCGATGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGc                                                                                                                 <  2:39228/136‑1 (MQ=255)
                cgcACGGTCGGGGTGGCCGATGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGcggtcg                                                                                                            <  4:91786/136‑1 (MQ=255)
                         ggggTGGCCGATGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGc                                                                                                    <  5:11375/135‑1 (MQ=255)
                          gggTGGCCGATGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCt                                                                                                  >  6:163958/1‑136 (MQ=255)
                             tGGCCGATGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGt                                                                                               <  1:317701/136‑1 (MQ=255)
                             tGGCCGATGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGt                                                                                               <  5:210471/136‑1 (MQ=255)
                                  gATGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCcgcg                                                                                          <  4:40369/136‑1 (MQ=255)
                                    tGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCgg                                                                                         >  2:283101/1‑135 (MQ=255)
                                    tGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGt                                                                                        <  4:242452/136‑1 (MQ=255)
                                                  aTGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACg                                                                          <  6:180755/136‑1 (MQ=255)
                                                  aTGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACg                                                                          <  3:343145/136‑1 (MQ=255)
                                                   tGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACAGACAATGGCATCGAACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGTAACg                                                                          <  3:395779/135‑1 (MQ=255)
                                                        ggAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAAt                                                                     >  6:170066/1‑135 (MQ=255)
                                                          aGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGcc                                                                  <  7:85954/136‑1 (MQ=255)
                                                          aGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGcc                                                                  <  6:217778/136‑1 (MQ=255)
                                                              aTCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTca                                                               <  1:302120/135‑1 (MQ=255)
                                                              aTCTGCCACGAAATCCCCCATTTGTCCCTGAACCAGCCGCAAACGCTGGCCTGGCGGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTca                                                               <  2:260371/135‑1 (MQ=255)
                                                               tCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTca                                                               >  7:198726/1‑134 (MQ=255)
                                                                 tGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGt                                                           >  8:316956/1‑136 (MQ=255)
                                                                       gAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTg                                                      >  2:263340/1‑135 (MQ=255)
                                                                        aaaTCCCCCATTTTTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGaa                                                    <  7:64043/136‑1 (MQ=255)
                                                                          aTCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAgg                                                  <  4:391595/136‑1 (MQ=255)
                                                                             ccccATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCg                                               <  5:208803/136‑1 (MQ=255)
                                                                             ccccATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCg                                               <  2:15265/136‑1 (MQ=255)
                                                                               ccATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCggg                                             <  1:28036/136‑1 (MQ=255)
                                                                                cATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGc                                            >  7:22174/1‑136 (MQ=255)
                                                                                cATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGc                                            >  5:309335/1‑136 (MQ=255)
                                                                                    tGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATAGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGccac                                         >  8:342153/1‑135 (MQ=255)
                                                                                    tGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATAGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGccac                                         >  8:347385/1‑135 (MQ=255)
                                                                                         cTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGa                                   <  5:213481/136‑1 (MQ=255)
                                                                                         cTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGa                                   <  4:26460/136‑1 (MQ=255)
                                                                                         cTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGa                                   <  8:145951/136‑1 (MQ=255)
                                                                                         cTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACAGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCTTTGa                                   <  8:145956/136‑1 (MQ=255)
                                                                                         cTGCACAAGCCGCAAACGGTGGCCTCGCCTTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTGCGCGGTACTGACCTGGAACGATTATGCCTCAGAGTGCTTGAAGGCCGGGCCACCATTGa                                   <  5:308238/136‑1 (MQ=255)
                                                                                            cACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAg                                  <  3:367504/134‑1 (MQ=255)
                                                                                             aCCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGcc                               >  7:101767/1‑136 (MQ=255)
                                                                                               cAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCTCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCa                             <  8:154571/136‑1 (MQ=255)
                                                                                                aGCCGAAAACGCGGTCCTAGCCTTCGTTACCGACAATTTCATCACACAAGCGTGCGGTCTCTGCCTGGTCCGAGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCAc                            <  8:372634/136‑1 (MQ=255)
                                                                                                   cGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGc                          >  4:182028/1‑135 (MQ=255)
                                                                                                    gCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGc                          >  4:280069/1‑134 (MQ=255)
                                                                                                    gCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCaa                        <  1:48718/136‑1 (MQ=255)
                                                                                                      ataCGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAgg                      <  6:33357/134‑1 (MQ=255)
                                                                                                      aaccGCTGGCCGCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGGCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAgg                      <  5:293362/133‑1 (MQ=255)
                                                                                                                          ttACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAAtt  <  2:333171/136‑1 (MQ=255)
                                                                                                                          ttACCGACAATAGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAAtt  <  7:347385/136‑1 (MQ=255)
                                                                                                                          ttACCGACAATAGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAAtt  <  7:342153/136‑1 (MQ=255)
                                                                                                                                     |                                                                                                                            
ACGCACGCTTGGCCGCCGCACGGTCGGGGTGGCCGATGGCTTCGATGAGGATGCGTGGAGTGATCTGCCACGAAATCCCCCATTTGTCCCTGCACCAGCCGCAAACGCTGGCCTCGCCGTCGTTACCGACAATGGCATCCCACAAGCGGTCGGTCTCTGCCTGGTCCGCGGTACTGACCTGGAACGAAAATGCCTCACAGTGCTTGAAGGCCGGGCCACCATTGAGGCCCACGCAAGGGATACCCATTACCGTGAATT  >  ACB143/4119045‑4119302

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