Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	621,944	C→A	intergenic (+29/‑242)	unknown → / → unknown	Excinuclease ABC subunit A/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	621,944	0	C	A	100.0%	106.5 / NA	39	intergenic (+29/‑242)	unknown/unknown	Excinuclease ABC subunit A/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (18/21);  total (18/21)

GGATTGTAGATCTCGGCCCTGAAGGCGGCAGCGGCGGCGGACAAATTATCGCGACCGGCACACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGCTCTGCACCCCAAAACCT‑‑‑‑‑‑GGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAG  >  1/621782‑622174                                                                                                                                                                   |                                                                                                                                                                                                                                             ggATTGTAGATCTCGGCCCTGAAGGCGGCAGCGGCGGCGGACAAATTATCGCGACCGGCACACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaa                                                                                                                                                       <  1:163099‑M2/250‑71 (MQ=255) ggATTGTAGATCTCGGCCCTGAAGGCGGCAGCGGCGGCGGACAAATTATCGCGACCGGCACACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaa                                                                                                                                                       <  1:792303‑M2/250‑71 (MQ=255) ggATTGTAGATCTCGGCCCTGAAGGCGGCAGCGGCGGCGGACAAATTATCGCGACCGGCACACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaa                                                                                                                                                       <  1:261979‑M2/250‑71 (MQ=255) ggATTGTAGATCTCGGCCCTGAAGGCGGCAGCGGCGGCGGACAAATTATCGCGACCGGCACACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaa                                                                                                                                                       <  1:211211‑M2/250‑71 (MQ=255)  gATTGTAGATCTCGGCCCTGAAGGCGGCAGCGGCGGCGGACAAATTATCGCGACCGGCACACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaac                                                                                                                                                      >  1:98741‑M2/1‑179 (MQ=255)  gATTGTAGATCTCGGCCCTGAAGGCGGCAGCGGCGGCGGACAAATTATCGCGACCGGCACACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaac                                                                                                                                                      >  1:470985‑M2/1‑179 (MQ=255)   aTTGTAGATCTCGGCCCTGAAGGCGGCAGCGGCGGCGGACAAATTATCGCGACCGGCACACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacag                                                                                                                                                    >  1:699774‑M2/1‑178 (MQ=255)         gATCTCGGCCCTGAAGGCGGCAGCGGCGGCGGACAAATTATCGCGACCGGCACACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaa                                                                                                                                                        <  1:896379‑M2/241‑70 (MQ=255)          aTCTCGGCCCTGAAGGCGGCAGCGGCGGCGGACAAATTATCGCGACCGGCACACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccg                                                                                                                                             <  1:492406‑M2/251‑81 (MQ=255)          aTCTCGGCCCTGAAGGCGGCAGCGGCGGCGGACAAATTATCGCGACCGGCACACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccg                                                                                                                                             >  1:625345‑M2/1‑171 (MQ=255)                  ccTGAAGGCGGCAGCGGCGGCGGACAAATTATCGCGACCGGCACACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggat                                                                                                                                     >  1:711596‑M2/1‑163 (MQ=255)                     gAAGGCGGCAGCGGCGGCGGACAAATTATCGCGACCGGCACACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggatttt                                                                                                                                  <  1:385843‑M2/251‑92 (MQ=255)                         gcggcAGCGGCGGCGGACAAATTATCGCGACCGGCACACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttggt                                                                                                                               >  1:891314‑M2/1‑156 (MQ=255)                                            aaTTATCGCGACCGGCACACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGaaact                                                                                                           <  1:211689/251‑3 (MQ=21)                                                             cacCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAGGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttggtatccgtgaggtggcgaaacttcatcagatttctcgtt                                                                                          <  1:275174‑M2/251‑132 (MQ=255)                                                             cacCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttggtatccgtgaggtggcgaaacttcatcagatttctcgtt                                                                                          <  1:215795‑M2/251‑132 (MQ=255)                                                                           gCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttggtatccgtgaggtggcgaaacttcatcagatttctcgttctcaagtcatttat                                                                            <  1:341581‑M2/251‑146 (MQ=255)                                                                           gCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttgg                                                                                                                                <  1:665739‑M2/199‑94 (MQ=255)                                                                               aaaTGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttggtatccgtgaggtggcgaaacttcatcagatttctcgttctcaagtcatttattgga                                                                        >  1:56730‑M2/1‑102 (MQ=255)                                                                                       ggCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttggtatccgtgaggtggcgaaacttcatcagatttctcgttctcaagtcatttattggaaaagcctt                                                                <  1:858432‑M2/251‑158 (MQ=255)                                                                                       ggCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttggtatccgtgaggtggcgaaacttcatcagatttctcgttctcaagtcatttattggaaaagcc                                                                  >  1:130955‑M2/1‑94 (MQ=255)                                                                                       ggCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcacc                                                                                                                                              >  1:464158‑M2/1‑94 (MQ=255)                                                                                             cacaCCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttggtatccgtgaggtggcgaaacttcatcagatttctcgttctcaagtcatttattggaaaagcc                                                                  <  1:68607‑M2/243‑156 (MQ=255)                                                                                              acacCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttggtatccgtgaggtggcgaaacttcatcagatttctcgttctcaagtcatttattggaaaagcctttcgtgaa                                                         >  1:580042‑M2/1‑87 (MQ=255)                                                                                                        cTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttggta                                                                                                                              >  1:728117‑M2/1‑77 (MQ=255)                                                                                                          tCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttggtatccgtgaggtggcgaaacttcatcagatttctcgttctcaagtcatttattggaaaagcctttcgtgaaagagggcttaat                                             <  1:26391‑M2/251‑177 (MQ=255)                                                                                                          tCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAg                                                                                                    >  1:110595/1‑196 (MQ=14)                                                                                                                                   cTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttggtatccgtgaggtggcgaaacttcatcagatttctcgttctcaagtcatttattggaaaagcctttcgtgaaagagggcttaatggcgtaaaatccccttatataaacc                    >  1:427356‑M2/1‑50 (MQ=255)                                                                                                                                   cTAAAAACGCTGCATTTTCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttggtatccgtgaggtggcgaaacttcatcagatttctcgttctcaagtcatttattggaaaagcctttcgtgaaagagggcttaatggcgtaaaatccccttatataaacc                    >  1:254574‑M2/1‑50 (MQ=255)                                                                                                                                                   ttCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttggtatccgtgaggtggcgaaacttcatcagatttctcgttctcaagtcatttattggaaaagcctttcgtgaaagagggcttaatggcgtaaaatccccttatataaaccctcaacgccctaaaat    >  1:828780‑M2/1‑34 (MQ=255)                                                                                                                                                   ttCCGACCGCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttggtatccgtgaggtggcgaaacttcatcagatttctcgttctcaagtcatttattggaaaagcctttcgtgaaagagggcttaatggcgtaaaatccccttatataaaccctcaacgccctaaaat    >  1:185103‑M2/1‑34 (MQ=255)                                                                                                                                                           gCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttggtatccgtgaggtggcgaaacttcatcagatttctcgttctcaagtcatttattggaaaagcctttcgtgaaagagggcttaatggcgtaaaatccccttatat                         <  1:503364‑M2/222‑197 (MQ=255)                                                                                                                                                           gCACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggat                                                                                                                                     <  1:300011‑M2/114‑89 (MQ=255)                                                                                                                                                            cACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttggtatccgtgaggtggcgaaacttcatcagatttctcgttctcaagtcatttattggaaaagcctttcgtgaaagagggcttaatggcgtaaaatccccttatataaaccctcaacgccctaaaatag  <  1:516002‑M2/244‑220 (MQ=255)                                                                                                                                                            cACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttggtatccgtgaggtggcgaaacttcatcagatttctcgttctcaagtcatttattggaaaagcctttcgtgaaagagggcttaatggcgtaaaatccccttatataaaccc                   >  1:9422‑M2/1‑25 (MQ=255)                                                                                                                                                            cACTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttacaaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttggtatcagtgaggtggcgaaacttcatcagatttctcgttctcaagtcatttattggaaaagcctttcgtgcaagagggcttaat                                             <  1:190524‑M2/201‑177 (MQ=255)                                                                                                                                                             aCTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttggtatccgtgaggtggcgaaacttcatcagatttctcgttctcaagtcatttattggaaaagcctttcgtgaaagagggcttaatggcgtaaaatccccttatataaacc                    >  1:829796‑M2/1‑24 (MQ=255)                                                                                                                                                              cTTTGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacaggtcaccgagcaggattttggtatccgtgaggtggcgaaacttcatcagatttctcgttctcaagtcatttattggaaaagcctttcgtgaaagagggc                                                  <  1:571814‑M2/194‑172 (MQ=255)                                                                                                                                                                ttGATCTGCACCCCAAAACCTggacacggacacctaatttgaggtgcagattatgtcttactcttaccaatttcgtttaaaaactatcaaacag                                                                                                                                                    <  1:817441‑M2/94‑74 (MQ=255)                                                                                                                                                                   |                                                                                                                                                                                                                                             GGATTGTAGATCTCGGCCCTGAAGGCGGCAGCGGCGGCGGACAAATTATCGCGACCGGCACACCGGAACAGGTTGCCAAAATGGAGGGCTCGCACACCGCCCGCTTCCTCAAAGAGATTTTGGCAAAAGGCTAAAAACGCTGCATTTTCCGACCGCACTTTGCTCTGCACCCCAAAACCT‑‑‑‑‑‑GGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTATCAATTTCGTTAAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGGAAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAG  >  1/621782‑622174

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