Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 123,050 A→T I12F (ATC→TTC)  aceE → pyruvate dehydrogenase E1 component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913123,0500AT97.3% 106.0 / NA 37I12F (ATC→TTC) aceEpyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (16/20);  minor base G (0/1);  total (16/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 = 1.00e+00

GCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGCACAGGTATCA  >  NC_000913/122912‑123188
                                                                                                                                          |                                                                                                                                          
gcgcGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGa                                                                                                                                        >  1:58612/1‑143 (MQ=255)
    ggCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACt                                                                                                                                    >  1:156941/1‑143 (MQ=255)
        aCTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTcgcg                                                                                                                                >  1:404811/1‑143 (MQ=255)
         cTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGa                                                                                                                               <  2:407391/143‑1 (MQ=255)
             aCGTCGAACCTGGCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTgg                                                                                                                           <  1:11053/143‑1 (MQ=255)
                         gctcaTTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATc                                                                                                               <  2:110762/139‑1 (MQ=255)
                              ttGAGCTTTCCCGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATc                                                                                                          <  1:424580/143‑1 (MQ=255)
                                 aGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGCAACTCGCGACTCGCAACAGGCGATCGACTCgct                                                                                                       >  2:95492/1‑141 (MQ=255)
                                  gCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTc                                                                                                      <  2:551164/143‑1 (MQ=255)
                                        cGGCGAGAGTTCAATGGGACAGTTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGCCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGt                                                                                                <  2:611657/143‑1 (MQ=255)
                                         ggCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTg                                                                                               >  1:168218/1‑143 (MQ=255)
                                                    aaTGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGtt                                                                                    <  1:344964/143‑1 (MQ=255)
                                                    aaTGGGACAGGTGCAAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGtt                                                                                    <  1:357078/143‑1 (MQ=255)
                                                       gggACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAACGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAg                                                                                 <  1:511360/143‑1 (MQ=255)
                                                            aGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGCCTGGCTCCCGGCGATCGACGCGGTCAGCCGTGAAGACGGCGTTGCGCGTGc                                                                            >  2:687018/1‑143 (MQ=255)
                                                            aGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGc                                                                            >  1:325215/1‑143 (MQ=255)
                                                                  cAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGACATCCGAAAGGAAGGTGTTGAGCGTGCTCAGTa                                                                      <  2:156941/143‑1 (MQ=255)
                                                                        aCTCAACATTATTAGATAGATAAGAGAGAACCCACGTCAGAACGGTGCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGAt                                                                <  1:319124/143‑1 (MQ=255)
                                                                         cTCAACGTTAGTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATc                                                               <  1:72887/143‑1 (MQ=255)
                                                                           cAACGTTATTGCATATATAAGGAATACCCCATGTCAGAACGTTTCCCAGATGGCGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGa                                                             <  2:319133/143‑1 (MQ=255)
                                                                             aCGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGAcc                                                           <  1:35378/143‑1 (MQ=255)
                                                                                ttattaGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAAc                                                        <  1:207554/143‑1 (MQ=255)
                                                                                ttatcaAATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAAc                                                        <  2:623301/143‑1 (MQ=255)
                                                                                       ataggtgAGGAATAGCCCATGTCGGAACGTTTCCCAAATGACCTGCGTCCTGTCGAAACGCGCGACTGGCTCCAGGCGATCGACGCGGTCAGCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACtgcttgc                                                 <  2:286716/143‑1 (MQ=255)
                                                                                            taAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAg                                            <  1:309704/143‑1 (MQ=255)
                                                                                                gAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTCCTTTCTGAACCCCg                                        >  2:109244/1‑143 (MQ=255)
                                                                                                            tCAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTTAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGgtgt                            >  2:424087/1‑143 (MQ=255)
                                                                                                             cAGAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGACCGAATCGGTCATACGTAAAGAAAGTATTGAACGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTa                           >  1:206908/1‑143 (MQ=255)
                                                                                                              aGAACGTTTCCCAAATGACGTGGATCCGTTCGAACCTCGCGCCTGGCCCCACGCGATCGAATCGGTCATCCGTGAAGAAGGTGTCGAGCGTGCCCAGTATCCGATCGGCCACCTGCCTGCTCACCCCCCCAAAGGcagcgtcc                          >  2:700657/1‑136 (MQ=255)
                                                                                                               gAACGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTaaa                         <  1:383615/143‑1 (MQ=255)
                                                                                                                aaCGTTTCCCAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACGGCTTGCTGAAGCCCGCAAAGGCGGTGTAAAc                        >  1:70236/1‑143 (MQ=255)
                                                                                                                        ccAAATGACGTGGATCCGTTCGAAGCTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGc                <  2:447978/143‑1 (MQ=255)
                                                                                                                        ccAAATGACGTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGc                >  1:19331/1‑143 (MQ=255)
                                                                                                                                cGTGGCTTCGTTCGAAGTGCGCGGCTGGCTCCAGTGGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTTAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGCACAg        <  2:608797/143‑1 (MQ=255)
                                                                                                                                 gTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAGAGGCGGTGTAGACGTAGCCGCGGGCACAgg       >  1:76716/1‑143 (MQ=255)
                                                                                                                                 gTGGATCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGCACAgg       >  1:337595/1‑143 (MQ=255)
                                                                                                                                      tCCGTTCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGCACAGGTATCa  >  1:33282/1‑143 (MQ=255)
                                                                                                                                          |                                                                                                                                          
GCGCGGCAACTAAACGTAGAACCTGTCTTATTGAGCTTTCCGGCGAGAGTTCAATGGGACAGGTTCCAGAAAACTCAACGTTATTAGATAGATAAGGAATAACCCATGTCAGAACGTTTCCCAAATGACGTGGATCCGATCGAAACTCGCGACTGGCTCCAGGCGATCGAATCGGTCATCCGTGAAGAAGGTGTTGAGCGTGCTCAGTATCTGATCGACCAACTGCTTGCTGAAGCCCGCAAAGGCGGTGTAAACGTAGCCGCAGGCACAGGTATCA  >  NC_000913/122912‑123188

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