Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	3,911,973	C→T	intergenic (‑174/‑149)	yifB ← / → ilvL	predicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	3,911,973	0	C	T	100.0%	161.9 / NA	43	intergenic (‑174/‑149)	yifB/ilvL	predicted bifunctional enzyme and transcriptional regulator/ilvG operon leader peptide
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (21/22);  total (21/22)

GCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAA  >  REL606/3911898‑3912045                                                                            |                                                                         gCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTtc                                                                          >  2:617858/1‑76 (MQ=255)    cGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTtcact                                                                       >  1:967801/1‑76 (MQ=255)      cAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTttactca                                                                     <  1:144832/76‑1 (MQ=255)       aaaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTttactcac                                                                    <  1:1605221/76‑1 (MQ=255)        aaaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTtactcact                                                                   >  2:1515529/1‑76 (MQ=255)         aaaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTACTCACta                                                                  >  2:723959/1‑76 (MQ=255)          aaTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTACTCACtat                                                                 <  1:681375/76‑1 (MQ=255)           aTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTACTCACtatt                                                                <  1:1286000/76‑1 (MQ=255)           aTGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATCTTTTTTACTCACtatt                                                                <  2:1575358/76‑1 (MQ=255)              cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTACTCACtatttta                                                             <  2:901628/76‑1 (MQ=255)              cAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTACTCACtatttta                                                             <  2:1192578/76‑1 (MQ=255)                   gACAAAGGATGAACTACGAGGAAGGGAACAACATTCATNCTGAAATTGAATTTTTTTTACTCACTATTTTAttttt                                                        <  1:1995717/76‑1 (MQ=255)                   gACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTACTCACTATTTTAttttt                                                        >  1:1857336/1‑76 (MQ=255)                       aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTACTCACTATTTTATTTTTaaaa                                                    >  2:853254/1‑76 (MQ=255)                       aaGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTACTCACTATTTTATTTTTaaaa                                                    <  1:967797/76‑1 (MQ=255)                         ggATGAACTACGAGGAAGGGAACAACATNGATACTGAAATTGAATTTTATNTACTCACTATTTTATTTTTaaaaaa                                                  >  1:288763/1‑76 (MQ=255)                           aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAaaca                                                >  2:710892/1‑76 (MQ=255)                           aTGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTACTCACTATTCTATTTTTAAAAaaca                                                >  1:1690683/1‑76 (MQ=255)                             gAACTACGAGGAAGGGAGCAACAGTTATACTGAAATTGANTTTTTATTACTCACTNTTTTATATTTAAAAaacaac                                              >  2:2025755/1‑76 (MQ=255)                                cTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAAt                                           >  1:987143/1‑76 (MQ=255)                                 tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAAtt                                          >  1:1943313/1‑76 (MQ=255)                                 tACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAAtt                                          >  2:1086582/1‑76 (MQ=255)                                    gAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTtat                                       >  2:972972/1‑76 (MQ=255)                                     aGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTtata                                      >  1:1129882/1‑76 (MQ=255)                                      ggAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTatat                                     >  2:637072/1‑76 (MQ=255)                                      ggAAGGGAACAACATTCATACTGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTatat                                     >  1:232542/1‑76 (MQ=255)                                        aaGGGAACAACATTCATACTGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTATATTg                                   <  1:451645/76‑1 (MQ=255)                                        aaGGCAAGACCATGCGTACTGCAATTGCATTTCTGTTACTTCCTAACTTTTTTTTAATGAACAACAATTTATATTg                                   <  1:2036494/76‑1 (MQ=255)                                           ggAACAACATTCATACTGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGaaa                                <  1:454282/76‑1 (MQ=255)                                             aacaacATTCATACTGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatt                              <  2:873901/76‑1 (MQ=255)                                             aacaacATTCATACTGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatt                              <  2:846555/76‑1 (MQ=255)                                             aacaacATTCATACTGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAatt                              <  2:1034653/76‑1 (MQ=255)                                               caacaTTCATACTGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattat                            >  2:702382/1‑76 (MQ=255)                                                aacaTTCATACTGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAattatt                           >  1:1932790/1‑76 (MQ=255)                                                 acaTTCATACTGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta                          <  1:1083299/76‑1 (MQ=255)                                                 acaTTCATACTGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAAttatta                          <  2:931894/76‑1 (MQ=255)                                                        tACTGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcat                   >  2:881704/1‑76 (MQ=255)                                                         aCTGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatc                  >  1:422410/1‑76 (MQ=255)                                                           tGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat                <  1:138135/76‑1 (MQ=255)                                                           tGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat                >  2:67812/1‑76 (MQ=255)                                                           tGAAATTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGcatcat                >  1:1242644/1‑76 (MQ=255)                                                             aaaTTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaa              <  1:156316/76‑1 (MQ=255)                                                              aaTTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa             <  2:546680/76‑1 (MQ=255)                                                              aaTTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa             <  2:1867668/76‑1 (MQ=255)                                                              aaTTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa             <  2:1242643/76‑1 (MQ=255)                                                              aaTTGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaa             <  1:765236/76‑1 (MQ=255)                                                                ttGAATTTTTTTTACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATaaaaa           >  2:1490514/1‑76 (MQ=255)                                                                      ttttttTANTCAATATTTTATTTTGAAGAAACACACATTTATATTGAAATTATTAAACGCATCATAAAAATCGGcc     <  2:349173/76‑1 (MQ=255)                                                                      ttttttTACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGcc     <  1:1030946/76‑1 (MQ=255)                                                                       ttttttACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCa    <  1:710893/76‑1 (MQ=255)                                                                         ttttACTCGCTATTTTATTTTTAAAAAACAACGATTTATATTGAAATTATTAAACGCATCANAAAAATCGGCCaaa  >  2:33980/1‑76 (MQ=255)                                                                         ttttACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaa  >  2:1559497/1‑76 (MQ=255)                                                                         ttttACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaa  >  1:1214262/1‑76 (MQ=255)                                                                         ttttACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaa  >  1:1142687/1‑76 (MQ=255)                                                                         ttttACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaa  >  1:195202/1‑76 (MQ=255)                                                                         ttttACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCaaa  >  1:266463/1‑76 (MQ=255)                                                                            |                                                                         GCTCGCAAAAATGCAGCGGACAAAGGATGAACTACGAGGAAGGGAACAACATTCATACTGAAATTGAATTTTTTTCACTCACTATTTTATTTTTAAAAAACAACAATTTATATTGAAATTATTAAACGCATCATAAAAATCGGCCAAA  >  REL606/3911898‑3912045

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