Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	2,150,261	+C	coding (766/960 nt)	SAUSA300_RS10960 ←	fructokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	2,150,257	1	.	C	96.0%	204.7 / 3.6	74	coding (770/960 nt)	SAUSA300_RS10960	fructokinase
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base C (4/67);  total (4/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

GCTAAAATTTCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑CCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGGCAGCACCATCTGCTCCCTTCGT  >  NC_007793/2150158‑2150351                                                                                                     |                                                                                                gCTAAAATTTCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTg                                                                                                  <  1:393694/100‑1 (MQ=255) gCTAAAATTTCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTg                                                                                                  <  1:2812709/100‑1 (MQ=255) gCTAAAATTTCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTg                                                                                                  <  1:1503133/100‑1 (MQ=255)  cTAAAATTTCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTga                                                                                                 <  1:5380740/100‑2 (MQ=255)  cTAAAATTTCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCAG‑c                                                                                                <  1:5719708/100‑1 (MQ=255)   tAAAATTTCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑cg                                                                                               <  1:4382703/100‑2 (MQ=255)   tAAAATTTCCTCACCTTCATTTTCAAATAACTGACTTAAATTTAACACATCTATTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑cc                                                                                               <  1:2496202/100‑1 (MQ=255)     aaaTTTCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑ccag                                                                                             <  1:796377/100‑3 (MQ=255)       aTTTCCTCACCTTCATTTTCAAATAACTGAGTTAAACTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑cctggt                                                                                           <  1:2776966/100‑5 (MQ=255)         ttCCTCACCTTCATTTTCAAATAACTGAGTTGAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑cccccagt                                                                                         <  1:2336497/100‑5 (MQ=255)         ttCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑cccccagt                                                                                         <  1:1979171/100‑5 (MQ=255)         ttCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑cccccagt                                                                                         <  1:5282197/100‑5 (MQ=255)         ttCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑cccccagt                                                                                         <  1:5210432/100‑5 (MQ=255)         gatcTCACCCTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑cccccagt                                                                                         <  1:1589560/97‑5 (MQ=255)          tCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑cccccagtt                                                                                        <  1:829723/100‑6 (MQ=255)          tCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑cccccagtt                                                                                        <  1:3018927/100‑6 (MQ=255)          tCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑cccccagta                                                                                        <  1:2089249/100‑6 (MQ=255)            ttCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGt                                                                                      <  1:3712490/99‑1 (MQ=255)            tctcccTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGt                                                                                      <  1:3197552/96‑1 (MQ=255)            gatctcTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGt                                                                                      <  1:5959268/95‑1 (MQ=255)            cttaCCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGt                                                                                      <  1:2993521/97‑1 (MQ=255)            cgatctttCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGt                                                                                      <  1:5483783/94‑1 (MQ=255)            cTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGt                                                                                      <  1:30401/100‑1 (MQ=255)            cTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGt                                                                                      <  1:1299077/100‑1 (MQ=255)            cTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGt                                                                                      <  1:5550097/100‑1 (MQ=255)            cTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGt                                                                                      <  1:2084760/100‑1 (MQ=255)            cTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCCGCCCCC‑AGTTGt                                                                                      <  1:1118232/100‑1 (MQ=255)            atcttcTTCGTTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGt                                                                                      <  1:3605410/95‑1 (MQ=255)             tcgccTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTa                                                                                     <  1:5524232/97‑1 (MQ=255)             tCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTa                                                                                     <  1:5730437/100‑1 (MQ=255)              cACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTac                                                                                    <  1:5717713/100‑2 (MQ=255)               aCCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATc                                                                                   <  1:3981318/100‑1 (MQ=255)               aCCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATc                                                                                   <  1:2312502/100‑1 (MQ=255)               aCCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATc                                                                                   <  1:897342/100‑1 (MQ=255)               aCCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATc                                                                                   <  1:403788/100‑1 (MQ=255)               aCCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATc                                                                                   <  1:3459937/100‑1 (MQ=255)                ccTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCa                                                                                  <  1:899744/100‑1 (MQ=255)                ccTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCa                                                                                  <  1:141404/100‑1 (MQ=255)                    cATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAg                                                                              >  1:1260864/1‑100 (MQ=255)                       tttCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGtt                                                                           <  1:6037201/100‑1 (MQ=255)                               aaCTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTATACCTta                                                                   >  1:4606839/1‑100 (MQ=255)                               aaCTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑CCCC‑AGTTGTATCAACAGGTTTTACCTtat                                                                  <  1:1107638/100‑1 (MQ=255)                                    aGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGcc                                                              <  1:5996634/100‑1 (MQ=255)                                    aGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGcc                                                              <  1:3076377/100‑1 (MQ=255)                                    aGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGcc                                                              <  1:1536781/100‑1 (MQ=255)                                    aGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGcc                                                              <  1:3199573/100‑1 (MQ=255)                                    aGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGcc                                                              <  1:1051998/100‑1 (MQ=255)                                      ttAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCAc                                                            <  1:850861/100‑1 (MQ=255)                                       tAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACt                                                           <  1:5593855/100‑1 (MQ=255)                                       tAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACt                                                           <  1:3744163/100‑1 (MQ=255)                                       tAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACt                                                           <  1:442497/100‑1 (MQ=255)                                       tAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACt                                                           <  1:1375919/100‑1 (MQ=255)                                       tAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACt                                                           <  1:3990264/100‑1 (MQ=255)                                        aaaTTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTa                                                          <  1:2110802/100‑1 (MQ=255)                                        aaaTTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTa                                                          <  1:4503940/100‑1 (MQ=255)                                        aaaTTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTa                                                          <  1:704605/100‑1 (MQ=255)                                        aaaTTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTa                                                          <  1:4851137/100‑1 (MQ=255)                                        aaaTTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTa                                                          <  1:3829274/100‑1 (MQ=255)                                           tttAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑CCCCTAGTTGTATCAACAGGTTTTACCTTATAGCCACTATgg                                                       <  1:920222/100‑2 (MQ=255)                                            ttAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGAt                                                      <  1:2677132/100‑1 (MQ=255)                                             tAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGacg                                                     <  1:1387336/100‑3 (MQ=255)                                             tAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATa                                                     <  1:5469452/100‑1 (MQ=255)                                                  caTCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATaacccc                                                <  1:4216874/100‑5 (MQ=255)                                                  caTCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGAt                                                <  1:3164523/100‑1 (MQ=255)                                                    tCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGAttt                                              <  1:4413539/100‑2 (MQ=255)                                                    tCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGAttg                                              <  1:5793263/100‑2 (MQ=255)                                                      ggTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAAt                                            <  1:4143453/100‑1 (MQ=255)                                                           cTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCag                                       <  1:2885474/100‑2 (MQ=255)                                                             aaaaTCCTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCAtta                                     <  1:5280055/100‑2 (MQ=255)                                                                  ccTGCTAATAACTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCaa                                <  1:2290315/100‑1 (MQ=255)                                                                           atcTGCACCGATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAatatggggc                       <  1:840999/98‑6 (MQ=255)                                                                           aaCTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgc                       <  1:575469/100‑1 (MQ=255)                                                                           aaCTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgc                       <  1:4804006/100‑1 (MQ=255)                                                                           aaCTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgc                       <  1:2960385/100‑1 (MQ=255)                                                                           aaCTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgc                       <  1:715353/100‑1 (MQ=255)                                                                             cTGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgcat                     <  1:53749/100‑2 (MQ=255)                                                                              tGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgcagc                    <  1:2910071/100‑1 (MQ=255)                                                                              tGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgcagc                    <  1:4306548/100‑1 (MQ=255)                                                                              tGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgcagc                    <  1:3084739/100‑1 (MQ=255)                                                                              tGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgcagc                    <  1:5928060/100‑1 (MQ=255)                                                                              tGCACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgcagc                    <  1:538232/100‑1 (MQ=255)                                                                               tgcccAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTTCCTTATAGCCACTATGATAATGATTAATACCATTTTTCGAATATACGgcagca                   <  1:6010649/97‑1 (MQ=255)                                                                               gCACCGATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgcacca                   >  1:717927/1‑100 (MQ=255)                                                                                cACCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGgcagcat                  <  1:3821450/100‑2 (MQ=255)                                                                                 aCCAATAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGGCAGCAcc                 <  1:1749664/100‑1 (MQ=255)                                                                                   tgctcttccgaTCTCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGGCAGCACCAt               <  1:2106400/90‑1 (MQ=255)                                                                                    aaTAAATGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGGCAGCACCATc              <  1:5408198/100‑1 (MQ=255)                                                                                        aaTGCATCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGGCAGCACCATCTGCt          <  1:227580/100‑1 (MQ=255)                                                                                             ctCCCCTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGGCAGCACCATCTGCTCCCtt     <  1:5915703/99‑1 (MQ=255)                                                                                               gatcTG‑CCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGGCAGCACCATCTGCTCCCTTCGt  <  1:858678/97‑1 (MQ=255)                                                                                               ccccTGCCCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGGCAGCACCATCTGCTCCCTTCg   >  1:5416833/1‑100 (MQ=255)                                                                                                     |                                                                                                GCTAAAATTTCCTCACCTTCATTTTCAAATAACTGAGTTAAATTTAACACATCGGTTGCTAAAATCCTGCTAATAACTGCACCAATAAATGCATCCCCTG‑CCCC‑AGTTGTATCAACAGGTTTTACCTTATAGCCACTATGATAATGATTAATACCATTTTTCAAATATACGGCAGCACCATCTGCTCCCTTCGT  >  NC_007793/2150158‑2150351

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