Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 353,359 A→G 8.6% Q218Q (CAA→CAG)  prpE → propionate‑‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913353,3590AG 8.6% 115.4 / 4.2 49Q218Q (CAA→CAG) prpEpropionate‑‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base A (20/21);  new base G (2/2);  total (25/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-01

AAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCA  >  NC_000913/353213‑353507
                                                                                                                                                  |                                                                                                                                                    
aaaaTCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAacac                                                                                                                                                   >  1:747698/1‑150 (MQ=255)
aaaaTCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAacac                                                                                                                                                   <  2:310656/150‑1 (MQ=255)
     cATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCACCCCATCgg                                                                                                                                              >  2:1293903/1‑150 (MQ=255)
      aTTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGgc                                                                                                                                             <  1:1343542/150‑1 (MQ=255)
         cccTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGgcgcg                                                                                                                                          >  1:621225/1‑150 (MQ=255)
         cccTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGgcgcg                                                                                                                                          >  1:38252/1‑150 (MQ=255)
            tataAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCgg                                                                                                                                       <  1:761672/150‑1 (MQ=255)
            tataAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCgg                                                                                                                                       >  1:1149700/1‑150 (MQ=255)
            tataAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCgg                                                                                                                                       >  1:939986/1‑150 (MQ=255)
             ataAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCggg                                                                                                                                      <  2:1205540/150‑1 (MQ=255)
                aaaaaTTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCATCACATCGGCGCGCGGGTAc                                                                                                                                   >  1:243105/1‑150 (MQ=255)
                aaaaaTTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTAc                                                                                                                                   >  2:1051099/1‑150 (MQ=255)
                aaaaaTTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTAc                                                                                                                                   <  2:749721/150‑1 (MQ=255)
                 aaaaTTGCGCGACGATGCGATAAGGCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAAGGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTAcc                                                                                                                                  <  1:1217622/150‑1 (MQ=255)
                        cTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGgtggcg                                                                                                                           >  1:506264/1‑150 (MQ=255)
                           gagtaTGCGATAAGTCAGGCGAAGGACCAGCCACGCCATGTTTTGCTGGTTGATAGCGGGCTGGCGAAAATGGCCCGCGGCAGCGGGCGGGATGGCGATTTCGCGTCGGTTCGCCATCAACACTTCGGCGCGCGGGTACCGgtggcgtgg                                                                                                                        <  1:855154/146‑1 (MQ=255)
                           gacgaTGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGgtggcgtgg                                                                                                                        >  2:1128048/1‑150 (MQ=255)
                           gacgaTGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTGTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGgtggcgtgg                                                                                                                        <  1:26077/150‑1 (MQ=255)
                             cgttGCGATAAGTCAGGCGCATCCCCAGCCACGCCATGTTTCGCTGGTGGATCGCGGGCTGGCGAAAAGGGCGCGCGTCAGCGGGCGGGATGTCGATATCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCt                                                                                                                      <  1:258569/147‑1 (MQ=255)
                               aTGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTgg                                                                                                                    <  2:1100808/150‑1 (MQ=255)
                                tGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGa                                                                                                                   >  2:284688/1‑150 (MQ=255)
                                 gCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGaa                                                                                                                  >  1:1113524/1‑150 (MQ=255)
                                  cGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAAt                                                                                                                 <  2:907508/150‑1 (MQ=255)
                                       aGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCACCACATCGGCGCGCGGGTACCGGTGGCGTGGGTGGAATCCAAc                                                                                                            >  1:382388/1‑150 (MQ=255)
                                       aGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAAc                                                                                                            <  2:115533/150‑1 (MQ=255)
                                           cggCGCAGCACCAGCCACGCCATGGTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGTCGGGATGTCGATGTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGaaa                                                                                                        <  2:344600/149‑1 (MQ=255)
                                                        gCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAGCACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATtct                                                                                           >  1:885206/1‑150 (MQ=255)
                                                        gCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAGCACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATtct                                                                                           >  1:889666/1‑150 (MQ=255)
                                                         ccACGCCATGGTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGAGCCATCAACACAGCGGCGCGCGGGTACCGGGGGCGTGGCTGGAATCCAACGAAACCTCCTGCATtctc                                                                                          <  2:1157347/150‑1 (MQ=255)
                                                          cACCCCATGTTTTGCGGGTGGATCGCGGGCTTGCGAAAATGGCGCGCGTCAGCGGGCGGGAAGTCGATTTCGCGTCCTTGCGCCATCACCACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTctct                                                                                         <  2:56188/150‑1 (MQ=255)
                                                                  gTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCCGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACCTCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGCAACCTCCTGCATTCTCTACACTTcc                                                                                 >  1:1097074/1‑150 (MQ=255)
                                                                           gtggATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGAcc                                                                        >  1:457499/1‑150 (MQ=255)
                                                                           gtggATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCAGTCTCTTCACTGCCGGCACGccc                                                                        >  1:87924/1‑150 (MQ=255)
                                                                              gATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGc                                                                     <  2:406289/150‑1 (MQ=255)
                                                                                   cgGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAAcc                                                                >  2:694485/1‑150 (MQ=255)
                                                                                       cTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTaaa                                                            >  2:659787/1‑150 (MQ=255)
                                                                                       cTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTaaa                                                            <  1:352855/150‑1 (MQ=255)
                                                                                        tGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAg                                                           <  1:983079/150‑1 (MQ=255)
                                                                                        tGGCGAAAATGGCCCGCGTCACCGAGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCCCGGCTACCGCTGTCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTACCGCCACGACCGGCAAACCTAAAg                                                           >  2:1346677/1‑150 (MQ=255)
                                                                                           cGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCg                                                        <  2:704819/150‑1 (MQ=255)
                                                                                                    cgcgcgTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAGCACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTg                                               <  2:889666/150‑1 (MQ=255)
                                                                                                    cgcgcgTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAGCACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTg                                               <  2:885206/150‑1 (MQ=255)
                                                                                                           cAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCCCGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTcgg                                        >  2:1219576/1‑150 (MQ=255)
                                                                                                                            aTTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTTCCGGTGGCGTGGCTGGAATCCAACGACACCTGCGGCAGTCTCTTCACTTCCGGCACGACCGGGACACCTAAGGGCGTGCAGCGTGACGTCGGCTGATAAGCGGTGgcgg                       >  1:856314/1‑149 (MQ=255)
                                                                                                                                cgcgTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGACACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGc                   <  2:1043183/150‑1 (MQ=255)
                                                                                                                                cgcgTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGc                   >  2:520469/1‑150 (MQ=255)
                                                                                                                                cgcgTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGGGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGc                   <  2:38252/150‑1 (MQ=255)
                                                                                                                                      gtTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTc             >  2:1076787/1‑150 (MQ=255)
                                                                                                                                         gcgcCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGAt          <  1:472591/150‑1 (MQ=255)
                                                                                                                                          cgcCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATg         <  1:520469/150‑1 (MQ=255)
                                                                                                                                            ccATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCACCGAAACCTCCTGCATTCTCTACACTTCCGGCACGCCCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGCCCTCGAGGGa       >  1:110001/1‑150 (MQ=255)
                                                                                                                                                 aaCACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCGCCTGCATTCTCTACACTTACGGCACGACCGGCAAACCTTACGGCGTGCAGCGTGACGACTGCGGCTCTGCGGTGGCGCTGGCGACCGCGATGGACCCCa  >  2:256775/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AAAATCATTCCCTATAAAAAATTGCTCGACGATGCGATAAGTCAGGCGCAGCACCAGCCACGCCATGTTTTGCTGGTGGATCGCGGGCTGGCGAAAATGGCGCGCGTCAGCGGGCGGGATGTCGATTTCGCGTCGTTGCGCCATCAACACATCGGCGCGCGGGTACCGGTGGCGTGGCTGGAATCCAACGAAACCTCCTGCATTCTCTACACTTCCGGCACGACCGGCAAACCTAAAGGCGTGCAGCGTGACGTCGGCGGATATGCGGTGGCGCTGGCGACCTCGATGGACACCA  >  NC_000913/353213‑353507

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