Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 87,417 G→T G21C (GGC→TGC)  ilvH → acetolactate synthase/acetohydroxybutanoate synthase, regulatory subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091387,4170GT98.9% 321.5 / ‑6.5 89G21C (GGC→TGC) ilvHacetolactate synthase/acetohydroxybutanoate synthase, regulatory subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base T (49/39);  minor base C (1/0);  total (50/39)
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.96e-01

TATCAGTCTTACTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTGGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCTGACCGTTGC  >  NC_000913/87370‑87463
                                               |                                              
tatCAGTCTTACTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGc                                              >  1:2950239/1‑50 (MQ=255)
 atCAGTCTTACTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGcc                                             >  1:6069280/1‑50 (MQ=255)
 atCAGTCTTACTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGcc                                             >  1:1027970/1‑50 (MQ=255)
 atCAGTCTTACTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGcc                                             >  1:90219/1‑50 (MQ=255)
 atCAGTCTTACTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGcc                                             >  1:3775391/1‑50 (MQ=255)
  tCAGTCTTACTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGctt                                            >  1:1854099/1‑48 (MQ=255)
  tCAGTCTTACTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGcc                                             <  1:1503509/49‑1 (MQ=255)
  tCAGTCTTACTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCt                                            >  1:2185707/1‑50 (MQ=255)
   cAGTCTTACTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCtt                                           >  1:4279587/1‑50 (MQ=255)
   cAGTCTTACTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCtt                                           >  1:4625129/1‑50 (MQ=255)
    aGTCTTACTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCtt                                           <  1:1851319/49‑1 (MQ=255)
     gTCTTACTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCtttt                                         <  1:2424492/50‑1 (MQ=255)
     gTCTTACTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCtttt                                         <  1:5653720/50‑1 (MQ=255)
     gTCTTACTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCtttt                                         >  1:2025047/1‑50 (MQ=255)
      tCTTACTAGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCttttt                                        <  1:4080653/50‑1 (MQ=38)
       cTTACTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCtttttt                                       >  1:5743663/1‑50 (MQ=255)
       cTTACTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCtttttt                                       >  1:5992227/1‑50 (MQ=255)
        ttACTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTc                                      >  1:2032668/1‑50 (MQ=255)
        ttACTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTc                                      >  1:5246765/1‑50 (MQ=255)
        ttACTCGAAAATGAATCAGGCGCGTAATCCCGCGAGATTCGCCTTTTTTc                                      >  1:3729790/1‑50 (MQ=255)
         tACTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTcc                                     >  1:2790441/1‑50 (MQ=255)
          aCTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCtttttt                                       <  1:2822062/47‑1 (MQ=255)
          aCTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCtttttt                                       <  1:2917654/47‑1 (MQ=255)
          aCTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTc                                      <  1:5082384/48‑1 (MQ=255)
          aCTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTc                                      <  1:3763109/48‑1 (MQ=255)
           cTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTccc                                    <  1:2228059/49‑1 (MQ=255)
           cTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTcc                                     >  1:2566128/1‑48 (MQ=255)
           cTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCa                                   <  1:1978398/50‑1 (MQ=255)
              gAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAg                                  <  1:3027112/48‑1 (MQ=255)
              gAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGc                                 <  1:2764734/49‑1 (MQ=255)
              gAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGc                                 <  1:5775891/49‑1 (MQ=255)
              gAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGc                                 <  1:4914632/49‑1 (MQ=255)
              gAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCg                                >  1:6087039/1‑50 (MQ=255)
              gAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCg                                >  1:885854/1‑50 (MQ=255)
              gAAAATGAATCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCg                                >  1:6573361/1‑50 (MQ=255)
              gAAAATGAATCAGGCGCGTTATCCCGCGCGATTTGCCTTTTTTCCCAGCg                                >  1:6707692/1‑50 (MQ=255)
                aaaTGAATCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTg                              >  1:279376/1‑50 (MQ=255)
                aaaTGAATCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTg                              >  1:1116408/1‑50 (MQ=255)
                 aaTGAATCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTgg                             >  1:2450040/1‑50 (MQ=255)
                  aTGAATCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGc                            >  1:6566069/1‑50 (MQ=255)
                   tGAATCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCt                           <  1:6293483/50‑1 (MQ=255)
                    gAATCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTa                          >  1:2116183/1‑50 (MQ=255)
                     aaTCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTac                         <  1:4230064/50‑1 (MQ=255)
                     aaTCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTa                          <  1:5377794/49‑1 (MQ=255)
                     aaTCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTa                          >  1:929055/1‑49 (MQ=255)
                     aaTCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTa                          <  1:2158732/49‑1 (MQ=255)
                     aaTCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTa                          >  1:3055177/1‑49 (MQ=255)
                       tCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTacaa                       <  1:3967095/50‑1 (MQ=255)
                       tCAGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTacaa                       >  1:588782/1‑50 (MQ=255)
                         aGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTacaaca                     >  1:5896022/1‑50 (MQ=255)
                         aGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTacaac                      <  1:1287231/49‑1 (MQ=255)
                         aGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTacaac                      <  1:612302/49‑1 (MQ=255)
                         aGGCGCGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTacaac                      <  1:2268718/49‑1 (MQ=255)
                           gcgcGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACAtt                   <  1:3680532/50‑1 (MQ=255)
                           gcgcGTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACAtt                   >  1:5883924/1‑50 (MQ=255)
                            cgcgTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTg                  <  1:2015016/50‑1 (MQ=255)
                             gcgTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTg                  >  1:6455298/1‑49 (MQ=255)
                              cgTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGaa                <  1:1404263/50‑1 (MQ=255)
                              cgTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGaa                <  1:2208626/50‑1 (MQ=255)
                               gTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGa                 >  1:2802590/1‑48 (MQ=255)
                               gTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGa                 >  1:5988389/1‑48 (MQ=255)
                               gTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGa                 >  1:6571096/1‑48 (MQ=255)
                               gTTATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGa                 >  1:1426021/1‑48 (MQ=255)
                                ttATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAg              <  1:1278371/50‑1 (MQ=255)
                                 tATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGc             >  1:143318/1‑50 (MQ=255)
                                 tATCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGc             >  1:3370038/1‑50 (MQ=255)
                                  aTCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGcc            >  1:4896084/1‑50 (MQ=255)
                                   tCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGcc            <  1:2845905/49‑1 (MQ=255)
                                   tCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCt           <  1:191182/50‑1 (MQ=255)
                                   tCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCt           <  1:6135371/50‑1 (MQ=255)
                                   tCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCt           <  1:6743315/50‑1 (MQ=255)
                                   tCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCt           <  1:341521/50‑1 (MQ=255)
                                   tCCCGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCt           >  1:2140387/1‑50 (MQ=255)
                                    cccGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCt           >  1:3602043/1‑49 (MQ=255)
                                     ccGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCTGa         <  1:2450437/50‑1 (MQ=255)
                                     ccGCGTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCTGa         >  1:314278/1‑50 (MQ=255)
                                       gcgTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCTGAcc       >  1:5752332/1‑50 (MQ=255)
                                       gcgTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCTGAcc       <  1:4439167/50‑1 (MQ=255)
                                       gcgTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCTGAcc       <  1:3288127/50‑1 (MQ=255)
                                       gcgTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCTGAcc       >  1:3797442/1‑50 (MQ=255)
                                        cgTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCTGACCg      >  1:2171537/1‑50 (MQ=255)
                                        cgTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCTGACCg      <  1:4794249/50‑1 (MQ=255)
                                        cgTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCTGACCg      >  1:6745156/1‑50 (MQ=255)
                                        cgTGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCTGACCg      >  1:1165474/1‑50 (MQ=255)
                                          tGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCTGACCGtt    >  1:2765838/1‑50 (MQ=255)
                                          tGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCTGACCGtt    >  1:3182370/1‑50 (MQ=255)
                                          tGATTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCTGACCGtt    >  1:716829/1‑50 (MQ=255)
                                            aTTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCTGACCGTTgc  <  1:6041092/50‑1 (MQ=255)
                                            aTTTGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCTGACCGTTg   <  1:564490/49‑1 (MQ=255)
                                               |                                              
TATCAGTCTTACTCGAAAATGAATCAGGCGCGTTATCCCGCGTGATTGGCCTTTTTTCCCAGCGTGGCTACAACATTGAAAGCCTGACCGTTGC  >  NC_000913/87370‑87463

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