Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,365,608	T→G	intergenic (+57/+94)	yhcE → / ← insH1	pseudogene fragment; interrupted by IS5/IS5 transposase and trans‑activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,365,608	0	T	G	100.0%	160.0 / NA	44	intergenic (+57/+94)	yhcE/insH1	pseudogene fragment; interrupted by IS5/IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (21/23);  total (21/23)

GCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTT  >  NC_000913/3365461‑3365700                                                                                                                                                    |                                                                                             gCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATgtct                                                                                           <  2:346091/151‑1 (MQ=255)    aGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATgtcttc                                                                                         >  2:36156/1‑150 (MQ=255)       gATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGt                                                                                     >  1:132280/1‑151 (MQ=255)            aTTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAg                                                                                <  2:414180/151‑1 (MQ=255)              tttCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGcc                                                                              >  1:177848/1‑151 (MQ=255)                 cAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGGGcc                                                                              <  2:907427/148‑1 (MQ=255)                  aGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGAtt                                                                          >  1:445914/1‑151 (MQ=255)                      agtaACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGAtttttt                                                                      <  2:622874/151‑1 (MQ=255)                         aaCAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTtct                                                                    <  1:534503/150‑1 (MQ=255)                            aaTCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTccag                                                                <  1:138635/151‑3 (MQ=255)                                tcGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTaa                                                             >  1:217315/1‑150 (MQ=255)                                 cGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGt                                                               <  1:84050/147‑1 (MQ=255)                                  gACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATg                                                          >  1:234343/1‑151 (MQ=255)                                        gCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGa                                                    <  2:179771/151‑1 (MQ=255)                                        gCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGa                                                    <  2:414597/151‑1 (MQ=255)                                          gTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAAt                                                  >  2:841503/1‑151 (MQ=255)                                            tAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATc                                                 <  2:706177/150‑1 (MQ=255)                                             aGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATc                                                 >  1:95703/1‑149 (MQ=255)                                                cAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCt                                            >  1:178016/1‑151 (MQ=255)                                                cAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCt                                            >  1:163247/1‑151 (MQ=255)                                                  gTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTAt                                          <  2:844883/151‑1 (MQ=255)                                                  gTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTAt                                          >  2:35950/1‑151 (MQ=255)                                                   tGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTAtt                                         >  1:269924/1‑151 (MQ=255)                                                     gCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTa                                       <  2:43163/151‑1 (MQ=255)                                                     gCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTa                                       <  2:678303/151‑1 (MQ=255)                                                       aaacaTGACAATACGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATAAGCCTAtt                                         >  2:183039/1‑147 (MQ=255)                                                        aaaaTGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGAc                                    <  2:165183/151‑1 (MQ=255)                                                           aTGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGt                                 >  2:340795/1‑151 (MQ=255)                                                              aCAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGttt                               <  2:60314/150‑1 (MQ=255)                                                              aCAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGAGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttc                              <  2:293481/151‑1 (MQ=255)                                                                aaTCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttctt                            >  2:283598/1‑151 (MQ=255)                                                                 aTCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTttcttc                           <  1:59977/151‑1 (MQ=255)                                                                    cGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGcc                        >  1:363915/1‑151 (MQ=255)                                                                     gggCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTCCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCa                       <  1:472001/151‑1 (MQ=255)                                                                         agcaCAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCAttt                    >  1:683477/1‑150 (MQ=255)                                                                           cacaGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTa                   >  1:61607/1‑149 (MQ=255)                                                                           cacaGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTa                   >  2:702049/1‑149 (MQ=255)                                                                                 aaCAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGtt            >  2:40306/1‑150 (MQ=38)                                                                                 aaCAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTa           <  2:158343/151‑1 (MQ=37)                                                                                  aCAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTa           <  2:618115/150‑1 (MQ=37)                                                                                  aCAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTAt          <  1:667337/151‑1 (MQ=37)                                                                                   cAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATc         <  1:603079/151‑1 (MQ=25)                                                                                      cATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcccc      <  2:418419/151‑1 (MQ=18)                                                                                       aTCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcccc      >  1:752088/1‑150 (MQ=18)                                                                                        tCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcccc      >  1:897787/1‑149 (MQ=17)                                                                                          aaTTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATGTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGtt  <  1:1346/151‑1 (MQ=14)                                                                                                                                                    |                                                                                             GCCAGCGATGCATTTTCAGGTAGTAACAATCTCGACCTGGCGTTAGCCAGTGGCAAAAATGACAATCCGGGCAGCACAGGAACAGCATCAATTAGGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTT  >  NC_000913/3365461‑3365700

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑