Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,481,633	G→A	intergenic (+284/‑725)	unknown → / → unknown	Zona occludens toxin‑like phage protein/Bacterial non‑heme ferritin (EC 1.16.3.2)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,481,633	0	G	A	100.0%	127.3 / NA	38	intergenic (+284/‑725)	unknown/unknown	Zona occludens toxin‑like phage protein/Bacterial non‑heme ferritin (EC 1.16.3.2)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (15/23);  total (15/23)

GACTGGCCGAAAACAAGATAACTGAAGGCGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTA‑TTTTCCGCTTGC‑GG‑AAAATGCAGGGCGAAGGCTTGCAGCCCGGTAAATAATGCTCGATGTTATTCGGGCATTTTTATTTGTTGCGCG‑CTGATGATTTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTTCATAACGTTTAACT  >  1/1481407‑1481815                                                                                                                                                                                                                                       |                                                                                                                                                                                       gtcTGGCTGAAAATAAAATAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAAt                                                                                                                                                                  >  1:850041/3‑251 (MQ=255)           aaataaaaTAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGccc                                                                                                                                                         >  1:712035/8‑250 (MQ=255)                aaaTAACGGAAGGTGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCaa                                                                                                                                                   <  1:689263/249‑1 (MQ=255)                             tgAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTAA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGc                                                                                                                                      <  1:484113/250‑1 (MQ=255)                             tgAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTGGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGc                                                                                                                                      <  1:575145/250‑1 (MQ=255)                                    ggTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGtaa                                                                                                                               >  1:875811/1‑251 (MQ=255)                                        tttCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGtaataa                                                                                                                            >  1:521567/1‑250 (MQ=255)                                        tttCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAAc                                                                                                                           <  1:41444/251‑1 (MQ=255)                                           cATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAActgc                                                                                                                        <  1:498703/251‑1 (MQ=255)                                                    cTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGAt                                                                                                               <  1:790901/251‑1 (MQ=255)                                                            gtttgatagttttttaacgaaattGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGAt                                                                                                       >  1:297496/24‑251 (MQ=255)                                                                       aaaaTCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGaa                                                                                            <  1:341661/251‑1 (MQ=255)                                                                                   tGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAAc                                                                                                                           >  1:400230/1‑208 (MQ=255)                                                                                      tACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGc                                                                                                          <  1:531319/222‑1 (MQ=255)                                                                                      tACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATAc                                                                                                     >  1:303813/1‑227 (MQ=255)                                                                                      tACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACg                                                                                                    <  1:613501/228‑1 (MQ=255)                                                                                       aCTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACAtttt                                                                            >  1:739068/1‑251 (MQ=255)                                                                                        cTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACAttttt                                                                           >  1:756241/1‑251 (MQ=255)                                                                                            gACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTaa                                                                                                                                         <  1:723725/185‑1 (MQ=255)                                                                                            gACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTa                                                                                        >  1:791023/1‑234 (MQ=255)                                                                                              cAGATATAGTGATTTTCCTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCAGGCAGAATGTTATTCGGGCATTTTTATTTTTTCCGCG‑TCCATAAATTAAATCGAGAATTGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATAc                                                                                                     <  1:226758/219‑1 (MQ=255)                                                                                                 atatAGTGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGAt                                                                                                               <  1:171824/206‑1 (MQ=255)                                                                                                       tGATTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGGTTGCAGCCCGATAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCGTCCCA‑AAATTTAATCAAGAAATAAAAATTTACCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATa                                                            <  1:711642/251‑1 (MQ=255)                                                                                                         aTTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAAc                                                                                                                           >  1:799159/1‑186 (MQ=255)                                                                                                         aTTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTc                                                                  <  1:426022/243‑1 (MQ=255)                                                                                                         aTTTTACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATaaa                                                          <  1:268366/251‑1 (MQ=255)                                                                                                          ttttACTCCAAGTGAAACTTGTCTAACTATTTTTGAGCAAGCTCA‑AAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAAt                                                         >  1:784907/1‑251 (MQ=255)                                                                                                                   aaGTGAAACTTGTCTAACTATTTTTGAGCAAGC‑TCAAAAATGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAAtattta                                                >  1:844796/1‑251 (MQ=255)                                                                                                                                   aactatttttgggcaagctcaaaaaTGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAgg                                  >  1:346280/22‑250 (MQ=255)                                                                                                                                   aactatttttgagcaagctcaaaaaTGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGc                                 <  1:236940/230‑1 (MQ=255)                                                                                                                                                   gctcaaaaaTGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATa                                                                                                                                                <  1:720040/119‑1 (MQ=255)                                                                                                                                                      caaaaaTGTAGGGCGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAg                                                                                                                                       <  1:233728/128‑1 (MQ=255)                                                                                                                                                                   cGAAGGCTTGCAGCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAAcc  <  1:367120/250‑2 (MQ=255)                                                                                                                                                                               gCCCGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACaa                                    <  1:828074/204‑1 (MQ=255)                                                                                                                                                                                cccGGTAAACCATGCTCGATGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAatat                                                   <  1:530159/188‑1 (MQ=255)                                                                                                                                                                                                  aTGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGc                                                                                                                                      >  1:791278/1‑87 (MQ=255)                                                                                                                                                                                                  aTGTTATTCGGGCATTTTTATTTTTTGCGCG‑TACCTAAATTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGc                                                                                     <  1:51778/136‑1 (MQ=255)                                                                                                                                                                                                  aTGTTATTCGGGCATTTTTATTTTTTGCGAA‑TACCTAAATTTAACCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGc                                                                                     <  1:785375/136‑1 (MQ=255)                                                                                                                                                                                                                                       |                                                                                                                                                                                       GACTGGCCGAAAACAAGATAACTGAAGGCGAAACAGGTTTTTCATTTCTTCCTAAAGGTTTTTTAACGGGAAAATCCAAAAATGATACTAAGACAGATATAGTGATTTTACTCCAAGTGAAACTTGTCTAACTA‑TTTTCCGCTTGC‑GG‑AAAATGCAGGGCGAAGGCTTGCAGCCCGGTAAATAATGCTCGATGTTATTCGGGCATTTTTATTTGTTGCGCG‑CTGATGATTTTAATCGAGAAATGAAAATTTCCCGCCCCTGCAAATACCCCTAAAGCCACGTAATAACTGCTGATTAGATAATGCGATACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTTCATAACGTTTAACT  >  1/1481407‑1481815

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