Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	58	3,930	G→A	intergenic (–/–)	– / –	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	58	3,930	0	G	A	100.0%	242.7 / NA	69	intergenic (–/–)	–/–	–/–
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (33/36);  total (33/36)

ACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTA  >  58/3834‑4024                                                                                                 |                                                                                               acGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGac                                                                                            <  2:417808/101‑1 (MQ=255)  cGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGaca                                                                                           <  2:20493/101‑1 (MQ=255)  cGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGaca                                                                                           <  1:1300889/101‑1 (MQ=255)   gCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGacac                                                                                          >  2:406248/1‑101 (MQ=255)    cTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACAcc                                                                                         <  2:11015/101‑1 (MQ=255)     tGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCa                                                                                        >  1:1278745/1‑101 (MQ=255)          ggCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGtt                                                                                   >  1:1422167/1‑101 (MQ=255)             aaCGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCg                                                                                >  2:1467801/1‑101 (MQ=255)               cGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTg                                                                              <  2:156174/101‑1 (MQ=255)                gTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTgg                                                                             <  2:874143/101‑1 (MQ=255)                gTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCATTTAAATGCACAGAAGACACCAGCGTTGCgggg                                                                             <  1:169005/101‑1 (MQ=255)                  aGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGt                                                                           >  1:1475828/1‑101 (MQ=255)                  aGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGt                                                                           >  2:948282/1‑101 (MQ=255)                  aGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGt                                                                           >  1:1488688/1‑101 (MQ=255)                    ccccAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGtat                                                                         <  2:1337953/101‑1 (MQ=255)                    ccccAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGtat                                                                         <  1:1387795/101‑1 (MQ=255)                     cccAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGtata                                                                        >  1:405907/1‑101 (MQ=255)                     cccAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGtata                                                                        >  2:140327/1‑101 (MQ=255)                        aGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGc                                                                     <  1:1455802/101‑1 (MQ=255)                           gcATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATg                                                                  >  1:30282/1‑101 (MQ=255)                            cATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGt                                                                 >  1:737078/1‑101 (MQ=255)                             aTGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTa                                                                >  2:681711/1‑101 (MQ=255)                              tGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTaa                                                               >  1:1168161/1‑101 (MQ=255)                                 tCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAAccc                                                            <  2:209886/101‑1 (MQ=255)                                  ccTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCt                                                           <  1:1450254/101‑1 (MQ=255)                                  ccTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCt                                                           >  2:415414/1‑101 (MQ=255)                                  ccTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCt                                                           >  1:219847/1‑101 (MQ=255)                                   cTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCtt                                                          <  1:1048114/101‑1 (MQ=255)                                   cTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCtt                                                          >  2:914488/1‑101 (MQ=255)                                   cTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCtt                                                          >  2:63789/1‑101 (MQ=255)                                    tCATCGTCAAATGAATAGCCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCttt                                                         <  1:676852/101‑1 (MQ=255)                                    tCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCttt                                                         <  2:1464733/101‑1 (MQ=255)                                     cATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTg                                                        >  2:1230226/1‑101 (MQ=255)                                      aTCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGc                                                       >  2:53745/1‑101 (MQ=255)                                      aTCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGc                                                       <  1:165875/101‑1 (MQ=255)                                       tCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCg                                                      <  1:373243/101‑1 (MQ=255)                                         gTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGt                                                    >  1:1091774/1‑101 (MQ=255)                                          tCAAATGAATATCCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTg                                                   <  1:1467807/101‑1 (MQ=255)                                          tCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTg                                                   <  1:598145/101‑1 (MQ=255)                                           cAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTgc                                                  >  2:744906/1‑101 (MQ=255)                                             aaTGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTgcag                                                <  2:1042194/101‑1 (MQ=255)                                             aaTGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTgcag                                                <  1:1498586/101‑1 (MQ=255)                                                    atccccGTTTGCGCGTCTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTc                                         <  2:1312772/99‑1 (MQ=255)                                                    aaCCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTc                                         >  1:499054/1‑101 (MQ=255)                                                    aaCCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTc                                         <  1:1280518/101‑1 (MQ=255)                                                     aCCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCg                                        >  1:1135956/1‑101 (MQ=255)                                                       cccGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTc                                      >  1:594204/1‑101 (MQ=255)                                                             tGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAACGTTGCGTGGGTATAGCCATGTAACCCTTTGCGGTGCAGAAACTTCGTCACCTGt                                >  2:1691275/1‑101 (MQ=255)                                                                gcgTTTGGGCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGttct                             <  1:954570/101‑1 (MQ=255)                                                                      gggCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTc                       <  1:1469604/101‑1 (MQ=255)                                                                       ggCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCCCCTGTGCTTCGCTCa                      <  2:54323/101‑1 (MQ=255)                                                                       ggCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCa                      <  1:1292211/101‑1 (MQ=255)                                                                       ggCGAGATCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCa                      <  2:680789/101‑1 (MQ=255)                                                                          gagaTCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGtt                   >  1:389547/1‑101 (MQ=255)                                                                           agaTCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTg                  >  2:205498/1‑101 (MQ=255)                                                                            gaTCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTgg                 >  2:53035/1‑101 (MQ=255)                                                                            gaTCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGAAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTgg                 <  1:1278011/101‑1 (MQ=255)                                                                             aTCTTCTTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATTGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTggg                <  2:1422161/101‑1 (MQ=255)                                                                              tcttctTTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGGGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGc               >  1:1027903/1‑101 (MQ=255)                                                                                ttcttTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTa             <  1:1389159/101‑1 (MQ=255)                                                                                ttcttTTAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTa             <  1:780736/101‑1 (MQ=255)                                                                                    tttAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTaaaaa         >  1:1428806/1‑101 (MQ=255)                                                                                     ttAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAg        >  1:1592625/1‑101 (MQ=255)                                                                                      tAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAgg       <  2:338629/101‑1 (MQ=255)                                                                                      tAAATGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAgg       <  2:1475822/101‑1 (MQ=255)                                                                                       aaaTGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGc      >  2:368412/1‑101 (MQ=255)                                                                                         aTGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCtt    <  1:407433/101‑1 (MQ=255)                                                                                         aTGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCtt    >  1:1652393/1‑101 (MQ=255)                                                                                         aTGCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCtt    >  1:958360/1‑101 (MQ=255)                                                                                           gCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCCCCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTa  <  2:1428589/101‑1 (MQ=255)                                                                                           gCACAGAAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTa  <  1:140327/101‑1 (MQ=255)                                                                                                 |                                                                                               ACGCTGCAAGGCAACGTAGCCCCAGCGCATGTTCCTCATCGTCAAATGAATAACCCCGTTTGCGCGTTTGGGCGAGATCTTCTTTTAAATGCACAGGAGACACCAGCGTTGCGTGGGTATAGGCATGTAACCCTTTGCGGTGCAGCAGCTTCGTCACCTGTTCTTCGCTCAGTTGGGCTAAAAAGGCTTTA  >  58/3834‑4024

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