Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,954,273 G→T 8.4% E266* (GAA→TAA) ‡ ilvD → dihydroxy‑acid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,954,2730GT 8.4% 139.2 / 3.5 58E266* (GAA→TAA) ‡ilvDdihydroxy‑acid dehydratase
Reads supporting (aligned to +/- strand):  ref base G (23/29);  new base T (3/2);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

GCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACA  >  NC_000913/3954133‑3954406
                                                                                                                                             |                                                                                                                                     
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)
                                                                                                                                             |                                                                                                                                     
GCTCGCTGCTGGCAACCCACGCCGACCGTAAGCAGCTGTTCCTTAATGCTGGTAAACGCATTGTTGAATTGACCAAACGTTATTACGAGCAAAACGACGAAAGTGCACTGCCGCGTAATATCGCCAGTAAGGCGGCGTTTGAAAACGCCATGACGCTGGATATCGCGATGGGTGGATCGACTAACACCGTACTTCACCTGCTGGCGGCGGCGCAGGAAGCGGAAATCGACTTCACCATGAGTGATATCGATAAGCTTTCCCGCAAGGTTCCACA  >  NC_000913/3954133‑3954406

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: