Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	AM260479	518,612	G→C	intergenic (+17/‑48)	h16_A0492 → / → pgam2	Phospholipid‑binding protein/phosphoglycerate mutase 2 protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260479	518,612	0	G	C	100.0%	154.9 / NA	40	intergenic (+17/‑48)	h16_A0492/pgam2	Phospholipid‑binding protein/phosphoglycerate mutase 2 protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (19/21);  total (19/21)

CGGCGCGCAGGTGCGCGAGGCGATCCAGGGCCATGTGCTGGCGCAGGCCAGCCTCACCGGCACTTATACGCTGAACCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGGGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGCCCCGGGCCGCATTCGCTGGCCTATACCCACCTGATCGTGATCCGCCACGGC  >  AM260479/518502‑518722                                                                                                               |                                                                                                               cGGCGCGCAGGTGCGCGAGGCGATCCAGGGCCATGTGCTGGCGCAGGCCAGCCTCACCGGCACTTATACGCTGAACCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCt                                                                                          >  1:931196/1‑133 (MQ=255)         aGGTGCGCGAGGCGATCCAGGGCCATGTGCTGGCGCAGGCCAGCCTCACCGGCACTTATACGCTGAACCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGaaa                                                                                 >  1:423607/1‑134 (MQ=255)                cgAGGCGATCCAGGGCCATGTGCTGGCGCAGGCCAGCCTCACCGGCACTTATACGCTGAACCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGaa                                                                          <  2:518036/134‑1 (MQ=255)                    gCGATCCAGGGCCATGTGCTGGCGCAGGCCAGCCTCACCGGCACTTATACGCTGAACCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGa                                                                       >  1:310274/1‑133 (MQ=255)                             ggCCATGTGCTGGCGCAGGCCAGCCTCACCGGCACTTATACGCTGAACCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGt                                                             <  1:1175907/134‑1 (MQ=255)                                      cTGGCGCAGGCCAGCCTCACCGGCACTTATACGCTGAACCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAAc                                                                    <  1:118787/118‑1 (MQ=255)                                      cTGGCGCAGGCCAGCCTCACCGGCACTTATACGCTGAACCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAAc                                                                    >  2:118787/1‑118 (MQ=255)                                      cTGGCGCAGGCCAGCCTCACCGGCACTTATACGCTGAACCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGcc                                                    >  2:729205/1‑134 (MQ=255)                                            cAGGCCAGCCTCACCGGCACTTATACGCTGAACCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAATGGCATGTCGCAAAGCCCCGGGc                                              <  2:1050477/134‑1 (MQ=255)                                             aGGCCAGCCTCACCGGCACTTATACGCTGAACCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCcgcgc                                                                                                     <  1:774811/78‑1 (MQ=255)                                             aGGCCAGCCTCACCGGCACTTATACGCTGAACCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCcgcgc                                                                                                     >  2:774811/1‑78 (MQ=255)                                                ccAGCCTCACCGGCACTTATACGCTGAACCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGCCCCGGGCCGCa                                          <  2:435668/134‑1 (MQ=255)                                                  aGCCTCACCGGCACTTATACGCTGAACCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGCCCCGGGCCGCAtt                                        <  2:310274/134‑1 (MQ=255)                                                  aGCCTCACCGGCACTTATACGCTGAACCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGCCCCGGGCCGCAtt                                        >  1:610604/1‑134 (MQ=255)                                                        aCCGGCACTTATACGCTGAACCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACg                                                                   <  1:214961/101‑1 (MQ=255)                                                        aCCGGCACTTATACGCTGAACCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACg                                                                   >  2:214961/1‑101 (MQ=255)                                                          cGGCACTTATACGCTGAACCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGCCCCGGGCCGCATTCGCTGGcc                                <  2:931196/134‑1 (MQ=255)                                                                      gCTGAACCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATc                                                                               <  1:1126444/75‑1 (MQ=255)                                                                      gCTGAACCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATc                                                                               >  2:1126444/1‑75 (MQ=255)                                                                           aCCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGCCCCGGGCCGCATTCGc                                     >  1:907844/1‑112 (MQ=255)                                                                           aCCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGCCCCGGGCCGCATTCGc                                     <  2:907844/112‑1 (MQ=255)                                                                              ccAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAAc                                                                    >  1:296837/1‑78 (MQ=255)                                                                              ccAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAAc                                                                    <  2:296837/78‑1 (MQ=255)                                                                              ccAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGCCCCGGGCCGCa                                          >  1:106425/1‑104 (MQ=255)                                                                              ccAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGCCCCGGGCCGCa                                          <  1:260605/104‑1 (MQ=255)                                                                              ccAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGCCCCGGGCCGCa                                          >  2:260605/1‑104 (MQ=255)                                                                              ccAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGCCCCGGGCCGCa                                          <  2:106425/104‑1 (MQ=255)                                                                              ccAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGCCCCGGGCCGCATTCGc                                     <  2:1080425/109‑1 (MQ=255)                                                                              ccAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGCCCCGGGCCGCATTCGc                                     >  1:1080425/1‑109 (MQ=255)                                                                              ccAATCTGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGCCCCGGGCCGCATTCGCTGGCCTATACCCACCTGATCGTGAt            >  1:563532/1‑134 (MQ=255)                                                                                    tGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCgcgg                                                                                                   >  1:384063/1‑41 (MQ=255)                                                                                    tGGCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCgcgg                                                                                                   <  2:384063/41‑1 (MQ=255)                                                                                      gCGAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGCCCCGGGCCGCATTCGCTGGCCTATACCCACCTGATCGTGATCCGc        <  2:563532/130‑1 (MQ=255)                                                                                        gAAATAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGCCCCGGGCCGCATTCGCTGGCCTATACCCACCTGATCGTGATCCGCCACGGc  <  1:729205/134‑1 (MQ=255)                                                                                          aaTAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGCCCCGGGCCGCATTCGCTGGCCTATACCCACCTGATCg                >  1:676386/1‑118 (MQ=255)                                                                                          aaTAACGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGCCCCGGGCCGCATTCGCTGGCCTATACCCACCTGATCg                <  2:676386/118‑1 (MQ=255)                                                                                               cGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGCCCCGGGCCGCATTCGCTGGCCTATACCCACCTGATCGTGATCCGCCACGGc  <  2:850017/127‑1 (MQ=255)                                                                                               cGTCTGGCCGGCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGCCCCGGGCCGCATTCGCTGGCCTATACCCACCTGATCGTGATCCGCCACGGc  >  1:850017/1‑127 (MQ=255)                                                                                                        ggCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAg                                                                         <  2:83117/47‑1 (MQ=255)                                                                                                        ggCGCCGCGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAg                                                                         >  1:83117/1‑47 (MQ=255)                                                                                                               |                                                                                                               CGGCGCGCAGGTGCGCGAGGCGATCCAGGGCCATGTGCTGGCGCAGGCCAGCCTCACCGGCACTTATACGCTGAACCCCAATCTGGCGAAATAACGTCTGGCCGGCGCCGGGCCAGCCGCGCGGCGCCGGCCTCAAAGGAAATCCCGAAGGAAACGGCATGTCGCAAAGCCCCGGGCCGCATTCGCTGGCCTATACCCACCTGATCGTGATCCGCCACGGC  >  AM260479/518502‑518722

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