Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007793	748,914	T→G	V298G (GTA→GGA)	SAUSA300_RS03615 →	aldo/keto reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007793	748,914	0	T	G	100.0%	252.3 / NA	96	V298G (GTA→GGA)	SAUSA300_RS03615	aldo/keto reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (90/6);  total (90/6)

GGAACAAGTCAGTTAAAGCGTATTGATCAAGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGTACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGCATATCATGACAAACGAAGATAAACGTTTCGAACAATTAA  >  NC_007793/748817‑748998                                                                                                  |                                                                                     ggAACAAGTCAGTTAAAGCGTATTGGTCAAGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGAc                                                                                    >  1:422271/1‑100 (MQ=255)  tgaCAAGTCAGTTAAAGCGTATTGATCAAGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACa                                                                                   >  1:1824995/3‑100 (MQ=255)    aCAAGTCAGTTAAAGCGTATTGATCAAGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAg                                                                                 >  1:3016654/1‑100 (MQ=255)    aCAAGTCAGTTAAAGCGTATTGATCAAGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAg                                                                                 >  1:2451425/1‑100 (MQ=255)      aaGTCAGTTAAAGCGTATTGATCAAGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGat                                                                               <  1:1022495/100‑1 (MQ=255)          cAGTTAAAGCGTATTGATCAAGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGatcaga                                                                           >  1:3147014/1‑96 (MQ=255)             ttAAAGCGTATTGATCAAGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGt                                                                        <  1:2887986/100‑1 (MQ=255)               aaaGCGTATTGATCAAGCCATCTCCGGGCTACCCCTTAATTTAGATTCTCCTTCGTGGTTTGACATTTCCACCTCTATTATCGGACACGatctgaagagc                                                                      >  1:1849252/1‑91 (MQ=255)                          tagaattaaTCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACt                                                           >  1:2360112/8‑100 (MQ=255)                          ctCAAGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACt                                                           >  1:2557756/2‑100 (MQ=255)                           tCAAGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACtt                                                          >  1:1071861/1‑100 (MQ=255)                           ccAAGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACtt                                                          >  1:907406/2‑100 (MQ=255)                           ataaGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACtt                                                          >  1:4796278/3‑100 (MQ=255)                           acAAGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACtt                                                          >  1:2851130/2‑100 (MQ=255)                            cAAGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACttt                                                         >  1:3484051/1‑100 (MQ=255)                            cAAGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACttt                                                         >  1:4615435/1‑100 (MQ=255)                            cAAGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACttt                                                         >  1:3823673/1‑100 (MQ=255)                            cAAGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACttt                                                         >  1:5478191/1‑100 (MQ=255)                            cAAGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACttt                                                         >  1:4353693/1‑100 (MQ=255)                            cAAGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACttt                                                         >  1:5222493/1‑100 (MQ=255)                            cAAGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACttt                                                         >  1:4467489/1‑100 (MQ=255)                            cAACCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACttt                                                         >  1:2419773/1‑100 (MQ=255)                            aaaGCAATCGAAGGGGTACAACTTAATTTAGATGATCAGTCGCGGCTTTACGTTGACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACttt                                                         >  1:4410885/2‑100 (MQ=255)                             caGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACtttt                                                        >  1:5596169/2‑100 (MQ=255)                             aaGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACtttt                                                        >  1:4415801/1‑100 (MQ=255)                             aaGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACtttt                                                        >  1:1464778/1‑100 (MQ=255)                             aaGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACtttt                                                        >  1:3705433/1‑100 (MQ=255)                             aaGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACtttt                                                        >  1:1650872/1‑100 (MQ=255)                             aaGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACtttt                                                        >  1:4271059/1‑100 (MQ=255)                             aaGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACtttt                                                        >  1:2608508/1‑100 (MQ=255)                              tgCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTa                                                       >  1:1897752/2‑100 (MQ=255)                              tctgtgcGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGCTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACttttg                                                       >  1:3116896/7‑99 (MQ=255)                              aGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTa                                                       >  1:1779471/1‑100 (MQ=255)                              aGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTa                                                       >  1:1844071/1‑100 (MQ=255)                              aGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTa                                                       >  1:3970202/1‑100 (MQ=255)                              aGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTa                                                       >  1:4517037/1‑100 (MQ=255)                               gCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTaa                                                      >  1:312082/1‑100 (MQ=255)                               gCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTaa                                                      >  1:282206/1‑100 (MQ=255)                                cAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTaaa                                                     >  1:1221987/1‑100 (MQ=255)                                 caTCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAAt                                                    >  1:4571431/2‑100 (MQ=255)                                  ctCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATc                                                   >  1:4849621/2‑100 (MQ=255)                                    tgtggggCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCAt                                                 >  1:1024407/5‑100 (MQ=255)                                    tctaGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCAt                                                 >  1:2631174/4‑100 (MQ=255)                                    cttaGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCAt                                                 >  1:888756/4‑100 (MQ=255)                                    cGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACATCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCAt                                                 >  1:3147458/1‑100 (MQ=255)                                    cGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCAt                                                 >  1:5552206/1‑100 (MQ=255)                                    cGAAGGGCTACAACTTAATTTAGATGATCAGTCGCGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAcatcat                                                 >  1:3366265/1‑100 (MQ=255)                                     gAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATa                                                <  1:3438018/100‑1 (MQ=255)                                      caGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATaa                                               >  1:1611431/2‑100 (MQ=255)                                      caGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATaa                                               >  1:4903900/2‑100 (MQ=255)                                      aaGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATaa                                               >  1:708876/1‑100 (MQ=255)                                      aaGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATaa                                               >  1:2440590/1‑100 (MQ=255)                                      aaGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATaa                                               >  1:2748890/1‑100 (MQ=255)                                      aaGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATaa                                               >  1:4221528/1‑100 (MQ=255)                                      aaGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATaa                                               >  1:2840549/1‑100 (MQ=255)                                      aaGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATaa                                               >  1:3896299/1‑100 (MQ=255)                                      aaGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATaa                                               >  1:5021678/1‑100 (MQ=255)                                      aaGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATaa                                               >  1:2380825/1‑100 (MQ=255)                                      aaGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTCAAATCagat                                               >  1:4800283/1‑97 (MQ=255)                                      aaGGGCCACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGCAAAATTATTTACTTTTAAATCATaa                                               >  1:4780378/1‑100 (MQ=255)                                       tgggCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATaaa                                              >  1:415421/2‑100 (MQ=255)                                       taggCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATaaa                                              >  1:2688889/3‑100 (MQ=255)                                       aGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATaaa                                              >  1:1087774/1‑100 (MQ=255)                                         tgCTCCAACTTAATTTAGATGTTCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGTTATTCCGTAAACTTATTTACTTTTAAATCATaaaag                                            >  1:453246/2‑98 (MQ=255)                                         ggCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAgg                                            >  1:2912787/1‑100 (MQ=255)                                         ggCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATCATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAga                                            >  1:4488153/1‑99 (MQ=255)                                          gCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAGATcggaagagc                                           >  1:2935213/1‑92 (MQ=255)                                          gCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGa                                           >  1:1213420/1‑100 (MQ=255)                                          acTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGa                                           >  1:4148857/2‑100 (MQ=255)                                           cTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAg                                          >  1:3003516/1‑100 (MQ=255)                                           cTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAg                                          >  1:2016119/1‑100 (MQ=255)                                            tACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAGGGAGc                                         <  1:1133550/100‑1 (MQ=255)                                            tACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGc                                         <  1:1241099/100‑1 (MQ=255)                                            tACAACTTAATTTAGATGATCAGTCGTGGATTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGc                                         <  1:1685057/100‑1 (MQ=255)                                             aCAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAgat                                        >  1:3767950/1‑98 (MQ=255)                                             aCAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGCa                                        >  1:4541202/1‑100 (MQ=255)                                             aCAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGCa                                        >  1:4403138/1‑100 (MQ=255)                                             aCAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGCa                                        >  1:1462632/1‑100 (MQ=255)                                             aCAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGCa                                        >  1:4176128/1‑100 (MQ=255)                                             aCAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGCa                                        >  1:5552784/1‑100 (MQ=255)                                             aCAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGCa                                        >  1:952281/1‑100 (MQ=255)                                             aCAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGCa                                        >  1:2905437/1‑100 (MQ=255)                                             aCAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGCa                                        >  1:1263597/1‑100 (MQ=255)                                             aCAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAAACATaagatcgga                                        >  1:837913/1‑93 (MQ=255)                                             aCAACTTAATTTAGATGATCAGTCGTGGTTCGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGCa                                        >  1:929164/1‑100 (MQ=255)                                                   tAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAgagatcgga                                  >  1:1319165/1‑92 (MQ=255)                                                      acaggatgatCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGCATATCATGAc                               >  1:4806189/5‑100 (MQ=255)                                                       ttAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGCAGATcggaaga                              >  1:4899976/1‑94 (MQ=255)                                                        tAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGCATATCATGAcga                             >  1:1044560/1‑98 (MQ=255)                                                          gatgatCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGCATATCATGACAAAc                           >  1:2387225/1‑100 (MQ=255)                                                             gatCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGCATATCATGAcgaacgaa                        >  1:2996774/1‑100 (MQ=255)                                                             gatCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGCATATCATGACaaagatc                        >  1:4014080/1‑96 (MQ=255)                                                             gatCAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGCATATCATGACAAACGaa                        >  1:98743/1‑100 (MQ=255)                                                                cAGTCGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGCATATCATGACAAACGAAGAt                     >  1:3803191/1‑100 (MQ=255)                                                                   gcGTGGTTTGACATTTACACCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGCATATCATGACGAACGAAGAtcgg                  >  1:1642164/2‑97 (MQ=255)                                                                                   acacCGCTATTATCGGACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGCATATCATGACAAACGAAGATAAACGTTTCGAACAAttta  >  1:648744/1‑98 (MQ=255)                                                                                                  |                                                                                     GGAACAAGTCAGTTAAAGCGTATTGATCAAGCAATCGAAGGGCTACAACTTAATTTAGATGATCAGTCGTGGTTTGACATTTACACCGCTATTATCGTACAAGATATTCCGTAAACTTATTTACTTTTAAATCATAAAGGAGCATATCATGACAAACGAAGATAAACGTTTCGAACAATTAA  >  NC_007793/748817‑748998

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