Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,172,353 A→T 7.3% L174Q (CTG→CAG)  gatD ← galactitol‑1‑phosphate 5‑dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,172,3530AT 7.3% 147.0 / 4.7 56L174Q (CTG→CAG) gatDgalactitol‑1‑phosphate 5‑dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (25/26);  new base T (2/2);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-02

CTGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGAC  >  NC_000913/2172211‑2172492
                                                                                                                                              |                                                                                                                                           
ctGCATTTGCGGCTCGCTCATCTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCcagca                                                                                                                                             <  2:350123/143‑1 (MQ=255)
        gcggcgCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCattc                                                                                                                                     <  2:928314/143‑4 (MQ=255)
        gcggcgCGCTCATTTCGCGACTGTTAAATGGTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATg                                                                                                                                     <  1:981219/143‑1 (MQ=255)
           gcgcgcTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCTGCTGGCCattcagt                                                                                                                                  <  2:468885/143‑7 (MQ=255)
              cgcTCATTTCGCTACTGTTAAACGCTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCg                                                                                                                               <  1:384941/143‑1 (MQ=255)
                cTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGc                                                                                                                             >  2:226732/1‑143 (MQ=255)
                cTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTCGCTCCCAGCGCGACAGCGCACTGCATGGCCAGCAGCCCAACGGCTCCGGc                                                                                                                             >  2:239014/1‑143 (MQ=255)
                 tCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCa                                                                                                                            >  1:501390/1‑143 (MQ=255)
                    tttCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATCGCCAGCAGGCCAATGGTTCCGGCAccc                                                                                                                         >  2:643655/1‑142 (MQ=255)
                        gCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGCCACCcataa                                                                                                                     >  2:248886/1‑143 (MQ=255)
                         ccaCTGGTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACACATCGCTCCTAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAACGGATCCGGCACCactaat                                                                                                                    >  2:1263624/3‑143 (MQ=255)
                         cTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATCGCAAGCCGGCCAATGGTTCCGGCACAactaat                                                                                                                    >  2:804497/1‑143 (MQ=255)
                               ttAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAAc                                                                                                              >  1:836416/1‑143 (MQ=255)
                                     gTTTGCATCGCACCGAAAGATTTTGCCAGTGCAAGTTATCTGGTATTAATGTGGATCGCCGTCACACTCTTTGTGGTCATCGCGATAGTGCACTGAATGTCGAGCAGTCTAATGTTTCCGGCAGCAATAAAAATAACGttttt                                                                                                        >  2:153408/1‑143 (MQ=255)
                                      tttGCATGGCACCGAATGATTTTGCCAGTGCCAGTGTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTa                                                                                                       <  2:707977/143‑1 (MQ=255)
                                         gCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAAGAACGTTTTTATtc                                                                                                    >  1:1165444/1‑143 (MQ=255)
                                          cATCGCACCGAACGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATtct                                                                                                   <  1:220376/143‑1 (MQ=255)
                                          cATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATgct                                                                                                   >  1:334570/1‑143 (MQ=255)
                                           aTCGCACCGAAAGAGTTTGCCAGTGCCAGTTGTTCTCAACTAATGTCGACCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATcctc                                                                                                  <  2:222957/143‑1 (MQ=255)
                                            tCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTAGATCGCCGTCATGTTCTGTGGTCACACACTTACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTca                                                                                                 >  1:295471/1‑143 (MQ=255)
                                             cGCACCGATAGAGTGTGCCCGTGCCAGTCTTTCTGCACTAATGTCGATCGCCGTCCCACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTcac                                                                                                <  2:613430/143‑1 (MQ=255)
                                              gCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCCCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGCTCCGGCACCACTAATAATCACGTTTTTATACCcaca                                                                                               >  2:84844/1‑143 (MQ=255)
                                               cACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACaa                                                                                              <  1:431544/143‑1 (MQ=255)
                                               cACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACaa                                                                                              <  1:818806/143‑1 (MQ=255)
                                               cACCGAAAGATTCTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACaa                                                                                              <  2:475556/143‑1 (MQ=255)
                                                      aGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTgcg                                                                                       <  1:766865/143‑1 (MQ=255)
                                                         ttttGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTa                                                                                    >  1:447552/1‑143 (MQ=255)
                                                                agtgccagtTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGaaa                                                                             <  1:261164/143‑1 (MQ=255)
                                                                 gtgccaggTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGaaaa                                                                            <  1:338736/143‑1 (MQ=255)
                                                                      ggtggtgcTGAACTAGTGGTGTGCGCAGTCACACGCGTTGCTGCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCGATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATg                                                                       <  2:507297/136‑1 (MQ=255)
                                                                           tgcTGAACTAATGCCGATCGCCGTCACACTCTTTGCCCCCAGCGGGTCAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGAc                                                                  <  2:836416/141‑1 (MQ=255)
                                                                              tGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCaa                                                               >  2:279219/1‑143 (MQ=255)
                                                                               gAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCAAAACCTTGCGCTAAATGAAACGCATGAAGACCCAc                                                              >  1:494903/1‑143 (MQ=255)
                                                                                      tGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATCCCCACAACCTTCGCCTAAATGAAACGCAGGAAGCCCAACGGaaacg                                                       >  2:869461/1‑141 (MQ=255)
                                                                                       gTCGATCGTCGTCGGACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATCGTTCCGGCACCAATAATAATAACGTTTTTAGTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCg                                                      <  2:917684/143‑1 (MQ=255)
                                                                                         cGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGATAATAGGc                                                    >  2:396508/1‑143 (MQ=255)
                                                                                          gATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCCCCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTCACTGAACACCTGCAGACCAACGGTAATCCGccc                                                  >  1:146134/1‑141 (MQ=255)
                                                                                                    acacTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGc                                         <  1:696175/143‑1 (MQ=255)
                                                                                                    acacTCCTTGCCCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGc                                         <  2:800856/143‑1 (MQ=255)
                                                                                                     cacTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACGAACGGTAATCGGCTCAATAAAAGcc                                        >  1:1129467/1‑143 (MQ=255)
                                                                                                      cctctTTGCTCCCAGCGCGACACCACACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGccc                                       <  2:501390/142‑1 (MQ=255)
                                                                                                      acTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGccc                                       >  2:749738/1‑143 (MQ=255)
                                                                                                       ctctTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGcccc                                      >  1:16512/1‑143 (MQ=255)
                                                                                                          tttGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCa                                                                     <  2:675617/109‑1 (MQ=255)
                                                                                                          tttGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCa                                                                     >  1:675617/1‑109 (MQ=255)
                                                                                                            tGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATACAAGCCCCATCCt                                 >  1:1182296/1‑143 (MQ=255)
                                                                                                                 ccggcgcGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATa                            <  2:295471/140‑1 (MQ=255)
                                                                                                                    gcgcGACAGCGCACTGAATGGCCAGCAGGCCAGTGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGc                         <  2:885582/143‑1 (MQ=255)
                                                                                                                            gcgcACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGt                 <  1:1143989/143‑1 (MQ=255)
                                                                                                                             ccaaCTGGATGCCCAGCATGTCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCATGCGCTAAATCAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCGCAATAGGCATATCCGTg                <  2:1109004/140‑1 (MQ=255)
                                                                                                                              gcACTGGATGGCCCGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTgg               <  2:164562/143‑1 (MQ=255)
                                                                                                                              gcACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTgg               <  1:298486/143‑1 (MQ=255)
                                                                                                                               cACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATATTAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGGCACATCCTCAATAGGCATATCCGTggg              >  1:105894/1‑143 (MQ=255)
                                                                                                                                   cattGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATGCGGGGGTaga          >  1:792492/4‑143 (MQ=255)
                                                                                                                                    attGGCCTGCTGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTagag         >  2:725312/3‑143 (MQ=255)
                                                                                                                                      tGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCa       <  2:930813/143‑1 (MQ=255)
                                                                                                                                         ccAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCACCCTCAATAGGCATATCCGCGCGTAGAGCAAAg    >  2:937599/1‑143 (MQ=255)
                                                                                                                                          cagcagGCCACTGGGTCCCGCACCAATAATAATAACGCTTTTACTCTCACAGCCTTGCGCTAAATGAAAAGCATGCAGCCCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGa   <  2:937236/143‑1 (MQ=255)
                                                                                                                                          cagcagGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATACGCATATCCCTCGGCAGAGCAAAGa   >  2:256562/1‑143 (MQ=255)
                                                                                                                                           agcagGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGCGGGTAGAGCAAAGAc  >  1:150358/1‑143 (MQ=255)
                                                                                                                                              |                                                                                                                                           
CTGCATTTGCGGCGCGCTCATTTCGCTACTGTTAAATGTTTGCATCGCACCGAAAGATTTTGCCAGTGCCAGTTTTTCTGAACTAATGTCGATCGCCGTCACACTCTTTGCTCCCAGCGCGACAGCGCACTGAATGGCCAGCAGGCCAATGGTTCCGGCACCAATAATAATAACGTTTTTATTCTCACAACCTTGCGCTAAATGAAAAGCATGCAGACCAACGGTAATCGGCTCAATAAAAGCCCCATCCTCAATAGGCATATCCGTGGGTAGAGCAAAGAC  >  NC_000913/2172211‑2172492

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