Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	612,643	G→A	intergenic (+52/‑16)	unknown → / → unknown	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	612,643	0	G	A	100.0%	56.5 / NA	31	intergenic (+52/‑16)	unknown/unknown	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (19/12);  total (19/12)

CACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAGATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAGCGTTTAAAAAG  >  1/612404‑612721                                                                                                                                                                                                                                                   |                                                                               cACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAg                                                                            >  1:516002/1‑244 (MQ=255)       gATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGAt                                                               >  1:342208/1‑251 (MQ=255)               ccccAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGTAAAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAaaag                                                        >  1:744744/1‑251 (MQ=255)               ccccAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAg                                                                            >  1:508810/1‑230 (MQ=255)                           gacacggacacCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGa                                                                >  1:629087/1‑230 (MQ=16)                                ggacacCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATc                                         >  1:879860/1‑250 (MQ=11)                                  acacCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATccc                                       <  1:901112/250‑1 (MQ=11)                                     ccTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAATTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGaaca                                 >  1:877292/1‑251 (MQ=9)                                     ccTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAaaagaa                                                      <  1:139465/230‑1 (MQ=9)                                     ccTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGaa                                   >  1:736451/1‑251 (MQ=11)                                      cTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAAACTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGaac                                  <  1:594780/251‑1 (MQ=11)                                        aaTTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGaacaa                                >  1:528141/1‑251 (MQ=11)                                        aaTTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGAGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAaaagaaag                                                    <  1:469041/229‑1 (MQ=9)                                         aTTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGaacaaaca                             >  1:829844/1‑251 (MQ=14)                                         aTTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGaacaaaca                             >  1:851005/1‑251 (MQ=14)                                                     gATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGGGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCa                 >  1:250326/1‑251 (MQ=9)                                                     gATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGaacaaact                             <  1:380148/241‑2 (MQ=11)                                                      aTTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTccc                   <  1:223326/248‑1 (MQ=14)                                                         aTGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAg             <  1:207233/251‑1 (MQ=14)                                                         aTGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAaaag                                                        <  1:404710/210‑1 (MQ=11)                                                           gTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGAt                                                               <  1:320903/199‑1 (MQ=12)                                                           gTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAGCg           >  1:450059/1‑251 (MQ=14)                                                           gTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAg             >  1:897128/1‑251 (MQ=11)                                                                actcttacCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAGCGTTTaa      >  1:366174/1‑251 (MQ=14)                                                                   cttacCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGAc                          <  1:656815/228‑1 (MQ=14)                                                                   cttacCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAGCGTTTaaaaa   >  1:732945/1‑251 (MQ=14)                                                                   cttacCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTCCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAGCGTTTaaaaa   >  1:297610/1‑251 (MQ=14)                                                                   cttacCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAg             >  1:249516/1‑241 (MQ=14)                                                                   cttacCAATTGCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATAAAATCTCCCATTCGCTCGTCATTTATTGGCTAAAAGCATTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAg                                                                            <  1:454698/180‑1 (MQ=11)                                                                   cttacCAACTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGa                                                                         <  1:779677/181‑1 (MQ=17)                                                                     tacCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCGAAACTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAAATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAGCGTTTAAAAAg  >  1:80582/1‑250 (MQ=14)                                                                                                                                                                                                                                                   |                                                                               CACTTTGATCTGCACCCCAAAACCTGGACACGGACACCTAATTTGAGGTGCAGATTATGTCTTACTCTTACCAATTTCGTTTAAAAACTATCAAACAGGTCACCGAGCAGGATTTTGGTATCCGTGAGGTGGCTAAATTTCATCAGATTTCTCGTTCTCAAGTCATTTATTGG‑‑AAAAGCCTTTCGTGAAAGAGGGCTTAATGGCGT‑AAAATCCCCTTATATAAACCCTCAACGCCCTAAGATAGTGAAGCCAAAGATGAAAAAGAAAGCGATTGAAATCCCGGAACAAACAGACTTTTCCCCAAAAGCGTTTAAAAAG  >  1/612404‑612721

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