Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,217,424 T→C Y106H (TAT→CAT)  aceA → isocitrate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,217,4240TC97.6% 133.6 / ‑4.5 41Y106H (TAT→CAT) aceAisocitrate lyase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (17/23);  total (18/23)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

AGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACT  >  NC_000913/4217287‑4217559
                                                                                                                                         |                                                                                                                                       
aGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGGGGCGGACGGCAACCTGGCGGCCAgcatgcat                                                                                                                                       >  2:367210/1‑140 (MQ=255)
  tCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCGGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATcc                                                                                                                                     <  2:774438/140‑1 (MQ=255)
     aaaaaaGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGa                                                                                                                                  >  1:48115/1‑140 (MQ=255)
              tACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGct                                                                                                                         <  1:829343/140‑1 (MQ=255)
                cATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGctct                                                                                                                       >  2:509209/1‑140 (MQ=255)
                cATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCCAACCTGGCGGCCAGCATGCATCCGGGTCAGTCGctct                                                                                                                       >  2:853102/1‑140 (MQ=255)
                cATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTACCGGCGGACTCTAAACTAGCGGCCCGGATGCGTCCGGATCAGGCGctct                                                                                                                       >  2:284556/1‑140 (MQ=255)
                            cGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACt                                                                                                           <  2:22784/140‑1 (MQ=255)
                             ggCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTc                                                                                                          <  1:617644/140‑1 (MQ=255)
                                gcgctgactgGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAACCTCGGt                                                                                                       <  2:851122/137‑1 (MQ=255)
                                   ctgactgGCGGTCAGGCCCTGCGACAGGCGAAAGCCGGTCTTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGcc                                                                                                    <  2:267533/140‑1 (MQ=255)
                                    tgactgGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCAGCCAGCATACATACGGGTCAGTCGCGCTATACGGCAAACTCGGTGccc                                                                                                   >  2:303911/1‑139 (MQ=255)
                                     gactgGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAg                                                                                                  <  1:551480/140‑1 (MQ=255)
                                                           cAGGCGAAAGCGGGTGTTGAAGCAGCCTAACTGTCGGAATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATcaacaa                                                                            <  1:416238/140‑1 (MQ=255)
                                                            aGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCTCTCTATCCCGCAGACTGCGTGCCAGGTGTGGCGGAGCGGATCaacaac                                                                           >  2:221576/1‑140 (MQ=255)
                                                               cGAAAGCGGGTATTGAATCAGTTTATCTGTCCGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACAcc                                                                        <  2:260107/140‑1 (MQ=255)
                                                               cGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACAcc                                                                        <  1:509209/140‑1 (MQ=255)
                                                                    gCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCcg                                                                   <  2:48115/140‑1 (MQ=255)
                                                                        gTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCAGATCAAAAACACCTTCcgtcgt                                                               >  1:260107/1‑140 (MQ=255)
                                                                          aTTGAAGCAGTCTATCTGTCGGTATTGCACGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGc                                                             <  1:284556/140‑1 (MQ=255)
                                                                          aTTGAAGCAGTCTATCTGTCGGTATGGCATGTAGCGGAGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGc                                                             <  1:303911/140‑1 (MQ=255)
                                                                               aGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATc                                                        <  2:586872/140‑1 (MQ=255)
                                                                               aGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATc                                                        <  2:771844/140‑1 (MQ=255)
                                                                                  aGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGTTGGGGCGGATCAACAATACCTCCCGTCGTGCCGATCAGa                                                     >  2:334993/1‑140 (MQ=255)
                                                                                    tCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATc                                                   >  1:496814/1‑140 (MQ=255)
                                                                                      tATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCa                                                 >  2:213972/1‑140 (MQ=255)
                                                                                            tCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTc                                           >  1:621635/1‑140 (MQ=255)
                                                                                            tCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTc                                           <  2:496814/140‑1 (MQ=255)
                                                                                              gggATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCcg                                         >  1:809644/1‑140 (MQ=255)
                                                                                              gggATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCcg                                         >  1:492570/1‑140 (MQ=255)
                                                                                                gATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCcgcg                                       >  1:209780/1‑140 (MQ=255)
                                                                                                  tGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCggg                                     >  2:700267/1‑140 (MQ=255)
                                                                                                    gCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCa                                   <  2:621635/140‑1 (MQ=255)
                                                                                                    gCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCa                                   <  1:213972/140‑1 (MQ=255)
                                                                                                          aGCGGCGGACGCTAATCTGGCGACCTGCAAGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGc                             <  1:367210/140‑1 (MQ=255)
                                                                                                            cggcggACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCg                           >  1:52627/1‑140 (MQ=255)
                                                                                                              gcggACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCACATCTAACGGTCCGCGGGCATTGAGCCggg                         >  1:335899/1‑140 (MQ=255)
                                                                                                               cggACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGc                             <  2:141790/135‑1 (MQ=255)
                                                                                                               cggACGCTAACCTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGc                             >  1:141790/1‑135 (MQ=255)
                                                                                                                      tAACCTGGCGGCCAGCATGCATCCGGATCCGTCGCTCTATCCGGCAAACTCGGTGCCAGTTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCcgc                 <  1:346008/140‑1 (MQ=255)
                                                                                                                         ccTGGCGGCCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCt              <  2:209780/140‑1 (MQ=255)
                                                                                                                         ccTGGCGGCCAGCATGCATCCGAATCAGTCGCTCTGTCCGGCAAAATCTGTGCCAGCTGTGGTGGAGCGGATCAACAACCCCTACCGTCGAGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCt              <  2:354637/140‑1 (MQ=255)
                                                                                                                            ggcggcCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATg           <  1:42031/140‑1 (MQ=255)
                                                                                                                              cggcCAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTc         >  2:692692/1‑140 (MQ=255)
                                                                                                                                  cAGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGAtt     >  2:836025/1‑140 (MQ=255)
                                                                                                                                   aGCATGCATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATCGAGCCGCGCGATCCGCGCTATGTCGATTa    >  1:320056/1‑140 (MQ=255)
                                                                                                                                     cATGCATCCGGATCACTCGCTCCATCCGGAAAACTCGCTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACt  <  1:221576/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
AGTCGAAAAAAGGCTACATCAACAGCCTCGGCGCACTGACTGGCGGTCAGGCGCTGCAACAGGCGAAAGCGGGTATTGAAGCAGTCTATCTGTCGGGATGGCAGGTAGCGGCGGACGCTAACCTGGCGGCCAGCATGTATCCGGATCAGTCGCTCTATCCGGCAAACTCGGTGCCAGCTGTGGTGGAGCGGATCAACAACACCTTCCGTCGTGCCGATCAGATCCAATGGTCCGCGGGCATTGAGCCGGGCGATCCGCGCTATGTCGATTACT  >  NC_000913/4217287‑4217559

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