Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	735,660	T→A	intergenic (‑401/‑308)	gltA ← / → sdhC	citrate synthase/succinate dehydrogenase cytochrome b556 large membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	735,660	0	T	A	100.0%	308.3 / NA	89	intergenic (‑401/‑308)	gltA/sdhC	citrate synthase/succinate dehydrogenase cytochrome b556 large membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (66/23);  total (66/23)

AATGATTGTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGTGACCTAGATCACTGTTCAGGATAAAACCCGACAAACTATATGTAGGT  >  REL606/735612‑735707                                                 |                                                aaTGATTGTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCgaga                                               <  1:2826713/51‑1 (MQ=255) aaTGATTGTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCgaga                                               >  1:2274996/1‑51 (MQ=255) aaTGATTGTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCgaga                                               >  1:3093394/1‑51 (MQ=255) aaTGATTGTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCgaga                                               <  1:4847873/51‑1 (MQ=255)  aTGATTGTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGAc                                              <  1:5154396/51‑1 (MQ=255)  aTGATTGTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGAc                                              >  1:3780365/1‑51 (MQ=255)   tGATTGTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGAcc                                             <  1:6199657/51‑1 (MQ=255)     atagTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTa                                           >  1:850039/4‑51 (MQ=255)     aTTGTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTa                                           >  1:128117/1‑51 (MQ=255)     aTTGTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTa                                           >  1:850041/1‑51 (MQ=255)     aTTGTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTa                                           >  1:1228228/1‑51 (MQ=255)     aTTGTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTa                                           >  1:7082079/1‑51 (MQ=255)     aTTGTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTa                                           >  1:4318576/1‑51 (MQ=255)      ttGTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTAg                                          >  1:1666722/1‑51 (MQ=255)      ttGTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTAg                                          >  1:1099000/1‑51 (MQ=255)      ttGTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTAg                                          >  1:4370455/1‑51 (MQ=255)       tGTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTAGa                                         >  1:7100694/1‑51 (MQ=255)        gTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTAGAt                                        >  1:3982713/1‑51 (MQ=255)        gTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTAGAt                                        <  1:2981025/51‑1 (MQ=255)        gTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTAGAt                                        <  1:401217/51‑1 (MQ=255)         tCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCt                                            >  1:959686/1‑46 (MQ=255)          cAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCa                                      >  1:2731072/1‑51 (MQ=255)          cAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCa                                      <  1:5357634/51‑1 (MQ=255)          cAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCa                                      >  1:88800/1‑51 (MQ=255)            aaTTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACt                                    <  1:1782810/51‑1 (MQ=255)            aaTTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACt                                    >  1:3135641/1‑51 (MQ=255)             aTTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTg                                   <  1:3900937/51‑1 (MQ=255)             aTTAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTg                                   >  1:4199072/1‑51 (MQ=255)              ttAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGt                                  >  1:6059318/1‑51 (MQ=255)              ttAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGt                                  >  1:2207487/1‑51 (MQ=255)              ttAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGt                                  >  1:3772895/1‑51 (MQ=255)              ttAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGt                                  >  1:1236850/1‑51 (MQ=255)               tAGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGtt                                 >  1:2840924/1‑51 (MQ=255)                aGATGATTAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTc                                >  1:3591454/1‑51 (MQ=255)                 gatgatTAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCa                               >  1:6432282/1‑51 (MQ=255)                 gatgatTAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCa                               <  1:1365814/51‑1 (MQ=255)                 gatgatTAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCa                               >  1:5795714/1‑51 (MQ=255)                   tgatTAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAgg                             >  1:3419639/1‑51 (MQ=255)                   tgatTAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAgg                             <  1:7215027/51‑1 (MQ=255)                   tgatTAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAgg                             >  1:102424/1‑51 (MQ=255)                    tatTAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGa                            >  1:6607760/2‑51 (MQ=255)                    gatTAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGa                            <  1:3877090/51‑1 (MQ=255)                      ttAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGATa                          >  1:5965358/1‑51 (MQ=255)                      ttAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGATa                          <  1:4871036/51‑1 (MQ=255)                      ttAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGATa                          <  1:34437/51‑1 (MQ=255)                      ttAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGATa                          >  1:2828675/1‑51 (MQ=255)                      ttAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGATa                          >  1:4027189/1‑51 (MQ=255)                       tAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGATaa                         <  1:4690107/51‑1 (MQ=255)                       tAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGATaa                         >  1:4567358/1‑51 (MQ=255)                       tAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGATaa                         >  1:222114/1‑51 (MQ=255)                       tAAAAATTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGATaa                         <  1:4198983/51‑1 (MQ=255)                        aaaaaTTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGATaaa                        <  1:3200767/51‑1 (MQ=255)                        aaaaaTTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGATaaa                        <  1:3189989/51‑1 (MQ=255)                         aaaaTTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGATaaaa                       >  1:4225039/1‑51 (MQ=255)                         aaaaTTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGATaaaa                       <  1:4478441/51‑1 (MQ=255)                         aaaaTTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGATaaaa                       >  1:1826318/1‑51 (MQ=255)                         aaaaTTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGATaaaa                       <  1:4970828/51‑1 (MQ=255)                          aaaTTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGATAAAAc                      >  1:4811124/1‑51 (MQ=255)                           aaTTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGATAAAAcc                     >  1:2472634/1‑51 (MQ=255)                           aaTTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGATAAAAcc                     <  1:4697450/51‑1 (MQ=255)                            aTTAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGATAAAAccc                    >  1:1169188/1‑51 (MQ=255)                             ttAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGATAAAACCCg                   >  1:4407300/1‑51 (MQ=255)                             ttAAATAAATGTTGTTATCGAGACCTAGATCACTGTTCAGGATAAAACCCg                   >  1:5276548/1‑51 (MQ=255)                              taaataaaTGTTGTTATCGAGACCTAGATCACTGTTCAGGATAAAACCCGa                  >  1:1873517/1‑51 (MQ=255)                               aaataaatGTTGTTATCGAGACCTAGATCACTGTTCAGGATAAAACCCGAc                 >  1:5996456/1‑51 (MQ=255)                                aataaatGTTGTTATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACa                >  1:6844182/1‑51 (MQ=255)                                 ataaatGTTGTTATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACaa               >  1:6522076/1‑51 (MQ=255)                                 ataaatGTTGTTATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACaa               >  1:4002657/1‑51 (MQ=255)                                  taaatGTTGTTATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACaaa              >  1:1019515/1‑51 (MQ=255)                                  taaatGTTGTTATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACaaa              >  1:2152621/1‑51 (MQ=255)                                  taaatGTTGTTATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACaaa              >  1:4661617/1‑51 (MQ=255)                                  taaatGTTGTTATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACaaa              >  1:6566273/1‑51 (MQ=255)                                  taaatGTTGTTATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACaaa              <  1:3356416/51‑1 (MQ=255)                                    aatGTTGTTATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACAAACt            >  1:5175183/1‑51 (MQ=255)                                     atGTTGTTATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACAAACta           >  1:6712786/1‑51 (MQ=255)                                     atGTTGTTATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACAAACta           >  1:1273236/1‑51 (MQ=255)                                      tgttgtTATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACAAACtat          <  1:2831101/51‑1 (MQ=255)                                        ttgttATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACAAACTatat        >  1:935793/1‑51 (MQ=255)                                        ttgttATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACAAACTatat        >  1:2759967/1‑51 (MQ=255)                                         tgttATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACAAACTATATg       >  1:4506832/1‑51 (MQ=255)                                         tgttATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACAAACTATATg       >  1:7268031/1‑51 (MQ=255)                                         tgttATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACAAACTATATg       >  1:2048168/1‑51 (MQ=255)                                         tgttATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACAAACTATATg       >  1:4356958/1‑51 (MQ=255)                                          gttATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACAAACTATATGt      >  1:5014870/1‑51 (MQ=255)                                           ttATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACAAACTATATGTa     <  1:3256876/51‑1 (MQ=255)                                           ttATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACAAACTATATGTa     >  1:6540681/1‑51 (MQ=255)                                           ttATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACAAACTATATGTa     >  1:2954745/1‑51 (MQ=255)                                            tATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACAAACTATATGTAg    <  1:3292315/51‑1 (MQ=255)                                            tATCGAGACCTAGATCACTGTTCAGGATAAAACCCGACAAACTATATGTAg    >  1:3482110/1‑51 (MQ=255)                                             aTCGAGACCTAGATCACTGTTCAGGATAAAACCCGACAAACTATATGTAgg   >  1:4514212/1‑51 (MQ=255)                                              tCGAGACCTAGATCACTGTTCAGGATAAAACCCGACAAACTATATGTAGGt  >  1:6767331/1‑51 (MQ=255)                                              tCGAGACCTAGATCACTGTTCAGGATAAAACCCGACAAACTATATGTAGGt  >  1:3753034/1‑51 (MQ=255)                                              tCGAGACCTAGATCACTGTTCAGGATAAAACCCGACAAACTATATGTAGGt  >  1:3220083/1‑51 (MQ=255)                                                 |                                                AATGATTGTCAAATTAGATGATTAAAAATTAAATAAATGTTGTTATCGTGACCTAGATCACTGTTCAGGATAAAACCCGACAAACTATATGTAGGT  >  REL606/735612‑735707

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