Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB143 970,158 C→T 7.1% R20R (CGC→CGT)  cysE → Serine acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB143970,1580CT 7.1% 186.9 / 4.2 56R20R (CGC→CGT) cysESerine acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (25/27);  new base T (2/2);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

TGCGCGGCATCCTCACCGAGACCCAGAAGGTGGTCCGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCTGTGGAGGCGTGATTTCAAATGGCTGGCCCGGCTGGTGTCGAACTTCGGTCGCTGGCTGAC  >  ACB143/970031‑970289
                                                                                                                               |                                                                                                                                   
tGCGCGGCATCCTCACCGAGACCCAGAAGGTGGTCCGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCt                                                                                                                              <  8:90469/135‑1 (MQ=255)
tGCGCGGCATCCTCACCGAGACCCAGAAGGTGGTCCGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCt                                                                                                                              <  3:123054/135‑1 (MQ=255)
 gcgcGGCATCCTCACCGAGACCCAGAAGGTGGTCCGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCtac                                                                                                                            >  3:100389/1‑134 (MQ=255)
 gcgcGGCATCCTCACCGAGACCCAGAAGGTGGTCCGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTc                                                                                                                            >  3:62988/1‑136 (MQ=255)
 gcgcGGCATCCTCACCGAGACCCAGAAGGTGGTCCGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTc                                                                                                                            >  7:289762/1‑136 (MQ=255)
      gCATCCTCACCGAGACCCAGAAGGTGGGCCGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGt                                                                                                                       <  6:377296/136‑1 (MQ=255)
         tCCTCACCGAGACCCAGAAGGTGGTCCGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCt                                                                                                                    >  1:289672/1‑136 (MQ=255)
              aCCGAGACCCAGAAGGTGGTCCGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCAc                                                                                                                 <  6:178888/134‑1 (MQ=255)
                  agaCCCAGAAGGTGGTCCGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTa                                                                                                           >  5:260300/1‑136 (MQ=255)
                  agaCCCAGAAGGTGGTCCGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTa                                                                                                           >  7:318861/1‑136 (MQ=255)
                         gAAGGTGGTCCGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCgg                                                                                                     >  8:187676/1‑135 (MQ=255)
                           aGGTGGTCCGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGt                                                                                                                       >  2:211084/1‑115 (MQ=255)
                           aGGTGGTCCGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGt                                                                                                                       <  1:211084/115‑1 (MQ=255)
                           aGGTGGTCCGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGgcat                                                                                                  <  6:376997/136‑1 (MQ=255)
                             gtggtCCGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGgcatgc                                                                                                <  8:260825/136‑1 (MQ=255)
                                gtCCGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGcatgcatg                                                                                             <  7:187676/136‑1 (MQ=255)
                                   cGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGTAATGCCTTCGAGGTGCTCACCTGCTACCCCGgcat                                                                                                  >  8:342399/1‑128 (MQ=255)
                                   cGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGTAATGCCTTCGAGGTGCTCACCTGCTACCCCGgcat                                                                                                  >  8:337220/1‑128 (MQ=255)
                                   cGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGTAATGCCTTCGAGGTGCTCACCTGCTACCCCGgcat                                                                                                  <  7:342399/128‑1 (MQ=255)
                                   cGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGTAATGCCTTCGAGGTGCTCACCTGCTACCCCGgcat                                                                                                  <  7:337220/128‑1 (MQ=255)
                                   cGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCa                                                                                          >  8:63395/1‑136 (MQ=255)
                                   cGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGACa                                                                                          >  5:46175/1‑136 (MQ=255)
                                   cGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATCCCa                                                                                          >  7:39870/1‑136 (MQ=255)
                                   cGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCCTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCa                                                                                          <  4:300831/136‑1 (MQ=255)
                                   cGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTACGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCa                                                                                          >  4:387604/1‑136 (MQ=255)
                                         gAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGc                                                                                    <  5:21098/136‑1 (MQ=255)
                                             cGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCa                                                                                >  6:377464/1‑136 (MQ=255)
                                             cGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCa                                                                                >  8:321987/1‑136 (MQ=255)
                                              gCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCAc                                                                               <  1:235742/136‑1 (MQ=255)
                                              gCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCAc                                                                               <  3:314422/136‑1 (MQ=255)
                                                aTAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCg                                                                             <  1:134948/136‑1 (MQ=255)
                                                aTAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCg                                                                             <  4:242860/136‑1 (MQ=255)
                                                   aaCGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCt                                                                          >  8:115128/1‑136 (MQ=255)
                                                          ggACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAAt                                                                   <  3:144605/136‑1 (MQ=255)
                                                          ggACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAAt                                                                   <  8:348968/136‑1 (MQ=255)
                                                          ggACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAAt                                                                   <  2:180574/136‑1 (MQ=255)
                                                                 cAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCtgt                                                            <  6:393322/136‑1 (MQ=255)
                                                                   gATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCtgtg                                                           >  8:272946/1‑135 (MQ=255)
                                                                             cGCCTGCGTGAAGATATTCAAAGCGTATTCCAGCGTGACGCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCTGTGGAGGCTTGAtt                                                <  7:240203/136‑1 (MQ=255)
                                                                                   cGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCTGTGGAGGCGTGATTTCAAAt                                          <  7:168112/136‑1 (MQ=255)
                                                                                   cGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCTGTGGAGGCGTGATTTCAAAt                                          <  7:173454/136‑1 (MQ=255)
                                                                                             ttCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCTGTGGAGGCGTGTGTTCAAATGGCTGGCCCg                                >  1:156764/1‑136 (MQ=255)
                                                                                             ttCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCTGTGGAGGCGTGATTTCAAATGGCTGGccc                                 <  2:331903/135‑1 (MQ=255)
                                                                                             ttCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCTGTGGAGGCGTGATTTCAAATGGCTGGccc                                 <  2:305981/135‑1 (MQ=255)
                                                                                             ttCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCTGTGGAGGCGTGATTTCAAATGGCTGGccc                                 <  2:130991/135‑1 (MQ=255)
                                                                                                aaaGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCTGTGGAGGCGTGATTTCAAATGGCTGGCCCGGCt                             >  1:277599/1‑136 (MQ=255)
                                                                                                     gTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCTGTGGAGGCGTGATTTCAAATGGCTGGCCCGGCTGgtgt                        >  6:189403/1‑136 (MQ=255)
                                                                                                         tCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCTGTGGAGGCGTGATTTCAAATGGCTGGCCCGGCTGGTGTCGaa                    >  7:80947/1‑136 (MQ=255)
                                                                                                          ccACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCTGTGGAGGCGTGATTTCAAATGGCTGGCCCGGCTGGGGTCGaa                    >  3:265073/1‑135 (MQ=255)
                                                                                                               gTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCTGTGGAGGCGTGATTTCAAATGGCTGGCCCGGCTGGTGTCGAACTTCgg              <  4:137847/136‑1 (MQ=255)
                                                                                                                 gACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCTGTGGAGGCGTGATTTCAAATGGCTGGCCCGGCTGGTGTCGAACTTCGGTc            >  4:68517/1‑136 (MQ=255)
                                                                                                                   cccGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCTGTGGAGGCGTGATTTCAAATGGCTGGCCCGGCTGGTGTCGAACTTCGGTCgc          >  5:303083/1‑136 (MQ=255)
                                                                                                                     cGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCTGTGGAGGCGTGATTTCAAATGGCTGGCCCGGCTGGTGTCGAACTTCGGTCgctg        <  8:261006/136‑1 (MQ=255)
                                                                                                                     cGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCTGTGGAGGCGTGATTTCAAATGGCTGGCCCGGCTGGTGTAGAACTTAGGTCgctg        >  4:11114/1‑136 (MQ=255)
                                                                                                                       gCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCTGTGGAGGCGTGATTTCAAATGGCTGGCCCGGCTGGTGTCGAACTTCGGTCgctggc      >  6:128160/1‑136 (MQ=255)
                                                                                                                           cACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCTGTGGAGGCGTGATTTCAAATGGCTGGCCCGGCTGGTGTCGAACTTCGGTCGCTGGCTGAc  <  5:237362/136‑1 (MQ=255)
                                                                                                                               |                                                                                                                                   
TGCGCGGCATCCTCACCGAGACCCAGAAGGTGGTCCGGGGCGAGCCGCATAAACGGAAGGACTGACAGATGTTCGAACGCCTGCGTGAAGATATTCAAAGCGTATTCCACCGTGACCCGGCTGCACGCAATGCCTTCGAGGTGCTCACCTGCTACCCCGGCATGCATGCCATCTGGCTGCACCGCCTGGGCAATGCGCTGTGGAGGCGTGATTTCAAATGGCTGGCCCGGCTGGTGTCGAACTTCGGTCGCTGGCTGAC  >  ACB143/970031‑970289

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