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%	196.3 / NA	50	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 (25/25);  total (25/25)

CGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGC  >  REL606/4273226‑4273370                                                                          |                                                                        cgTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTc                                                                       <  2:393932/76‑1 (MQ=255) cgTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTc                                                                       <  1:422741/76‑1 (MQ=255)   ttattaTGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCAt                                                                     <  1:545014/76‑1 (MQ=255)   ttactaTGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCAt                                                                     >  1:1031440/1‑76 (MQ=255)        atGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCg                                                                >  1:1075459/1‑76 (MQ=255)          gCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGtt                                                              >  1:852334/1‑76 (MQ=255)            aTTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGNTCATAAGCGTTaa                                                            >  2:1273874/1‑76 (MQ=255)              tAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaa                                                          >  1:574466/1‑76 (MQ=255)                aTGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaaa                                                        <  2:522504/76‑1 (MQ=255)                aTGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaaa                                                        <  2:431081/76‑1 (MQ=255)                aTGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTaaaaaa                                                        <  1:425805/76‑1 (MQ=255)                   cAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATc                                                     >  2:299860/1‑76 (MQ=255)                    aGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCt                                                    >  1:970385/1‑76 (MQ=255)                     gCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTa                                                   <  2:510414/76‑1 (MQ=255)                      cGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTAc                                                  >  2:1400430/1‑76 (MQ=255)                       gAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGTGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACa                                                 <  1:359582/76‑1 (MQ=255)                        aaaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACaa                                                >  1:1032550/1‑76 (MQ=255)                         aaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACaaa                                               >  2:291383/1‑76 (MQ=255)                         aaaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACaaa                                               <  2:528685/76‑1 (MQ=255)                         aaaGCTCTGTTGNTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACaaa                                               >  1:21436/1‑76 (MQ=255)                          aaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAAc                                              >  1:770235/1‑76 (MQ=255)                          aaGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAAc                                              >  1:961437/1‑76 (MQ=255)                           aGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAAcc                                             <  1:889097/76‑1 (MQ=255)                               ctGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACg                                         >  1:1248324/1‑76 (MQ=255)                                 gttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCa                                       >  1:325436/1‑76 (MQ=255)                                 gttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCa                                       <  2:320750/76‑1 (MQ=255)                                 gttgttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCa                                       <  1:1031330/76‑1 (MQ=255)                                    gttAAAGGGTTGCGGAACATATCGCGCAAATGATACTGTTCATAAGCGTTAAATAAATATACAAACCAACGCAaca                                    <  2:215517/76‑1 (MQ=255)                                    gttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAaca                                    <  2:943386/76‑1 (MQ=255)                                    gttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAaca                                    >  1:383571/1‑76 (MQ=255)                                     ttAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCaacacaacac                                   <  1:86185/76‑1 (MQ=255)                                      tAAAGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAAcaca                                  >  2:428015/1‑76 (MQ=255)                                        naGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAANCAACGCAACACAAt                                >  2:199927/2‑76 (MQ=255)                                        aaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACCCAAt                                <  2:1383190/76‑1 (MQ=255)                                        aaGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAt                                >  2:599510/1‑76 (MQ=255)                                         aGGGTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAtt                               >  1:860025/1‑76 (MQ=255)                                          gggTTGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTc                              <  1:100956/76‑1 (MQ=255)                                              tGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAtt                               <  1:414730/71‑1 (MQ=255)                                              tGCGCAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAAtt                               >  2:414730/1‑71 (MQ=255)                                                 gcAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACAGAATTCATGCCCt                       <  1:316815/76‑1 (MQ=255)                                                 gcAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCt                       >  1:594502/1‑76 (MQ=255)                                                 gcAACATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCt                       >  1:1179193/1‑76 (MQ=255)                                                     cATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGcn                   <  1:644220/76‑2 (MQ=255)                                                     cATACCGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATNCATGCCCTGGCa                   <  1:12555/76‑1 (MQ=255)                                                       tACCGCGCAAATGATACTGTTCATAAGCGTAAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCCGGCAGt                 >  1:1277767/1‑76 (MQ=255)                                                         ccGCGCAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTAt               <  2:834721/76‑1 (MQ=255)                                                             gcAAATGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACAGAATTCATGCCNTGGCAGTATgacg           >  1:184708/1‑73 (MQ=255)                                                                aaTGATACTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGt        <  2:251501/76‑1 (MQ=255)                                                                    aTACTGTTCATAAGCGTTAAAAAAAGCTNCAAACCAACTAAGCACAACTCATGCCCTGGCAGTATGTCACGTtctc    <  1:132190/76‑1 (MQ=255)                                                                     tACTGTTCATAAGCGTTAAAAAAATCTACAAANCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTNTcg   >  2:1301701/1‑76 (MQ=255)                                                                      aCTGTTCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTcgc  <  2:753435/76‑1 (MQ=255)                                                                          |                                                                        CGTTATTATGCATTAATGCAGCGAAAAGCTCTGTTGTTAAAGGGTTGCGCAACATACCGCGCAAATGATACTGATCATAAGCGTTAAAAAAATCTACAAACCAACGCAACACAATTCATGCCCTGGCAGTATGTCACGTTCTCGC  >  REL606/4273226‑4273370

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