Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,273,299	A→T	intergenic (+227/‑115)	nrfG → / → gltP	heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/glutamate/aspartate:proton symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,273,299	0	A	T	100.0%	150.9 / NA	39	intergenic (+227/‑115)	nrfG/gltP	heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfG/glutamate/aspartate:proton symporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (19/20);  total (19/20)

CGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCG  >  REL606/4273226‑4273369                                                                          |                                                                       cgTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTc                                                                      <  2:1019339/76‑1 (MQ=255) cgTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTc                                                                      <  1:906141/76‑1 (MQ=255)    tattatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATa                                                                   >  2:659395/1‑76 (MQ=255)     attatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATaa                                                                  >  1:892741/1‑76 (MQ=255)      ttatGCATTAATGGAGCGAANAGCTCTGTTCTTAAANGGTTGCGCAACATCCGGCGCAAATGATACTGTTCATAAg                                                                 <  2:1205120/76‑1 (MQ=255)      ttatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAg                                                                 >  1:499042/1‑76 (MQ=255)      ttatGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAg                                                                 <  1:582471/76‑1 (MQ=255)        atGCATTAATGCAGCGAAAAGCTCTGTTGTTGAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCg                                                               >  2:555726/1‑76 (MQ=255)         tGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGt                                                              >  1:886178/1‑76 (MQ=255)           cATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGGTGCGCANCATACCGCGCAAATGATANTGNTCATAAGCGTTa                                                            >  1:987607/1‑76 (MQ=255)              tAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaa                                                         <  1:233527/76‑1 (MQ=255)               aaTGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaa                                                        >  2:514767/1‑76 (MQ=255)                        aaaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTNAAAAAAATCTACaa                                               <  1:5276/76‑1 (MQ=255)                          aaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAAc                                             >  1:354062/1‑76 (MQ=255)                           aGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAAcc                                            >  1:410521/1‑76 (MQ=255)                              tctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAAc                                         <  1:1098140/76‑1 (MQ=255)                               ctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACg                                        >  2:903326/1‑76 (MQ=255)                               ctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACg                                        >  2:1063262/1‑76 (MQ=255)                                 gttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCa                                      <  1:612303/76‑1 (MQ=255)                                 gttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCa                                      >  2:714643/1‑76 (MQ=255)                                  ttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCaa                                     <  2:310295/76‑1 (MQ=255)                                      tAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAAcaca                                 <  1:202688/76‑1 (MQ=255)                                         aGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACAGAAtt                              <  2:585930/76‑1 (MQ=255)                                         aGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAtt                              >  2:1165272/1‑76 (MQ=255)                                         aGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAtt                              >  2:536628/1‑76 (MQ=255)                                            gTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCAt                           <  2:170132/76‑1 (MQ=255)                                                 gcAACATACCGCGCNAATGATACTGTTCATAAGCGTTCAAAAAATCTACAAACCAACGCAACACAATTCATGCCCt                      >  1:5415/1‑76 (MQ=255)                                                  cAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTg                     <  2:797158/76‑1 (MQ=255)                                                   aaNATACCNCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCNgg                    >  2:989812/1‑76 (MQ=255)                                                   aaCATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTgg                    <  2:753504/76‑1 (MQ=255)                                                      aTACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAg                 <  1:495007/76‑1 (MQ=255)                                                       tACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGt                >  2:425116/1‑76 (MQ=255)                                                       tACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGt                <  1:523303/76‑1 (MQ=255)                                                        aCCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTa               <  2:895529/76‑1 (MQ=255)                                                        aCCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTa               <  1:835809/76‑1 (MQ=255)                                                          cgcgCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATg             >  2:667989/1‑76 (MQ=255)                                                                aaTGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGt       >  2:441956/1‑76 (MQ=255)                                                                 aTGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGtt      <  2:897050/76‑1 (MQ=255)                                                                    aTACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTtctc   >  2:563043/1‑76 (MQ=255)                                                                     tACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTcg  <  2:892741/76‑1 (MQ=255)                                                                          |                                                                       CGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCG  >  REL606/4273226‑4273369

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