Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,232,858	T→C	intergenic (+38/+117)	ubiA → / ← plsB	4‑hydroxybenzoate octaprenyltransferase/glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,232,858	0	T	C	100.0%	267.4 / NA	63	intergenic (+38/+117)	ubiA/plsB	4‑hydroxybenzoate octaprenyltransferase/glycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/37);  total (26/37)

TTTTTAGGGCTGGCAATGAGTTACTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCTGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACATCGCAATTTATTGAATTTGCAGATTAT  >  REL606/4232785‑4232927                                                                          |                                                                      tttttAGGGCTGGCAATGAGTTACTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCgg                                                                     >  1:615957/1‑76 (MQ=255)  ttttAGGGCTGGCAATGAGTTACTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCggg                                                                    <  2:1524129/76‑1 (MQ=255)  ttttAGGGCTGGCAATGAGTTACTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCggg                                                                    <  2:1674518/76‑1 (MQ=255)  ttttAGGGCTGGCAATGAGTTACTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCggg                                                                    <  2:1276676/76‑1 (MQ=255)  ttttAGGGCTGGCAATGAGTTACTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCggg                                                                    <  2:301711/76‑1 (MQ=255)    ttAGGGCTGGCAATGAGTTACTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGtt                                                                  <  2:154558/76‑1 (MQ=255)    ttAGGGCTGGCAATGAGTTACTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGtt                                                                  <  2:603260/76‑1 (MQ=255)    ttAGGGCTGGCAATGAGTTACTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGtt                                                                  <  2:1598913/76‑1 (MQ=255)    ttAGGGCTGGCAATGAGTTACTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGtt                                                                  <  2:1322259/76‑1 (MQ=255)    ttAGGGCTGGCAATGAGTTACTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGtt                                                                  <  2:461902/76‑1 (MQ=255)    ttAGGGCTGGCAATGAGTTACTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGtt                                                                  <  2:1092324/76‑1 (MQ=255)    ttAGGGCTGGCAATGAGTTACTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGtt                                                                  <  2:341665/76‑1 (MQ=255)      aGGGCTGGCAATGAGTTACTGGCATTTCTGATGATGTAAAAAAGCCGGATGATNATCCGGCTTTCTTCCGGGTTGc                                                                <  2:240826/76‑1 (MQ=255)       gggCTGGCAATGAGTTACTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGcc                                                               <  2:1053750/76‑1 (MQ=255)           tGGCAATGAGTTACTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGAt                                                           <  2:825551/76‑1 (MQ=255)             gCAATGAGTTACTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATgc                                                         <  1:409407/76‑1 (MQ=255)               aaTGAGTTACTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATgcgc                                                       >  1:905035/1‑76 (MQ=255)                      tACTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCtt                                                <  2:1369622/76‑1 (MQ=255)                        cTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTtct                                              >  1:394067/1‑76 (MQ=255)                        cTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTtct                                              <  2:353550/76‑1 (MQ=255)                        cTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTtct                                              >  1:15298/1‑76 (MQ=255)                         tGGCATTTCTGATGATGTAAANGAGCTGGATGATCATCAGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTtctc                                             <  2:279769/76‑1 (MQ=255)                         tGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTtctc                                             >  1:543188/1‑76 (MQ=255)                         tGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTtctc                                             <  2:1713536/76‑1 (MQ=255)                         tGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTtctc                                             <  2:1261563/76‑1 (MQ=255)                         tGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCCtctc                                             >  1:1216347/1‑76 (MQ=255)                         tGGCAGTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTtctc                                             <  1:339779/76‑1 (MQ=255)                          ggCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCa                                            <  2:1685004/76‑1 (MQ=255)                          ggCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCa                                            <  1:1315573/76‑1 (MQ=255)                          ggCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCa                                            <  2:558105/76‑1 (MQ=255)                          ggCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCa                                            <  2:545527/76‑1 (MQ=255)                             aTTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGc                                         >  1:398746/1‑76 (MQ=255)                                tCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTa                                      >  2:9749/1‑76 (MQ=255)                                 cTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTac                                     <  1:114409/76‑1 (MQ=255)                                  tgatgaTGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTaca                                    >  2:1691631/1‑76 (MQ=255)                                   gatgatGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTacac                                   >  2:1392716/1‑76 (MQ=255)                                   gatgatGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGAAGCGCGGCGCTTCTCAGGCCTacac                                   >  2:562353/1‑76 (MQ=255)                                    atgaggTAAAAAAGCCGGATGATCTTCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTacaca                                  <  2:766970/76‑1 (MQ=38)                                         gTAAAAAAGCCGGATGATCCTCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTacacaacaca                             <  1:142288/76‑1 (MQ=255)                                         gTAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTacacaacaca                             <  1:1492514/76‑1 (MQ=255)                                          tAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTacacaacacac                            >  1:913462/1‑75 (MQ=255)                                          tAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACAt                            >  1:511864/1‑76 (MQ=255)                                          tAAAAAAGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCAGCGCTTCTCAGGCCTACACAACACAt                            <  1:546616/76‑1 (MQ=255)                                            aaaaaGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACATCg                          >  2:341005/1‑76 (MQ=255)                                             aaaaGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACATCGc                         >  2:1572950/1‑76 (MQ=255)                                             aaaaGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACATCGc                         >  2:89842/1‑76 (MQ=255)                                              aaaGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACATCGCa                        >  2:350346/1‑76 (MQ=255)                                              aaaGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACATCGCa                        >  2:63309/1‑76 (MQ=255)                                              aaaGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACATCGCa                        <  1:317474/76‑1 (MQ=255)                                              aaaGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACATCGCa                        >  2:1231476/1‑76 (MQ=255)                                              aaaGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACATCGCa                        >  2:1418146/1‑76 (MQ=255)                                               aaGCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACATCGCaa                       >  2:184881/1‑76 (MQ=255)                                                 gCCGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACATCGCAAtt                     <  2:338308/76‑1 (MQ=255)                                                  ccGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACATCGCAAttt                    <  2:821831/76‑1 (MQ=255)                                                  ccGGATGATCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACAACGCAAttt                    <  2:516742/76‑1 (MQ=255)                                                     gatgatCATCCGGCTTTC‑TCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACATCGCAATTTATTg                <  1:1622179/76‑1 (MQ=38)                                                        gatCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACATCGCAATTTATTGaa              <  1:1310042/76‑1 (MQ=255)                                                        gatCATCCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACATCGCAATTTATTGaa              >  1:1131599/1‑76 (MQ=255)                                                          tcatcCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACATCGCAATTTATTGAAtt            <  1:40259/76‑1 (MQ=255)                                                          tcatcCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACATCGCAATTTATTGAAtt            >  2:570917/1‑76 (MQ=255)                                                          tcatcCGGCTTTCTTCCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACATCGCAATTTATTGAAtt            >  2:1058141/1‑76 (MQ=255)                                                                    ttcttcCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACATCGCAATTTATTGAATTTGCAGATTAt  >  2:1539194/1‑76 (MQ=255)                                                                    ttcttcCGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACATCGCAATTTATTGAATTTGCAGATTAt  >  2:1578861/1‑76 (MQ=255)                                                                          |                                                                      TTTTTAGGGCTGGCAATGAGTTACTGGCATTTCTGATGATGTAAAAAAGCCGGATGATCATCCGGCTTTCTTCTGGGTTGCCTGATGCGCGGCGCTTCTCAGGCCTACACAACACATCGCAATTTATTGAATTTGCAGATTAT  >  REL606/4232785‑4232927

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