Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 124,241 C→T Q409* (CAG→TAG)  aceE → pyruvate dehydrogenase E1 component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913124,2410CT96.9% 80.6 / NA 33Q409* (CAG→TAG) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (13/18);  minor base G (0/1);  total (14/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACG  >  NC_000913/124108‑124376
                                                                                                                                     |                                                                                                                                       
tCACGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGa                                                                                                                                <  1:712552/143‑1 (MQ=255)
   cGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAATGGCGCCGCCACAGCAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGaaaa                                                                                                                             <  1:666573/143‑1 (MQ=255)
       cagaggaaAATCAACGCTGAATTCAAGAAAGCGCAGGAAACCAAAGGGAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGaa                                                                                                                         <  2:303883/138‑1 (MQ=255)
        cGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAAc                                                                                                                        >  1:234244/1‑143 (MQ=255)
         gaagaaAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAAGGAACa                                                                                                                       >  2:451007/1‑143 (MQ=255)
              aaaTCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGAc                                                                                                                  <  1:36254/143‑1 (MQ=255)
                aTCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATACTTGCACACAGCATTTAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACgg                                                                                                                >  1:575174/1‑143 (MQ=255)
                            ttCAAGGAAGCGCAGGAAACCAGAGGCAACGCGACAGTAAGCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCatat                                                                                                    <  2:812448/143‑1 (MQ=255)
                            ttCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAACTGAACATGGACGGTCTGCGTCatat                                                                                                    >  1:381830/1‑143 (MQ=255)
                             tCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATc                                                                                                   >  1:128159/1‑143 (MQ=255)
                               aaGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCcg                                                                                                 >  1:229737/1‑143 (MQ=255)
                                 gAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGATTACGGCATGGGCGACGCGGCTGCAGGTAAAAACATCGCTCACTAGGTTAAGAAAATGAACATGGCCGGTGTGCGTCATATCcgcg                                                                                               >  1:887161/1‑143 (MQ=255)
                                       gcAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCTTATCCGCGACCGtt                                                                                         >  1:158975/1‑143 (MQ=255)
                                             agccAAACGCAAAGCGCCAGGAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACAGCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAAtg                                                                                   <  1:133110/141‑1 (MQ=255)
                                                cAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGTATTTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGc                                                                                >  1:343499/1‑143 (MQ=255)
                                                       aaaGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAAtgt                                                                                  >  1:529061/1‑134 (MQ=255)
                                                       aaaGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAAtgt                                                                                  <  2:529061/134‑1 (MQ=255)
                                                               aGTAATCCTTGCTCATACCATTAAAGGTTACGGCAGGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGACGGTGTCTGATGCAg                                                                 >  1:872767/1‑143 (MQ=255)
                                                                gTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTCCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGAAGa                                                                >  1:353388/1‑143 (MQ=255)
                                                                       ttggccATACCATTAAGGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACGTCGCGCACTAGGCTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGaa                                                         <  2:448906/138‑1 (MQ=255)
                                                                              ggcgataaaaGGTTACGGCAGGGGCGACGCGGCTGAAGGTAAAAACAGCGCGCACGAGGTTAAGAAAATGAACATGGACGGTGTGCGGCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGc                                                  <  2:311162/136‑1 (MQ=255)
                                                                                    tAAAGGTCACGGCATGGGCGACGCGGCTGAAGGCAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACa                                            <  1:699048/143‑1 (MQ=255)
                                                                                     aaaGGTTCCGGCACGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGGGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACAt                                           <  1:695813/143‑1 (MQ=255)
                                                                                        ggTGACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGTACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCAc                                        <  2:71386/143‑1 (MQ=255)
                                                                                                  tGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGaa                              <  2:229737/143‑1 (MQ=255)
                                                                                                  tGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGaa                              <  1:895808/143‑1 (MQ=255)
                                                                                                     gCGACGCGGCTGAAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCACGACCGTATCAATGTGCCGGTGTCTGATGCAGATACCGAAAACCTGCCGGACATCGCCTTCCAGGACGGt                           >  2:500258/1‑143 (MQ=255)
                                                                                                             gagtaaGGTAAAATCATCGTGCACTAGTTTAAGAAAATGAACATGGACGGCGTGCGTCATATCCGCGACCGTTCCCATGTGCCGGTGTCTTATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAaga                   <  2:88022/139‑1 (MQ=255)
                                                                                                              gtgcgagTAAAAACCTCGCGAACTAGGTTAAGATAATGAACATGGACGTTGTGCGGAATATCAGCGACCGCTTCAATGTGCCGGTGCCTGATGTAGATATTGAAAAGCTGACGTACACCACCTTCCCGGAAGGTTCTGAagag                  <  2:874634/137‑1 (MQ=255)
                                                                                                              cTGCAGGTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAagag                  <  2:618132/143‑1 (MQ=255)
                                                                                                              cTGAAGGTAAAAACATCGCGCACTAGCTTAAGAAAATGAAAAAAAACAAAGGGGGGCATATCCGCGATCTTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAagag                  <  1:170741/143‑1 (MQ=255)
                                                                                                                   ggTAAAAACATCGCGCACTAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGACGGTTCTGAAGAGCATAc             >  1:756091/1‑143 (MQ=255)
                                                                                                                              ggcgcACTAGGTTAAGAAAAGGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCATACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACg  <  2:658838/142‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                       
TCACGATCCGAAGAAAATCTACGCTGCATTCAAGAAAGCGCAGGAAACCAAAGGCAAAGCGACAGTAATCCTTGCTCATACCATTAAAGGTTACGGCATGGGCGACGCGGCTGAAGGTAAAAACATCGCGCACCAGGTTAAGAAAATGAACATGGACGGTGTGCGTCATATCCGCGACCGTTTCAATGTGCCGGTGTCTGATGCAGATATCGAAAAACTGCCGTACATCACCTTCCCGGAAGGTTCTGAAGAGCATACCTATCTGCACG  >  NC_000913/124108‑124376

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: