Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	USA300TCH1516_ALE	629,681	C→A	L84I (CTA→ATA)	USA300HOU_RS02990 →	3‑hexulose‑6‑phosphate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	USA300TCH1516_ALE	629,681	0	C	A	100.0%	172.0 / NA	49	L84I (CTA→ATA)	USA300HOU_RS02990	3‑hexulose‑6‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (25/24);  total (25/24)

AACGAAGGTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTGATTATATTG  >  USA300TCH1516_ALE/629543‑629825                                                                                                                                           |                                                                                                                                                 aaCGAAGGTTTACCAGCAGTAAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTg                                                                                                                                         <  1:129544/148‑1 (MQ=255)  aCGAAGGTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAGTGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCCCGGATGTAATTACAATAATAGGTGTTGc                                                                                                                                        >  1:552038/1‑148 (MQ=255)     aaGGTTTACCAGCAGTTAAACATAAGGCCTATGTTGTGAGTAATGTAAAAGTATTAGCCGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGGGCGGATGGAATTACAATAATAGGGGTTGCgga                                                                                                                                     >  1:301779/1‑148 (MQ=255)        gTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCagaaga                                                                                                                                  <  1:1022271/148‑1 (MQ=255)        gTTTACCAGCAGTTAAACATATGGAATACAACATTAGTAATGAAAAAGTATTAGGAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCagaaga                                                                                                                                  <  1:378447/148‑1 (MQ=255)           tACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGTTGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGc                                                                                                                               <  1:917682/148‑1 (MQ=255)           tACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGc                                                                                                                               <  1:510831/148‑1 (MQ=255)               agcagTTAAACATATGGCAGACAACAATATTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCa                                                                                                                           <  1:828339/148‑1 (MQ=255)                             tGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTAt                                                                                                             >  1:466163/1‑148 (MQ=255)                               gCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGATGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTg                                                                                                           <  1:509715/148‑1 (MQ=255)                                cAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTga                                                                                                          <  1:101823/148‑1 (MQ=255)                                cAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTga                                                                                                          >  1:158083/1‑148 (MQ=255)                                cAGACAACATTAGTAATGAAAAAGGATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTga                                                                                                          <  1:1008021/148‑1 (MQ=255)                                     aacaTTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGaagaag                                                                                                     <  1:827616/148‑1 (MQ=255)                                     aacaTTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGaagaag                                                                                                     <  1:827455/148‑1 (MQ=255)                                     aacaTCAGTAATGTAAAAACATCAGCAGACATGAAAATTATAGAAACAACTGATTATGAAGCTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGaagaag                                                                                                     <  1:705721/148‑1 (MQ=255)                                         ttAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGAGGCATCAATTAAAGCAGCTATTGAAGAAGCTCa                                                                                                 >  1:485922/1‑148 (MQ=255)                                            gTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAGGCTCATaa                                                                                              >  1:986188/1‑148 (MQ=255)                                              aaTGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATaaaa                                                                                            >  1:457207/1‑148 (MQ=255)                                                 gTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCGATTGAAGAAGCTCATAAAaata                                                                                         <  1:929995/148‑1 (MQ=255)                                                     aaGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAAtaataa                                                                                     >  1:999448/1‑148 (MQ=255)                                                        tATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACa                                                                                  <  1:917237/148‑1 (MQ=255)                                                         aTTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACaa                                                                                 >  1:636139/1‑148 (MQ=255)                                                                        aaTTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATg                                                                  >  1:1085705/1‑148 (MQ=255)                                                                          ttATAGATGCAGCTGATTATGAGCTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGAt                                                                >  1:622809/1‑148 (MQ=255)                                                                          ttATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGAt                                                                >  1:1070575/1‑148 (MQ=255)                                                                          ttATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGAt                                                                >  1:770221/1‑148 (MQ=255)                                                                           tataGATGCAGCTGATTATGAAGTTAGCAAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGAtt                                                               <  1:848700/148‑1 (MQ=255)                                                                           tataGATGCAGCTGATTATGAAGTAAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGAtt                                                               <  1:461605/148‑1 (MQ=255)                                                                                    aGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATCAACAATTACCAGTTGATATGATTGCTGTTCaa                                                      <  1:629354/148‑1 (MQ=255)                                                                                                aGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGaaaaa                                          >  1:276765/1‑148 (MQ=255)                                                                                                aGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATAGGATTGCTGTTCAAGATTTAGaaaaa                                          >  1:308279/1‑148 (MQ=255)                                                                                                    aGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTg                                      <  1:1059358/148‑1 (MQ=255)                                                                                                        aaGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAATAACGTGCaaa                                  <  1:905240/148‑1 (MQ=255)                                                                                                        aaGCAATTAAATTTGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCaaa                                  <  1:431396/148‑1 (MQ=255)                                                                                                                    ttGGCGCGGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGa                      <  1:430136/148‑1 (MQ=255)                                                                                                                        cgcgGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATg                  >  1:692271/1‑148 (MQ=255)                                                                                                                        cgcgGATGTAATTACAATAATAGGAGCTGCAGAAGATGCAAAACTTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATg                  >  1:795119/1‑148 (MQ=255)                                                                                                                        cgcgGATGTAAGTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGAATTAGAAAAACGTGCAAAAGAACTAGATGAAATg                  >  1:853146/1‑148 (MQ=255)                                                                                                                        cgcgGATGTAAATACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGACTTAGAAAAACGTGCAAAAGAACTAGATGAAATg                  >  1:853141/1‑148 (MQ=255)                                                                                                                         gcgGATGTAATTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATgg                 >  1:406268/1‑148 (MQ=255)                                                                                                                                 aaTTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTGAt         <  1:524860/148‑1 (MQ=255)                                                                                                                                 aaTTACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTGAt         <  1:608963/148‑1 (MQ=255)                                                                                                                                   ttACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATACACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTGATta       >  1:790718/1‑148 (MQ=255)                                                                                                                                   ttACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTGAttt       >  1:723638/1‑147 (MQ=255)                                                                                                                                   ttACAATAATAGGTGTTGCAGAAAATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTGATta       >  1:746556/1‑148 (MQ=255)                                                                                                                                    tACAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTGATtat      <  1:421565/148‑1 (MQ=255)                                                                                                                                     aCAATAATAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAAAGGGTGATGATtata     >  1:883505/1‑148 (MQ=255)                                                                                                                                     aCAATAATAGGTGTTGCAGAAGATGCAACAATTAAAGCAGCTATTGAAGAAGATCATAAAAATAATAAACAGTTACTAGTGGATATGATTGCTGTGCAAGATATAGAAAAACGGGCAAAAGAACTAGATGAAATGGTTGCTCATtata     >  1:881385/1‑148 (MQ=255)                                                                                                                                        ataataGGTGTTGCAGAAGATGCGTCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGCAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTGATTATATTg  <  1:794590/148‑1 (MQ=255)                                                                                                                                        ataataGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTGATTATATTg  <  1:525295/148‑1 (MQ=255)                                                                                                                                        ataataGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTGATTATATTg  <  1:771319/148‑1 (MQ=255)                                                                                                                                        ataataGGTGTTGAAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTGATTATATTg  <  1:1109972/148‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                 AACGAAGGTTTACCAGCAGTTAAACATATGGCAGACAACATTAGTAATGTAAAAGTATTAGCAGACATGAAAATTATAGATGCAGCTGATTATGAAGTTAGCCAAGCAATTAAATTTGGCGCGGATGTAATTACAATACTAGGTGTTGCAGAAGATGCATCAATTAAAGCAGCTATTGAAGAAGCTCATAAAAATAATAAACAATTACTAGTTGATATGATTGCTGTTCAAGATTTAGAAAAACGTGCAAAAGAACTAGATGAAATGGGTGCTGATTATATTG  >  USA300TCH1516_ALE/629543‑629825

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