Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	193	2,088	A→C	intergenic (–/–)	– / –	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	193	2,088	0	A	C	100.0%	212.8 / NA	59	intergenic (–/–)	–/–	–/–
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (27/32);  total (27/32)

CAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGCCAGCGTGATGTTGCTGCACTGGATGAAAAATACACGAAGGAGTTAGCTAATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAG  >  193/1990‑2183                                                                                                   |                                                                                                cAAAGAGCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGcaa                                                                                               <  2:859550/101‑1 (MQ=255)  aaagagcAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGcaaa                                                                                              >  1:1415083/7‑101 (MQ=255)   aagagcAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGcaaaa                                                                                             >  1:1333600/6‑101 (MQ=255)    agagcAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAAt                                                                                            <  1:56709/97‑1 (MQ=255)      agcAGCGCGACAAAAATGCCAGAGAACTGAAACTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATac                                                                                          >  1:652392/3‑101 (MQ=255)        cAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATacac                                                                                        <  2:88839/101‑1 (MQ=255)         aGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACg                                                                                       >  2:1012983/1‑101 (MQ=255)         aGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACg                                                                                       <  2:1067594/101‑1 (MQ=255)          gcgcGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGa                                                                                      <  1:82877/101‑1 (MQ=255)          gcgcGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGa                                                                                      <  1:1131324/101‑1 (MQ=255)           cgcgACAAAAATTCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGaa                                                                                     <  2:1240601/101‑1 (MQ=255)           cgcgACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGaa                                                                                     <  2:1453681/101‑1 (MQ=255)            gcgACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAg                                                                                    <  2:632572/101‑1 (MQ=255)            gcgACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAg                                                                                    >  1:613868/1‑101 (MQ=255)             cgACAAAAATGCCAGAGAACTGAAGCTGGCGAACGGGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAgg                                                                                   >  1:1060999/1‑101 (MQ=255)             cgACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAgg                                                                                   <  2:1399585/101‑1 (MQ=255)              gACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGa                                                                                  >  1:749800/1‑101 (MQ=255)              gACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGa                                                                                  <  1:11726/101‑1 (MQ=255)                cAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGt                                                                                <  1:977941/101‑1 (MQ=255)                 aaaaaTGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGtt                                                                               >  2:382533/1‑101 (MQ=255)                     aTGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCt                                                                           >  1:1171949/1‑101 (MQ=255)                        ccAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGAt                                                                        <  2:41328/101‑1 (MQ=255)                         cAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATg                                                                       <  2:939813/101‑1 (MQ=255)                          agagAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGc                                                                      <  1:853953/101‑1 (MQ=255)                               aCTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAg                                                                 <  1:1002716/101‑1 (MQ=255)                                 tGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCt                                                               <  2:1060994/101‑1 (MQ=255)                                  gAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTg                                                              <  2:1032595/101‑1 (MQ=255)                                       tGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAat                                                         <  2:952897/101‑1 (MQ=255)                                         gCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAatga                                                       <  1:1487094/101‑1 (MQ=255)                                          cGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAatgat                                                      <  2:749795/101‑1 (MQ=255)                                           gAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAatgatg                                                     >  1:445418/1‑101 (MQ=255)                                           gAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAatgatg                                                     >  2:1459213/1‑101 (MQ=255)                                           gAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAatgatg                                                     <  2:213803/101‑1 (MQ=255)                                                 ggCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGc                                               <  1:266040/101‑1 (MQ=255)                                                  gCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCg                                              >  1:1665463/1‑101 (MQ=255)                                                     aTTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGtga                                           >  1:386130/1‑101 (MQ=255)                                                     aTTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGtga                                           >  2:385116/1‑101 (MQ=255)                                                     aTTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGtga                                           >  2:45601/1‑101 (MQ=255)                                                       tACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGtgatg                                         >  2:406977/1‑101 (MQ=255)                                                       tACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGtgatg                                         >  1:79589/1‑101 (MQ=255)                                                         cTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTgatgat                                       >  2:1526334/1‑101 (MQ=255)                                                         cTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTgatgat                                       <  2:1402754/101‑1 (MQ=255)                                                         cTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTgatgat                                       >  1:901114/1‑101 (MQ=255)                                                              aTGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGc                                  >  1:890601/1‑101 (MQ=255)                                                              aTGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGc                                  <  2:1451265/101‑1 (MQ=255)                                                                        gtcAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGgtcgt                        >  2:1162683/3‑101 (MQ=255)                                                                        gtcAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGgtcgt                        >  1:506450/3‑101 (MQ=255)                                                                        gtcAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGgtcgt                        >  2:666957/3‑101 (MQ=255)                                                                          cAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGgtcgtcg                      <  1:247443/101‑1 (MQ=255)                                                                          cAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGgtcgtcg                      >  2:1038581/1‑101 (MQ=255)                                                                           aGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGgtcgtcgt                     >  2:865566/1‑101 (MQ=255)                                                                            gCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGgtcgtcgtc                    <  1:1012988/101‑1 (MQ=255)                                                                              gTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCgg                  >  1:1378284/1‑101 (MQ=255)                                                                               tGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGt                 >  2:737692/1‑101 (MQ=255)                                                                                 aTGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTg               <  1:170385/101‑1 (MQ=255)                                                                                 aTGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTg               >  1:1584787/1‑101 (MQ=255)                                                                                      gctgcGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACAt          <  2:618686/101‑1 (MQ=255)                                                                                        tgcGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCa        <  2:972121/101‑1 (MQ=255)                                                                                            ctcgATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAAGGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGCCGGTTGCACATCAAAGc    <  2:42208/98‑1 (MQ=255)                                                                                              cgATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAg  <  1:1495765/100‑1 (MQ=255)                                                                                              cgATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAg  <  1:791275/100‑1 (MQ=255)                                                                                              cgATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAg  >  1:621289/2‑101 (MQ=255)                                                                                                   |                                                                                                CAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGCCAGCGTGATGTTGCTGCACTGGATGAAAAATACACGAAGGAGTTAGCTAATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAG  >  193/1990‑2183

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: ‑