Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,954,274 A→C 7.2% E266A (GAA→GCA) ‡ ilvD → dihydroxy‑acid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,954,2740AC 7.2% 128.1 / 2.4 56E266A (GAA→GCA) ‡ilvDdihydroxy‑acid dehydratase
Reads supporting (aligned to +/- strand):  ref base A (25/27);  new base C (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 = 1.80e-01

GCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTAAA  >  NC_000913/3954133‑3954416
                                                                                                                                              |                                                                                                                                              
gCTCGCTGCTGGCACCCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGaa                                                                                                                                               <  2:568647/143‑1 (MQ=255)
gCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGaa                                                                                                                                               <  1:219848/143‑1 (MQ=255)
     ctgctgGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGc                                                                                                                                          <  1:360651/143‑1 (MQ=255)
      tgctgGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTTCAAACGcc                                                                                                                                         <  1:607505/143‑1 (MQ=255)
        cttgCAACTTACGCCGACCGTAACTAGCTGTTCCCTAATGGTGGTAAACGCATTGTCGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTTCAAACGccgc                                                                                                                                       <  2:174714/140‑3 (MQ=255)
        ctgGCAAGCCACGACGACCGTAAGCAGCTGTTCCTTAATGCTGGGAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCAt                                                                                                                                       <  1:563656/143‑1 (MQ=255)
        ctgGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTGAGGCGGCGTTTGAAAACGCCag                                                                                                                                       >  2:613529/1‑142 (MQ=255)
         tgGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAGCGTTATTACGAGCAAAACGACGAAAGTGCACTGACGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCagg                                                                                                                                      >  1:372090/1‑141 (MQ=255)
          ggCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTTAAAACGCCATGa                                                                                                                                     >  1:356680/1‑143 (MQ=255)
          ggCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGa                                                                                                                                     <  1:107308/143‑1 (MQ=255)
              aCCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGGATGGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGTGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCt                                                                                                                                 <  1:54865/143‑1 (MQ=255)
                   cGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTCAAAACGCCATGACGCTGCata                                                                                                                            >  2:190263/1‑143 (MQ=255)
                       gACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGCTATTACGAGCAAAACGACGAAAGTGCAGTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATcgc                                                                                                                        <  1:357245/143‑1 (MQ=255)
                         ccGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGa                                                                                                                      <  1:366596/143‑1 (MQ=255)
                          cGTAAGCAGCTGTTCATTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAACCGCCATGACGCTGGATATCGCGAt                                                                                                                     >  2:526386/1‑143 (MQ=255)
                            tAAGGAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATgg                                                                                                                   <  2:348436/143‑1 (MQ=255)
                            tAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATgg                                                                                                                   <  2:296295/143‑1 (MQ=255)
                                         cTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTCAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTaa                                                                                                      >  1:327964/1‑143 (MQ=255)
                                          ttAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGAGGAACGACTaaa                                                                                                     >  1:324897/1‑142 (MQ=255)
                                          ttAATGCTGGTAAACGCATTGTAGAAATGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAac                                                                                                     >  2:546615/1‑143 (MQ=255)
                                             gtGCTGGTAAACGCATTGTTTAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGTCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACAcc                                                                                                  <  2:589923/142‑1 (MQ=255)
                                                cTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTa                                                                                               <  1:344023/143‑1 (MQ=255)
                                                 tGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTAc                                                                                              >  1:306013/1‑143 (MQ=255)
                                                  ggTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACt                                                                                             >  1:332372/1‑143 (MQ=255)
                                                   gTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACtt                                                                                            >  1:140039/1‑143 (MQ=255)
                                                    tAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCCAGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTCCTTc                                                                                           >  1:139769/1‑143 (MQ=255)
                                                    caaaCGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTc                                                                                           <  1:330642/142‑1 (MQ=255)
                                                            ttgttgAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACctgctg                                                                                   >  1:110084/1‑143 (MQ=255)
                                                               ttgAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTggcc                                                                                >  2:197675/1‑142 (MQ=255)
                                                                 gAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTggcggc                                                                              >  1:569302/1‑143 (MQ=255)
                                                                      gACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGgcgc                                                                         <  1:360671/143‑1 (MQ=255)
                                                                       aCCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCa                                                                        >  1:386739/1‑143 (MQ=255)
                                                                       aCCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGAGTAAAAGCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGCTCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCa                                                                        >  1:171323/1‑143 (MQ=255)
                                                                        ccgaaCGTTTTTACGAGTAAAACGATGAAAGTGCACTGTCGCGTAATATCGCCTGTAAGGCGGCGGTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAg                                                                       <  1:396206/140‑1 (MQ=255)
                                                                               cttataCGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCgg                                                                <  2:162095/139‑1 (MQ=255)
                                                                                      caGCAAAACGACGAAGGGGCACTGCCGCGTAATACCGCCAGGAAGGCGGCGTTGGAAGACGCCATGACGGTGGCTATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGa                                                         <  1:80007/142‑1 (MQ=255)
                                                                                        gCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAAAACCGTACTTCGCCTGCTGGCGGCGGGGCAGGAGGAGGAAATCGACt                                                       >  1:48362/1‑143 (MQ=255)
                                                                                           aaaCGACGAAAGTGCACTGCCGCGTGATATCGCCAGTAAGGCGGCGGTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCa                                                    >  2:218201/1‑143 (MQ=255)
                                                                                                acgaAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATg                                               >  1:568579/1‑143 (MQ=255)
                                                                                                 cgaAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGa                                              <  2:388792/143‑1 (MQ=255)
                                                                                                 cgaAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGa                                              <  2:246313/143‑1 (MQ=255)
                                                                                                   aaaggacACTGCCGCGTGATAGCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGt                                            <  1:528750/137‑1 (MQ=255)
                                                                                                     aGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCAAGAGTgg                                          >  1:110150/1‑142 (MQ=255)
                                                                                                        gCACTGCCGCGTAATCTACACCAGATGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTCACACCGTACTTCACCTGCTGGCGGCGGCGCAGGCAGCGGAAAGCGACTTCAGCATGAGAGatag                                       >  2:379532/1‑142 (MQ=255)
                                                                                                          aCTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGAGGCTGGATATCGCGATGGGTGGCTCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCCGGAAGCGGAAAGCGACGTCCCCCTGCGTGTTATCg                                     >  1:434606/1‑143 (MQ=255)
                                                                                                          aCTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCg                                     >  1:197348/1‑143 (MQ=255)
                                                                                                            tGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGAt                                   <  2:110084/143‑1 (MQ=255)
                                                                                                            tGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGAt                                   >  1:473505/1‑143 (MQ=255)
                                                                                                            gtcctcgTACTATCGCGCGTAAGGCTGCGGTTGAAGACGCCATGAGGCTGGACAGCGCGATGGGTGGATCGACTAACACCGTATTTCACCTCCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGAt                                   <  2:241046/138‑1 (MQ=255)
                                                                                                              ccGCGTAATATCGCCAGTAGGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATaa                                 <  2:491815/143‑1 (MQ=255)
                                                                                                              ccGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATaa                                 <  1:403112/143‑1 (MQ=255)
                                                                                                                gcgTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGc                               <  1:565586/143‑1 (MQ=255)
                                                                                                                            ccAGCAAGGCGGCGTTGGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAg                    <  2:233847/143‑1 (MQ=255)
                                                                                                                               gcgtcatggcgTTTTCAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGtt                 >  2:292429/8‑143 (MQ=255)
                                                                                                                               gTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGtt                 <  1:515487/143‑1 (MQ=255)
                                                                                                                                 aaGGCGGCGTTTGATAGCGCCAGGACGCAGGATATCGCGATGCGTGGATTCACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTcc               <  2:69300/143‑1 (MQ=255)
                                                                                                                                 aaGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTcc               <  1:307799/143‑1 (MQ=255)
                                                                                                                                    atggcgTTTTCAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCCCCATGCGAGATATCGATAAGCTCTCCCCCACGCTTCcaca            >  2:370787/3‑143 (MQ=255)
                                                                                                                                              aaaaCGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGATGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTaaa  <  2:306013/143‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                              
GCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACAGCTGTGTAAA  >  NC_000913/3954133‑3954416

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: