Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,285,403	G→C	G42G (GGC→GGG)	unknown ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,285,403	0	G	C	100.0%	101.2 / NA	37	G42G (GGC→GGG)	unknown	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (21/16);  total (21/16)

GCTAATTACGTAAGAGCGGTAATTATATAAAACGAATTTAAATAAACAGAAAATTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGATATAATGCAATGATTTATAAATGAATAAAAATGATATTGTTTATTCTCTACCGTCATTTTTACCGTCACCAGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGGCCCATAGACGCTTTAACACTTTTGGGGAATTTACCTTGTTTAATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAATATTTCGCGGATCTTGTCGCCTGATTTCTGCTCAATGAAATCATTTATCACTTTATA  >  1/1285168‑1285651                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                         gCTAATTACGTAAGAGCGGTAATTATATAAAACGAATTTAAATAAACAGAAAATTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTa                                                                                                                                                                                                                                           >  1:267202/1‑251 (MQ=255)                                                    aaTTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCg                                                                                                                                                                                        >  1:353600/1‑251 (MQ=255)                                                                    gAAGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACtt                                                                                                                                                                        >  1:534059/1‑251 (MQ=255)                                                                     aaGTATCCCACTCGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTcatcat                                                                                                                                                                                                           >  1:573166/1‑215 (MQ=255)                                                                                 cGATTATTTATCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTg                                                                                                                                                            >  1:811968/1‑250 (MQ=255)                                                                                      atttatCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAg                                                                                                                                                                                  <  1:186895/223‑1 (MQ=255)                                                                                      atttatCGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTa                                                                                                                                                       >  1:113098/1‑250 (MQ=255)                                                                                            cGGCGACATAATGCAATGATTTATCAATGAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTtagta                                                                                                                                                >  1:787709/1‑251 (MQ=255)                                                                                                                        gAATAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGtctcatca                                                                                                                                                                                                            <  1:558209/163‑8 (MQ=255)                                                                                                                         aaTAAAATTGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTtctc                                                                                                                   <  1:791741/251‑1 (MQ=255)                                                                                                                                ttGATATTGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAAttt                                                                                                                                                                                                                          >  1:787692/2‑141 (MQ=255)                                                                                                                                       tGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACGCGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAAtt                                                                                                                                                                                                 <  1:844659/159‑1 (MQ=255)                                                                                                                                       tGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTAcc                                                                                                                                                                                                                       >  1:374983/1‑137 (MQ=255)                                                                                                                                       tGTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCtt                                                                                                                                                         <  1:78960/199‑1 (MQ=255)                                                                                                                                        gTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTtagtag                                                                                                                                               <  1:208337/208‑1 (MQ=255)                                                                                                                                        gTTTATTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAg                                                                                                    >  1:825922/1‑251 (MQ=255)                                                                                                                                            aTTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGTTTTTATCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCa                                                                                                >  1:412077/1‑251 (MQ=255)                                                                                                                                            aTTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCtt                                                                                                                                                         <  1:730566/194‑1 (MQ=255)                                                                                                                                            aTTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCa                                                                                                >  1:340052/1‑251 (MQ=255)                                                                                                                                            aTTCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCAATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCt                                                                                                                                                                                >  1:361858/1‑171 (MQ=255)                                                                                                                                             ttCTATACCGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCAt                                                                                               >  1:488936/1‑251 (MQ=255)                                                                                                                                                    ccGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTgggg                                                                                                                                                                                                                               <  1:468394/116‑1 (MQ=255)                                                                                                                                                     cGTCATTTTTACCGTCTTTGGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGc                                                                                       <  1:412717/251‑1 (MQ=255)                                                                                                                                                          tttttaccgtctttggaagaaAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGtt                                                                                  >  1:142387/16‑251 (MQ=255)                                                                                                                                                          tttttaccgtctttggaagaaAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGtt                                                                                  >  1:640421/16‑251 (MQ=255)                                                                                                                                                                   tctttggaagaaAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTGCCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTGCTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCa                                                                         <  1:249440/245‑1 (MQ=255)                                                                                                                                                                                    cATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAatat                                                        >  1:336520/1‑251 (MQ=255)                                                                                                                                                                                     aTCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACtt                                                                                                                                                    <  1:884412/158‑1 (MQ=255)                                                                                                                                                                                     aTCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAATAtt                                                       >  1:287794/1‑251 (MQ=255)                                                                                                                                                                                       cAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAAttgttg                                                                                                                                                                                             >  1:683844/1‑115 (MQ=255)                                                                                                                                                                                          gAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAATAttt                                                      <  1:15111/247‑1 (MQ=255)                                                                                                                                                                                                     aaGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAATATTTCGCGGATCTTGTCGc                                       <  1:417985/251‑1 (MQ=255)                                                                                                                                                                                                       gCCAAGTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAAt                                                                                                                                                                                                  <  1:846544/94‑1 (MQ=255)                                                                                                                                                                                                            gTATCAATATCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTc                                                                                            <  1:397612/191‑1 (MQ=255)                                                                                                                                                                                                                  atatCAGATTTTCTCCACGCTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCACGATTTAATATTTCGCGGATCTTGTCGCCTGATTTCTGCt                           >  1:246589/1‑250 (MQ=255)                                                                                                                                                                                                                                     cTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAATATTTCGCGGATCTTGTCGCCTGATTTCTGCTc                          <  1:671454/232‑1 (MQ=255)                                                                                                                                                                                                                                     cTACACGCCCCATGGACGCTTTAACACTTTTGGGGAATTTACCTTGTCTCATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTG‑‑‑‑‑CTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAATATTTCGCGGATCTTGTCGCCTGATTTCTGCTCAATGAAATCATTTATCACTTtata  >  1:68068/1‑251 (MQ=37)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                         GCTAATTACGTAAGAGCGGTAATTATATAAAACGAATTTAAATAAACAGAAAATTATATGAAATGGGGAAGTATCCCACTCGATTATTTATCGGCGATATAATGCAATGATTTATAAATGAATAAAAATGATATTGTTTATTCTCTACCGTCATTTTTACCGTCACCAGAAGAAAATGCCATCAAGAACGTTCAGCAAGCCAAGTATCAATATCAGATTTTCTCCACGCTACACGGCCCATAGACGCTTTAACACTTTTGGGGAATTTACCTTGTTTAATCATACGCCAAATTGTTGTGTCGCTCAAGCTGGTGACTTCAATCACCTTTGCTTTACTTAGTAGCGCTGGTTGTTGTTCCATTATTCTTCTCCATTAAAACAAAAAGCCCATTTCTGGGCTTGTTATCTAACCAATATCGCGATTTAATATTTCGCGGATCTTGTCGCCTGATTTCTGCTCAATGAAATCATTTATCACTTTATA  >  1/1285168‑1285651

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