Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,953,437 A→G 8.3% intergenic (+24/‑41) ilvE → / → ilvD branched‑chain amino acid aminotransferase/dihydroxyacid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,953,4370AG 8.3% 148.4 / 2.8 48intergenic (+24/‑41)ilvE/ilvDbranched‑chain amino acid aminotransferase/dihydroxyacid dehydratase
Reads supporting (aligned to +/- strand):  ref base A (25/19);  new base G (2/2);  total (27/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

GGTTGGCGAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATT  >  NC_000913/3953299‑3953585
                                                                                                                                          |                                                                                                                                                    
ggTTGGCGAAGGCCGTTTTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGTGCGTGCCGTCCCAtt                                                                                                                                            <  3:229698/149‑1 (MQ=255)
       gAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACgaga                                                                                                                                     <  4:291119/149‑1 (MQ=255)
           gCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGagacaga                                                                                                                                 >  4:493407/1‑149 (MQ=255)
                tgtgGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGACGTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGGTCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGACACAGACACTg                                                                                                                            <  4:600112/149‑1 (MQ=255)
                     cccGGTTACCAAACGCATACAACAAGACTTCTTCGGCCTCTTAACTTGAGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGGTTCACAATACATACAATAAATGGAACAGCACGCACCGTCCCATTTACGAGACAGACACCGGGAGt                                                                                                                       >  2:58084/1‑149 (MQ=255)
                      ccGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGACACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAATAATGGGACGGCACCCACCGTCCCATTTACGAGACAGACACGGGGAgca                                                                                                                      >  1:646800/1‑147 (MQ=255)
                       cGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGtaa                                                                                                                     >  4:407182/1‑149 (MQ=255)
                       cGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGtaa                                                                                                                     >  3:320688/1‑149 (MQ=255)
                        ggTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGtaaa                                                                                                                    <  4:614778/149‑1 (MQ=255)
                               aaaCGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGTGCGTGCCGTCCCa                                                                                                                                              >  4:363706/1‑116 (MQ=255)
                               aaaCGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGAGCAAGTTAATCAATAAATACAAAAAATGGGACGGTGCGTGCCGTCCCa                                                                                                                                              <  3:363706/116‑1 (MQ=255)
                                aaCGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTAt                                                                                                            >  1:208488/1‑149 (MQ=255)
                                 aCGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATg                                                                                                           >  1:140335/1‑149 (MQ=255)
                                          gCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTAc                                                                                                  >  2:108089/1‑149 (MQ=255)
                                             aGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGagacaga                                                                                                                                 >  1:497540/1‑115 (MQ=255)
                                             aGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGagacaga                                                                                                                                 <  2:497540/115‑1 (MQ=255)
                                               ccTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTgg                                                                                                                           <  1:213081/119‑1 (MQ=255)
                                               ccTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTgg                                                                                                                           >  2:213081/1‑119 (MQ=255)
                                                  tcttcGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGc                                                                                          <  1:108089/149‑1 (MQ=255)
                                                        gCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGTGCGTGCCGTcc                                                                                                                                                >  4:528179/1‑89 (MQ=39)
                                                        gCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTcc                                                                                            >  2:153659/1‑141 (MQ=255)
                                                        gCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTcc                                                                                            <  1:153659/141‑1 (MQ=255)
                                                           tctTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTAc                                                                                                  <  3:226257/132‑1 (MQ=255)
                                                           tctTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTAc                                                                                                  >  4:226257/1‑132 (MQ=255)
                                                                 cTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATgg                                                                           <  3:564544/149‑1 (MQ=255)
                                                                  tGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAg                                                                                                     <  4:463776/122‑1 (MQ=255)
                                                                  tGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAg                                                                                                     >  3:463776/1‑122 (MQ=255)
                                                                   ggCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTc                                                                         <  3:493407/149‑1 (MQ=255)
                                                                            gAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGGCGTAATATgg                                                                >  2:181543/1‑149 (MQ=255)
                                                                                     tGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTc                                                       <  3:407182/149‑1 (MQ=255)
                                                                                      ggggCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCg                                                      >  1:528324/1‑149 (MQ=255)
                                                                                         gCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTc                                                                         >  2:377372/1‑127 (MQ=255)
                                                                                         gCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTc                                                                         <  1:377372/127‑1 (MQ=255)
                                                                                             gTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCATAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCtg                                               >  4:593844/1‑149 (MQ=255)
                                                                                              ttAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCATCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCtgt                                              <  4:115890/149‑1 (MQ=255)
                                                                                                  aTCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGT‑CCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGgcgc                                         >  2:483323/1‑149 (MQ=255)
                                                                                                        ttAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTACGCCACCACCACTCATGGTCGTAATATGGCGGGTTCTCGAGCGCTGTGGCGCGCCAc                                    >  4:674706/1‑149 (MQ=255)
                                                                                                            tCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGa                                >  2:614775/1‑149 (MQ=255)
                                                                                                                    taCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATg                                                                            >  2:101041/1‑97 (MQ=255)
                                                                                                                    taCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATg                                                                            <  1:101041/97‑1 (MQ=255)
                                                                                                                     aCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATg                                                                            <  2:138611/96‑1 (MQ=255)
                                                                                                                     aCAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATg                                                                            >  1:138611/1‑96 (MQ=255)
                                                                                                                       aaaaaaTGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTAt                                                                                                            <  2:338644/62‑1 (MQ=255)
                                                                                                                       aaaaaaTGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTAt                                                                                                            >  1:338644/1‑62 (MQ=255)
                                                                                                                           aaTGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTCCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGAt                 <  4:320688/149‑1 (MQ=255)
                                                                                                                              gggACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACCCCGATTTc              >  3:56320/1‑149 (MQ=255)
                                                                                                                               ggACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCg             <  2:528324/149‑1 (MQ=255)
                                                                                                                                gACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCgg            >  4:460129/1‑149 (MQ=255)
                                                                                                                                        gcacCGTTCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGa    <  3:593844/149‑1 (MQ=255)
                                                                                                                                         cacCGTCCCATTTACGGGACAGACACTGGGAGTAAACAAAGTTTGCCTAAGCACCGTTGCGCCACCACCACTCATGGTCGTAATACGGCGGGTGCTTGTGCGCTGTGTCGCGCCAACGGAATGAGCGACTCTGATTTCGGTAAGcggtt   <  2:646800/149‑5 (MQ=255)
                                                                                                                                         cacCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGAt   <  4:396788/149‑1 (MQ=255)
                                                                                                                                          acCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGTGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGAtt  <  2:140335/149‑1 (MQ=255)
                                                                                                                                          acCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGAtt  >  1:99487/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                    
GGTTGGCGAAGGCCGTTGTGGCCCGGTTACCAAACGCATTCAGCAAGCCTTCTTCGGCCTCTTCACTGGCGAAACCGAAGATAAATGGGGCTGGTTAGATCAAGTTAATCAATAAATACAAAAAATGGGACGGCACGCACCGTCCCATTTACGAGACAGACACTGGGAGTAAATAAAGTATGCCTAAGTACCGTTCCGCCACCACCACTCATGGTCGTAATATGGCGGGTGCTCGTGCGCTGTGGCGCGCCACCGGAATGACCGACGCCGATTTCGGTAAGCCGATT  >  NC_000913/3953299‑3953585

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑