Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,403,593 C→T 7.1% S37L (TCG→TTG)  acuI → putative acryloyl‑CoA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,403,5930CT 7.1% 173.7 / 3.6 56S37L (TCG→TTG) acuIputative acryloyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base C (22/30);  new base T (2/2);  total (24/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

AATTATTTTCCAATCCATCATGCACATGAGGACCACTTATGCAGGCGTTACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTGTACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACT  >  NC_000913/3403446‑3403741
                                                                                                                                                   |                                                                                                                                                    
aaTTATTTTCCAATCCATCATGCACATGAGGACCACTTATGCAGGCGTTACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCg                                                                                                                                                     >  1:134823/1‑149 (MQ=255)
aaTTATTTTCCAATCCATCATGCACATGAGGACCACTTATGCAGGCGTTACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCg                                                                                                                                                     >  1:250417/1‑149 (MQ=255)
aaTTATTTTCCAATCCATCATGCACATGAGGACCACTTATGCAGGCGTTACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCg                                                                                                                                                     >  4:330032/1‑149 (MQ=255)
 attattTTCCAATCCATCATGCACATGAGGACCACTTATGCAGGCGTTACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGa                                                                                                                                                    <  1:361529/149‑1 (MQ=255)
     ttttCCAATCCATCATGCACATGAGGACCACTTATGCAGGCGTTACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCt                                                                                                                                                >  3:288430/1‑149 (MQ=255)
      tttCCAATCCATCATGCACATGAGGACCACTTATGCAGGCGTTACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTg                                                                                                                                               <  2:250417/149‑1 (MQ=255)
       ttCCAATCCATCATGCACATGAGGACCACTTATGCAGGCGTTACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGa                                                                                                                                              <  4:295450/149‑1 (MQ=255)
        tCCAATCCATCATGCACATGAGGACCACTTATGCAGGCGTTACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGaa                                                                                                                                             <  1:494458/149‑1 (MQ=255)
          cAATCCATCATGCACATGAGGACCACTTATGCAGGCGTTACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACt                                                                                                                                           >  3:290837/1‑149 (MQ=255)
            aTCCATCATGCACATGAGGACCACTTATGCAGGCGTTACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACtat                                                                                                                                         >  1:382676/1‑149 (MQ=255)
                    tGCACATGAGGACCACTTATGCAGGCGTTACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATgc                                                                                                                                 >  4:34535/1‑149 (MQ=255)
                    tGCACATGAGGACCACTTATGCAGGCGTTACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATgc                                                                                                                                 <  3:34535/149‑1 (MQ=255)
                           gAGGACCACTTATGCAGGCGTTACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCg                                                                                                                          <  4:567010/149‑1 (MQ=255)
                                ccACTTATGCAGGCGTTACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTAc                                                                                                                     >  2:399492/1‑149 (MQ=255)
                                 cACTTATGCAGGCGTTACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTAcc                                                                                                                    >  1:444398/1‑149 (MQ=255)
                                  aCTTATGCAGGCGTTACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCg                                                                                                                   >  4:557936/1‑149 (MQ=255)
                                      aTGCAGGCGTTACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTaa                                                                                                               <  1:451965/149‑1 (MQ=255)
                                         cAGGCGTTACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAggg                                                                                                            >  4:541933/1‑149 (MQ=255)
                                               ttACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTaa                                                                                                               <  4:305922/140‑1 (MQ=255)
                                               ttACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTaa                                                                                                               >  3:305922/1‑140 (MQ=255)
                                                  cTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatcat                                                                                                   <  1:98839/149‑1 (MQ=255)
                                                   ttttAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatcatc                                                                                                  >  4:23415/1‑149 (MQ=255)
                                                   ttttAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatcatc                                                                                                  <  1:403693/149‑1 (MQ=255)
                                                    tttAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATcc                                                                                                 <  3:557936/149‑1 (MQ=255)
                                                    tttAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATcc                                                                                                 >  3:444095/1‑149 (MQ=255)
                                                     ttAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCg                                                                                                <  4:417250/149‑1 (MQ=255)
                                                              cagGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTcc                                                                                       >  3:132402/1‑149 (MQ=255)
                                                               agGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCg                                                                                      >  3:557360/1‑149 (MQ=255)
                                                                      cAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTc                                                                               >  4:334742/1‑149 (MQ=255)
                                                                      cAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCAGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTc                                                                               >  2:520332/1‑149 (MQ=255)
                                                                        aaaCTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCt                                                                             <  1:399492/149‑1 (MQ=255)
                                                                        aaaCTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCt                                                                             <  4:221735/149‑1 (MQ=255)
                                                                        aaaCTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCt                                                                             <  1:568120/149‑1 (MQ=255)
                                                                                gCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGa                                                                     <  4:544592/149‑1 (MQ=255)
                                                                                       tACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGcc                                                              <  3:330032/149‑1 (MQ=255)
                                                                                           gcctctGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGAtttt                                                                 <  1:236630/140‑1 (MQ=255)
                                                                                           gACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGAtttt                                                                 >  2:236630/1‑142 (MQ=255)
                                                                                               ctGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTGt                                                      >  3:489063/1‑149 (MQ=255)
                                                                                                     gAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGcc                                                              <  2:331745/135‑1 (MQ=255)
                                                                                                     gAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGcc                                                              >  1:331745/1‑135 (MQ=255)
                                                                                                                    gAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTGTACGCACCAGCGAAGATCcgcg                                 <  2:134823/149‑1 (MQ=255)
                                                                                                                       ggCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTGTACGCACCAGCGAAGATCCGCGttt                              <  4:290837/149‑1 (MQ=255)
                                                                                                                         cGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTGTACGCACCAGCGAAGATCCGCGTTTTc                            <  4:272451/149‑1 (MQ=255)
                                                                                                                         cGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTGTACGCACCAGCGAAGATCCGCGTTTTc                            <  4:444095/149‑1 (MQ=255)
                                                                                                                          gATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTGTACGCACCAGCGAAGATCCGCGTTTTCa                           <  4:489063/149‑1 (MQ=255)
                                                                                                                           aTGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTGTACGCACCAGCGAAGATCCGCGTTTTCAt                          <  2:244294/149‑1 (MQ=255)
                                                                                                                             gTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTGTACGCACCAGCGAAGATCCGCGTTTTCATGc                        <  4:395938/149‑1 (MQ=255)
                                                                                                                               cACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTGTACGCACCAGCGAAGATCCGCGTTTTCATGCCg                      <  3:334742/149‑1 (MQ=255)
                                                                                                                                aCGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTGTACGCACCAGCGAAGATCCGCGTTTTCATGCCgg                     <  2:77463/149‑1 (MQ=255)
                                                                                                                                   gTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTGTACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCa                  >  1:379957/1‑149 (MQ=255)
                                                                                                                                       aTGTTCACTGGTTGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatcat                                                                                                   >  1:472715/1‑64 (MQ=255)
                                                                                                                                       aTGTTCACTGGTTGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatcat                                                                                                   <  4:270591/64‑1 (MQ=255)
                                                                                                                                       aTGTTCACTGGTTGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatcat                                                                                                   <  2:472715/64‑1 (MQ=255)
                                                                                                                                       aTGTTCACTGGTTGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAatcat                                                                                                   >  3:270591/1‑64 (MQ=255)
                                                                                                                                             aCTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTGTACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTa        <  4:288430/149‑1 (MQ=255)
                                                                                                                                                  tCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTGTACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCAc   <  2:382676/149‑1 (MQ=255)
                                                                                                                                                   cGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAAATTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTGTACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCAct  <  4:152306/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AATTATTTTCCAATCCATCATGCACATGAGGACCACTTATGCAGGCGTTACTTTTAGAACAGCAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAAGATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTGTACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACT  >  NC_000913/3403446‑3403741

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