Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 125,387 C→T Q791* (CAG→TAG)  aceE → pyruvate dehydrogenase E1 component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913125,3870CT98.0% 139.7 / NA 51Q791* (CAG→TAG) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (19/30);  minor base G (1/0);  total (20/31)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAA  >  NC_000913/125247‑125526
                                                                                                                                            |                                                                                                                                           
tcaTCCTTGCGCTAGTTTTCGGCTTAGGTTCTTTCGTTTATATCTTTACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAg                                                                                                                                           <  1:119432/141‑1 (MQ=255)
     cTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTGAGGtgag                                                                                                                                      >  1:805668/1‑142 (MQ=255)
       ggCGAAAGATGATGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCAACGAGCTGCCGCTTGATGGTCAGGATTGTGAACGCTGGAACATTCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGtgatga                                                                                                                                    <  2:1098886/143‑1 (MQ=255)
           aaaGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacga                                                                                                                                <  1:670206/143‑1 (MQ=255)
            aaGATTACGTCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacgac                                                                                                                               >  1:807070/1‑143 (MQ=255)
            aaGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacgac                                                                                                                               >  1:848885/1‑143 (MQ=255)
             aGATTGCGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacgacg                                                                                                                              <  2:498143/143‑1 (MQ=255)
             aGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAacgacg                                                                                                                              <  2:494573/143‑1 (MQ=255)
                     gCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCa                                                                                                                      >  1:702687/1‑143 (MQ=255)
                       gTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGt                                                                                                                    <  1:19914/143‑1 (MQ=255)
                                   gTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGGGAGGTGCGGTGGTCCCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGa                                                                                                        >  1:1256443/1‑143 (MQ=255)
                                    tttATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGAc                                                                                                       >  2:1093695/1‑143 (MQ=255)
                                    tttATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGAc                                                                                                       <  1:426660/143‑1 (MQ=255)
                                    tttATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGCGGCATCTACCGAc                                                                                                       >  1:244176/1‑143 (MQ=255)
                                     ttATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGCCAGTGGCACCTACCggcg                                                                                                      >  1:297628/1‑140 (MQ=255)
                                       ataGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGGATATCGCTTCGGTGCTGAACGACGCTCCGGCAGTGGCATCTACCGACtct                                                                                                    >  2:473228/1‑141 (MQ=255)
                                        taGCGTGACCTCCTTCAACGAGCTGGCGCGGGATGGGCAGGAGTGGGCGCGCTGGCCCATGCTGCACCCGGTGGAAACTCCGCGCGTGCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCAGCTACCGACtata                                                                                                   <  1:243784/143‑1 (MQ=255)
                                           cGTGACGTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAGCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGa                                                                                                <  2:173261/143‑1 (MQ=255)
                                                  tcctTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGtt                                                                                         <  2:250231/143‑1 (MQ=255)
                                                              cTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGGATATCGCGTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGt                                                                             >  1:1039153/1‑143 (MQ=255)
                                                                 gcgcGTGATGGTCAGGATTGTGATCGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCGTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCg                                                                          <  2:1128198/143‑1 (MQ=255)
                                                                 gcgcGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACGGTATGAGACTGCTCGCTGAGCAGGGCCg                                                                          >  2:933/1‑143 (MQ=255)
                                                                   gcgTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTACgac                                                                        >  1:279928/1‑141 (MQ=255)
                                                                   gcgTGATGGCCAGGATCGTGACCGCCGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATAGCGCTTAGGTGATGACCGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTa                                                                        <  1:1213088/143‑1 (MQ=255)
                                                                       gATGGTCAGGAGTGTGAACGCTGGAACAGGCTGCGCCCGCTGGAAACTCCGCGCGGTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTa                                                                    <  2:319909/143‑1 (MQ=255)
                                                                         tGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACg                                                                  <  1:1236689/143‑1 (MQ=255)
                                                                            tCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAGCGACGCGCCGGCAGTGGCATCTACCGACTGTGTGGAACGATCCGCTGAGCGGGTCCGTTCTCACGTAc                                                               >  2:311149/1‑143 (MQ=255)
                                                                              aGGATTGTGACGGCTGGAACATGCTGCACGCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCg                                                             <  2:124788/143‑1 (MQ=255)
                                                                              aGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCg                                                             <  1:915437/143‑1 (MQ=255)
                                                                                 aTTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTGCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGCGGCATCTACCGACTATATGAAACTGTCCGCTGAGCAGCTCCCTACTTACGTACCgaat                                                          >  2:877429/1‑140 (MQ=255)
                                                                                    gtgAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGCCGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTTCCGTCTtac                                                       >  1:259376/1‑143 (MQ=255)
                                                                                                   tgctgcCCCCGCTGGAAACTCCGCTCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAACCTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTg                                        <  1:51835/143‑1 (MQ=255)
                                                                                                    gctgcACCCGCTGGAACCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGCCAGTGTCTTCTACTGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTgg                                       <  1:990081/143‑1 (MQ=255)
                                                                                                    gctgcACCCGCTGGAACCCCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAACCTGTGCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTgg                                       <  1:740746/143‑1 (MQ=255)
                                                                                                     ctgcACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGCCTACCGCGTACTggg                                      >  1:765711/1‑143 (MQ=255)
                                                                                                            cGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGAt                               <  1:1144555/143‑1 (MQ=255)
                                                                                                            cGCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGAt                               <  2:864417/143‑1 (MQ=255)
                                                                                                             gCTGGAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTCCTGGGTACTGATg                              <  2:702687/143‑1 (MQ=255)
                                                                                                              cgggACACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGGACCGGCTGACGACTACCGCGTACTGGGTACTGATgg                             <  1:819860/141‑1 (MQ=255)
                                                                                                                ggAAACTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTGTATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCt                           >  1:475086/1‑143 (MQ=255)
                                                                                                                   aaCTCCGCGCTTTCCTTATATCGCTTAGGTGATGCACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCg                        <  2:119645/143‑1 (MQ=255)
                                                                                                                   aaCTCCGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCg                        <  1:1093695/143‑1 (MQ=255)
                                                                                                                       ccGCGCGTTCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCg                    >  2:1259609/1‑143 (MQ=255)
                                                                                                                         gcgcgTTCGGTAGCTGCCTCAGGGGATGATCGACGCTCCGGCAGTCGCAGCTACCGACTATATGACACGGGCCGCTGAGCAGGTCCGTACTTACGTACCGGCTGCCGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGtt                  <  2:965135/143‑1 (MQ=255)
                                                                                                                         gcgcgTTCCGTATATCTCTTAGCGGATGAACGACGCTCCCGCAGTGGCATCTACCGACTATATGAAACCGTTCGCTGAGCTGGTCCGCACCTACGTACCGGCTGACGACGTCCCCGTACCCGGTACACATAGCTGCGGCCgct                  >  2:761224/1‑141 (MQ=255)
                                                                                                                             gTTCCGTATAGCTCTTAGGTTATGAACGGCGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGa              <  2:48419/143‑1 (MQ=255)
                                                                                                                              ttCCGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGAGGGCTTCGGTCGTTCCGAc             >  1:304085/1‑143 (MQ=255)
                                                                                                                               tCCGTATATCGCTTAGTTGTTGAATGACGTTCCGGCTGTGGCATCTAGCGACTATATGAAGCTGTTCGCTGAGGAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACa            <  1:538464/143‑1 (MQ=255)
                                                                                                                                ccGTATATCGCTTAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAg           <  1:1140015/143‑1 (MQ=255)
                                                                                                                                    atatCGCCTAGGTCATGTACGACGCTCCGGTGGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGGACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGt       <  2:289927/143‑1 (MQ=255)
                                                                                                                                         gCTTAGTTGATGAACGACGCTCCGGCAGTGGCAGCATGCGACTATATAAAACTTTGCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGaa  <  1:346916/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                           
AGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCTGGAAACTCCGCGCGTTCCGTATATCGCTCAGGTGATGAACGACGCTCCGGCAGTGGCATCTACCGACTATATGAAACTGTTCGCTGAGCAGGTCCGTACTTACGTACCGGCTGACGACTACCGCGTACTGGGTACTGATGGCTTCGGTCGTTCCGACAGCCGTGAGAA  >  NC_000913/125247‑125526

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: