Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,852,657 T→C Y43C (TAC→TGC)  ilvB ← acetohydroxy acid synthase I subunit IlvB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,852,6570TC97.6% 145.6 / ‑3.7 42Y43C (TAC→TGC) ilvBacetohydroxy acid synthase I subunit IlvB
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (9/32);  total (9/33)
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.45e-01

CTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCCCTGCTGTTCCAGGAAATGAACGATAAAT  >  NC_000913/3852581‑3852734
                                                                            |                                                                             
cTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAAc                                                                           <  1:698516/81‑1 (MQ=255)
   aGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAAcagg                                                                        <  1:698517/81‑1 (MQ=255)
   aGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAAcag                                                                         <  2:698518/80‑1 (MQ=255)
      gATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAAcaggcag                                                                     >  1:698585/1‑81 (MQ=255)
        tAAAGCCCGCGCCCTGTTCATGACGGGCCAGGATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGa                                                                   <  2:698519/81‑1 (MQ=255)
        tAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGa                                                                   <  1:698520/81‑1 (MQ=255)
          aaGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGAt                                                                  <  1:698521/80‑1 (MQ=255)
          aaGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCCTAAGGCATCGCAAACAGGCAGGATa                                                                 <  1:698522/81‑1 (MQ=255)
           aGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAg                                                                <  2:698523/81‑1 (MQ=255)
              ccGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAAc                                                             <  2:698524/81‑1 (MQ=255)
                 cgcCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGc                                                          >  2:698587/1‑81 (MQ=255)
                  gcCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGcc                                                         <  2:698525/81‑1 (MQ=255)
                   cccTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGccc                                                        <  2:698527/81‑1 (MQ=255)
                   cccTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGccc                                                        <  2:698526/81‑1 (MQ=255)
                     cTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCgg                                                      <  1:698528/81‑1 (MQ=255)
                        ttCATGGCGCGCCAGAATATGACGGATTTGCGTGCTTTGGCTTAACGCGTCGTAAACGGGGAGGATCGAACCACCGGGAAt                                                   <  2:892674/81‑1 (MQ=255)
                        ttCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGa                                                     <  2:698530/79‑1 (MQ=255)
                        ttCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAAt                                                   <  2:698529/81‑1 (MQ=255)
                         tCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATg                                                  >  1:698594/1‑81 (MQ=255)
                          cATGACGGGCCAGAATATGGCGGAGTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATg                                                  <  1:698531/80‑1 (MQ=255)
                           aTGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATGc                                                 <  1:698532/80‑1 (MQ=255)
                            tGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATGCCt                                               >  2:698602/1‑81 (MQ=255)
                              aCGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATGCCTGt                                             <  1:698533/81‑1 (MQ=255)
                                        atatGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTa                                   <  2:698534/81‑1 (MQ=255)
                                         tatGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTaa                                  <  1:698535/81‑1 (MQ=255)
                                             gCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGcc                              >  1:698603/1‑81 (MQ=255)
                                             gCGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGc                               <  1:698536/80‑1 (MQ=255)
                                              cGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGccc                             <  1:698537/81‑1 (MQ=255)
                                              cGGATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGccc                             >  1:698606/1‑81 (MQ=255)
                                               ggATTTGCGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCCct                            <  1:698538/81‑1 (MQ=255)
                                                      cGTGCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCCCTGCTGTTc                     <  1:698539/81‑1 (MQ=255)
                                                         gCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCCCTGCTGTTCCAg                  >  2:698598/1‑81 (MQ=255)
                                                         gCTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCCCTGCTGTTCCAg                  >  2:698600/1‑81 (MQ=255)
                                                          cTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCCCTGCTGTTCCAgg                 <  2:698540/81‑1 (MQ=255)
                                                          cTTTGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCCCTGCTGTTCCAgg                 <  2:698541/81‑1 (MQ=255)
                                                             tGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCCCTGCTGTTCCAGGa                <  2:698542/79‑1 (MQ=255)
                                                             tGGCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCCCTGCTGTTCCAGGa                <  2:698543/79‑1 (MQ=255)
                                                              ggCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCCCTGCTGTTCCAGGAAAt             >  1:698590/1‑81 (MQ=255)
                                                               gCTTAAGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCCCTGCTGTTCCAGGAAAt             <  1:698544/80‑1 (MQ=255)
                                                                   aaGGCATCGCAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCCCTGCTGTTCCAGGAAATGAAc         <  1:698545/80‑1 (MQ=255)
                                                                       cATCGCAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCCCTGCTGTTCCAGGAAATGAACGATaa    <  1:698546/81‑1 (MQ=255)
                                                                         tCGCAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCCCTGCTGTTCCAGGAAATGAACGATAAAt  <  1:698547/81‑1 (MQ=255)
                                                                            |                                                                             
CTGAGCGATAAAGCCCGCGCCCTGTTCATGACGGGCCAGAATATGGCGGATTTGCGTGCTTTGGCTTAAGGCATCGTAAACAGGCAGGATAGAACCGCCCGGAATGCCTGTCACAATCTTAATGCCCTGCTGTTCCAGGAAATGAACGATAAAT  >  NC_000913/3852581‑3852734

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