Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	113	4,960	C→A	intergenic (–/–)	– / –	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	113	4,960	0	C	A	100.0%	175.1 / NA	50	intergenic (–/–)	–/–	–/–
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (20/30);  total (20/30)

GACCTGGGCAGTAATGGTGCTTAAGCGCCTGCTGCCGGGGCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCCATGTCAAAGAAAATATCGACAACAGAGAGCGACTCCATGTCCGTAGAAAATATTGTCAACATTAACGAATCTAACCTGCAACAGGTTCTTGAACAGT  >  113/4867‑5057                                                                                              |                                                                                                  gACCTGGGCAGTAATGGTGCTTAAGCGCCTGCTGCCGGGGCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgta                                                                                            >  1:789484‑M1/1‑98 (MQ=255)  aCCTGGGCAGTAATGGTGCTTAAGCGCCTGCTGCCGGGGCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtac                                                                                           >  2:1111436‑M1/1‑97 (MQ=255)   ccTGGGCAGTAATGGTGCTTAAGCGCCTGCTGCCGGGGCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtact                                                                                          <  1:1668895‑M1/101‑6 (MQ=255)    cTGGGCAGTAATGGTGCTTAAGCGCCTGCTGCCGGGGCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactg                                                                                         >  1:325457‑M1/1‑95 (MQ=255)      gggCAGTAATGGTGCTTAAGCGCCTGCTGCCGGGGCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgac                                                                                       <  2:1631737‑M1/101‑9 (MQ=255)        gCAGTAATGGTGCTTAAGCGCCTGCTGCCGGGGCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTGAGCCCCAATGTgtactgaccc                                                                                     >  2:1340927‑M1/1‑91 (MQ=255)        gCAGTAATGGTGCTTAAGCGCCTGCTGCCGGGGCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccc                                                                                     >  1:888799‑M1/1‑91 (MQ=255)         cAGTAATGGTGCTTAAGCGCCTGCTGCCGGGGCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgacccc                                                                                    >  2:542547‑M1/1‑90 (MQ=255)                    cTTAAGCGCCTGCTGCCGGGGCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttgga                                                                         >  2:1357830‑M1/1‑79 (MQ=255)                     ttAAGCGCCTGCTGCCGGGGCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggac                                                                        <  1:1713054‑M1/101‑24 (MQ=255)                     ttAAGCGCCTGCTGCCGGGGCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggac                                                                        <  1:970327‑M1/101‑24 (MQ=255)                     ttAAGCGCCTGCTGCCGGGGCGCGAGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggac                                                                        >  2:176257‑M1/1‑78 (MQ=255)                       aaGCGCCTGCTGCCGGGGCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacag                                                                      >  2:615514‑M1/1‑76 (MQ=255)                         gcgcCTGCTGCCGGGGCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagtt                                                                    <  1:293842‑M1/101‑28 (MQ=255)                             ctgctgCCGGGGCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaac                                                                >  2:1190890‑M1/1‑70 (MQ=255)                             ctgctgCCGGGGCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaac                                                                >  1:427281‑M1/1‑70 (MQ=255)                             ctgctgCCGGGGCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaac                                                                <  1:1280361‑M1/101‑32 (MQ=255)                              tgctgcCGGGGCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaaca                                                               >  2:437171‑M1/1‑69 (MQ=255)                              tgctgcCGGGGCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaaca                                                               <  1:1340933‑M1/101‑33 (MQ=255)                                     ggggCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggca                                                        <  2:336621‑M1/101‑40 (MQ=255)                                     ggggCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGACCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggca                                                        >  2:521211‑M1/1‑62 (MQ=255)                                         cgcgTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatata                                                    >  1:1586828‑M1/1‑58 (MQ=255)                                          gcgTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatatag                                                   <  2:1155007‑M1/101‑45 (MQ=255)                                           cgTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatatagg                                                  >  1:357090‑M1/1‑56 (MQ=255)                                             tGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtc                                                <  2:715476‑M1/101‑48 (MQ=255)                                             tGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggcc                                                >  1:1334847‑M1/1‑54 (MQ=255)                                              gATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtct                                               <  1:1064732‑M1/101‑49 (MQ=255)                                              gATGGACAAAATATTGCAGGGGAGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggact                                               >  2:1284822‑M1/1‑53 (MQ=255)                                               aTGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctg                                              <  2:1483896‑M1/101‑50 (MQ=255)                                               aTGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctg                                              <  1:209993‑M1/101‑50 (MQ=255)                                               aTGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctg                                              <  1:1371494‑M1/101‑50 (MQ=255)                                                  gACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgatt                                           <  2:1143959‑M1/101‑53 (MQ=255)                                                   aCAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattc                                          >  2:1701722‑M1/1‑48 (MQ=255)                                                    cAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattcc                                         <  1:615519‑M1/101‑55 (MQ=255)                                                         tatTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatat                                    <  1:583970‑M1/101‑60 (MQ=255)                                                         tatTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatat                                    <  2:888794‑M1/101‑60 (MQ=255)                                                               aGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaatt                              <  1:1223616‑M1/101‑66 (MQ=255)                                                                  ggTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattgga                           <  1:1280305‑M1/101‑69 (MQ=255)                                                                   gTGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggag                          <  1:1515870‑M1/101‑70 (MQ=255)                                                                    tGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagt                         <  1:1199135‑M1/101‑71 (MQ=255)                                                                    tGAGTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagt                         >  1:668389‑M1/1‑31 (MQ=255)                                                                       gTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcag                      >  1:980265‑M1/1‑28 (MQ=255)                                                                       gTTGAAGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcag                      <  2:585817‑M1/101‑74 (MQ=255)                                                                            aGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagacctt                 <  2:244260‑M1/101‑79 (MQ=255)                                                                            aGCGCGCGCTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagacctt                 <  1:1190895‑M1/101‑79 (MQ=255)                                                                             gcgcgcgcTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttt                <  2:1082481‑M1/101‑80 (MQ=255)                                                                              cgcgcgcTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagacctttt               <  2:333898‑M1/101‑81 (MQ=255)                                                                              cgcgcgcTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagacctttt               <  1:599157‑M1/101‑81 (MQ=255)                                                                               gcgcgcTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagacctttta              <  1:986567‑M1/101‑82 (MQ=255)                                                                                  cgcTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttttaatt           >  2:694781‑M1/1‑17 (MQ=255)                                                                                    cTTAAGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttttaatttc         <  1:1206873‑M1/101‑87 (MQ=255)                                                                                       aaGCCCCAATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttttaatttcagg      <  1:101563‑M1/101‑90 (MQ=255)                                                                                          ccccaATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttttaatttcaggcta   >  2:1303951‑M1/1‑9 (MQ=255)                                                                                          ccccaATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttttaatttcaggcta   >  2:41168‑M1/1‑9 (MQ=255)                                                                                           cccaATGTgtactgaccccaaaaagttggacagttaaacacgaggcatataggtctgattccgatattcaattggagtcagaccttttaatttcaggctaa  >  2:953805‑M1/1‑8 (MQ=255)                                                                                              |                                                                                                  GACCTGGGCAGTAATGGTGCTTAAGCGCCTGCTGCCGGGGCGCGTGATGGACAAAATATTGCAGGGGTGAGTTGAAGCGCGCGCTTAAGCCCCCATGTCAAAGAAAATATCGACAACAGAGAGCGACTCCATGTCCGTAGAAAATATTGTCAACATTAACGAATCTAACCTGCAACAGGTTCTTGAACAGT  >  113/4867‑5057

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