Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_tpiA 264,098 G→T 100% L50L (CTC→CTA ykfH ← uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA264,0980GT100.0% 92.2 / ‑6.7 38L50L (CTC→CTAykfHuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (24/12);  total (25/13)
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATTT  >  NC_000913_3_pae_tpiA/263961‑264243
                                                                                                                                         |                                                                                                                                                 
aGGCGCAGGGTAATCTCCCGCTGCAGGCCCCTGGTAGGATAGCTCATTACAGTTCTCCACTATTGTCAGGTCAGGGGGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCGGGc                                                                                                                                        >  1:514901/1‑149 (MQ=255)
     cAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGatat                                                                                                                                   >  1:177974/1‑149 (MQ=255)
            aTCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTATACAATGCTGTCAGTTCAGGGTGATGCTATTCAGGCAGGTATAGGGCCCATTGCGGTCCTGGGGGCATTCGGCGTATACCGCTAGTAGTCACGCGATATCCGGaat                                                                                                                            >  2:482619/1‑148 (MQ=255)
             tctcCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTGCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACg                                                                                                                           >  1:86025/1‑149 (MQ=255)
             tctcCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCACATCAGGCAGGTATAGGGCCCATTGCGGTCCTCGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCAATATCCGGAACg                                                                                                                           >  1:308998/1‑149 (MQ=255)
                cccGTTGCCGGCCCCTGCCATGGCTGGTCATGAGTGTGTTCGACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTcc                                                                                                                        <  1:705058/149‑1 (MQ=255)
                cccGCTCGAGGCTCCTGGTAGGGTTGCCAATTACAGTTCTCCTCTATTGGCAGTTTAGGGTGGTGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTcc                                                                                                                        <  2:380763/149‑1 (MQ=255)
                 ccGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTCTTGTCAGTTCAGGGTTATGCTCATCAGGCAGGTCTAGGGCCCATTGCGGTCCTGGCGGCGTTCGCCGTATACCGCTCAGACTCCGGCGGGACCCGGAACGTccc                                                                                                                       >  2:604802/1‑149 (MQ=255)
                            cccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCGGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTaa                                                                                                            >  1:411731/1‑149 (MQ=255)
                            cccTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTaa                                                                                                            <  1:641846/149‑1 (MQ=255)
                                             aTTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACGGACGCTACCGTg                                                                                           >  1:190474/1‑149 (MQ=255)
                                               tACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGACCACCACAGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGcc                                                                                         <  1:722963/149‑1 (MQ=21)
                                                    ttCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCGTTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTgg                                                                                    >  1:91822/1‑149 (MQ=255)
                                                    ttCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGCGTAATGACAGACGCTACCGTCCCACTgg                                                                                    >  2:542394/1‑149 (MQ=255)
                                                       tcCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTAt                                                                                 >  1:612358/1‑149 (MQ=255)
                                                       tcCACTATTGTCAGTTCAAGGTGATGCTCATCAGGCAGGTAACGCGCCCGTTGAGGTCCTGGCGGCGTTCGGCGTATACCCCTAGGACTCCGGCGATAACCGGACCGTCCCTGCCGGTGTGCTCACCGCCGCTACCGCCCCATTGGTAt                                                                                 >  2:304288/1‑149 (MQ=255)
                                                         cACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTAGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTAttt                                                                               >  1:367856/1‑149 (MQ=255)
                                                             tgtgtgAGTGCAGGGTGAGGCTAGTCGGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCg                                                                           <  2:647554/147‑1 (MQ=21)
                                                               tGTCAGTTCAGGGTGATGCTCATCAGGCAGGTCTTGCATGCATGTCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGCACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGt                                                                         >  1:356312/1‑149 (MQ=21)
                                                                 tCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGc                                                                       >  1:49554/1‑149 (MQ=255)
                                                                 tCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGACAGTGc                                                                       >  2:196193/1‑149 (MQ=255)
                                                                         gggTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGCCGATATTCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGa                                                               <  1:468368/149‑1 (MQ=255)
                                                                                 cTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACGCCGGCGATATCCGGAACGTCCCTGCCGGGGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCg                                                       >  1:69908/1‑149 (MQ=255)
                                                                                              tataGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGc                                          >  1:661873/1‑149 (MQ=255)
                                                                                                tcgggCCCTTTGCGGTCCTGGCGGTGTTGGGCGCATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTg                                        <  1:160309/147‑1 (MQ=17)
                                                                                                   ggCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCg                                     >  1:733352/1‑149 (MQ=255)
                                                                                                   ggCCCATTGCGGTCCTGGCGCCGTTCGGAGTATCCCGGGAGGACTCAGGCGATATCCGGAACGTACCGGCCGGTGTAGTGACAGACGCTAGCGTGACACAGGCATTGGCCGGCGCAGTAGCGAAAGATTCGGGGCGCAGGGTGCTGGca                                     >  2:584594/1‑148 (MQ=255)
                                                                                                          tGCGGTCCTGGCGGCGTTCGGCGTAGACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCAGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATc                              <  1:663343/149‑1 (MQ=255)
                                                                                                               tCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCAt                         >  1:432828/1‑149 (MQ=255)
                                                                                                               tCCTGGCGGCGTTCGGCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCAt                         >  1:447237/1‑149 (MQ=255)
                                                                                                                     cggggTTCGGCGTATACCGCTAGGCCGCCGGCGATATCCGTAACGCCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATCCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCt                   <  2:571999/149‑1 (MQ=17)
                                                                                                                            cggcggATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAACGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTa            <  1:29138/149‑1 (MQ=255)
                                                                                                                             ggCGTATACCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTAc           >  1:475392/1‑149 (MQ=255)
                                                                                                                              gggTACACCGCTAGGTCTCCGGTGAGCTCCGGAACGGCCCTGCCGGTGTAATGACATACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACt          <  2:656072/147‑1 (MQ=14)
                                                                                                                                gTATGCCGCGGGGACGCCGGCGCTATCCGGAACGTCCCTGCCGGTGTAAGGACAGACGCTACAGTGCAGCTGGTACTTCCCGGCCTGGTGTTGAAAGATGCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGa        <  1:353036/149‑1 (MQ=14)
                                                                                                                                   tcccGCTAGGACTGCGGCGAAATCCGGAACGTCCCGCCCGGTGTAATGACAGAGGCTACCGTGCCTCTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATaa     <  1:161779/147‑1 (MQ=14)
                                                                                                                                    aCCGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCACGATGCAGGGCGCCTATCCTCATTGCCCAGCGTCTACTGATAAt    >  2:285150/1‑149 (MQ=14)
                                                                                                                                     gctagaggagTCGGGGGTTATGCGGGACGTTCCTGCCGGTGTAATGACAGACGCTACCGTGCCGATGGTATTTGACGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAtt   <  2:104845/145‑1 (MQ=14)
                                                                                                                                     ccGCTAGGACTCCGGCGATATCCGGAACGTCCCTGCAGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGCTAAt  >  2:559175/1‑148 (MQ=255)
                                                                                                                                      cgctAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGACAGATTCGGGACTGAGGATACTCGCGGTATACCGTTACCGCCCTGCGGTGCCTGATAGttt  >  2:299629/1‑149 (MQ=17)
                                                                                                                                      cgcaAGCACTGCGGCGAGATCCGTAACGTCCCTGCCGGTGCAATGACACACGCTCCCGTGCCACTGGTATTTTCCCGTGCAGTAGCGAAAGATTCGTGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAAttt  <  1:406570/149‑1 (MQ=14)
                                                                                                                                         |                                                                                                                                                 
AGGCGCAGTGTAATCTCCCGCTGCAGGCCCCTGGTAGGGTTGCTCATTACAGTTCTCCACTATTGTCAGTTCAGGGTGATGCTCATCAGGCAGGTATAGGGCCCATTGCGGTCCTGGCGGCGTTCGGCGTATACCGCGAGGACTCCGGCGATATCCGGAACGTCCCTGCCGGTGTAATGACAGACGCTACCGTGCCACTGGTATTTGCCGGTGCAGTAGCGAAAGATTCGGGACTCAGGATGCTGGCGGTATATCGTCATTGCCCTGCGTTTACTGATAATTT  >  NC_000913_3_pae_tpiA/263961‑264243

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