Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	410,922:1	+G	6.2%	coding (779/909 nt)	mak →	manno(fructo)kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	410,922	1	.	G	6.2%	263.4 / 13.8	65	coding (779/909 nt)	mak	manno(fructo)kinase
Reads supporting (aligned to +/- strand):  ref base . (32/29);  new base G (2/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

TCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATT  >  NC_000913/410775‑411067                                                                                                                                                     |                                                                                                                                                  tCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGTTCGGGATGAGCAATGTAGACCGTTt‑a                                                                                                                                                  <  3:120116/149‑1 (MQ=255)          ttGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGACGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACg                                                                                                                                         >  2:224350/1‑149 (MQ=255)             aagtaAGCGATCCGGTAGCGGATCTGGCATTGCGTCGGTAAGAGCTGCGGCTGGCAAAATCGCTGGCAGATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGt                                                                                                                                       <  6:38705/148‑1 (MQ=255)             aagaaTGCGATCCGGTAGCGTAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACAGtt                                                                                                                                      <  4:205931/149‑1 (MQ=255)              agaaAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTg                                                                                                                                     <  1:280554/149‑1 (MQ=255)              agaaAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTg                                                                                                                                     <  2:230446/149‑1 (MQ=255)               gaaAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTTGATAT‑AAACGGtcta                                                                                                                                    <  6:351015/149‑4 (MQ=255)               gaaAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACgg                                                                                                                                        <  6:222987/145‑1 (MQ=255)               gaaAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACgg                                                                                                                                        >  5:222987/1‑145 (MQ=255)               gaaAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCAGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTgg                                                                                                                                    >  4:168831/1‑149 (MQ=255)                   gCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAg                                                                                                                                >  5:272870/1‑149 (MQ=255)                       tCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGa                                                                                                                            <  3:319628/149‑1 (MQ=255)                       tCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGa                                                                                                                            <  4:314381/149‑1 (MQ=255)                         cGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGtt                                                                                                                              >  7:185228/1‑145 (MQ=255)                         cGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGtt                                                                                                                              <  8:185228/145‑1 (MQ=255)                                  aagtGGCATTGCGTCGCTACGAGCTGCGGATGGCAAAATCGCTGGCACATGTCTTGAATATTATCGATCCGGATGTGATTGTCCTGGGCGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAAttt                                                                                                                 <  5:188242/146‑1 (MQ=255)                                  aaCTGGCATCGCGTCGCTACGAGCCGCGTCTGGCAAAATCGCTGGCCCCTGTCGCGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGGTGATTAAACAAttt                                                                                                                 <  4:47220/149‑1 (MQ=255)                                     tGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTc                                                                                                              >  2:85671/1‑149 (MQ=255)                                                 gCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGa                                                                                                   <  2:96498/148‑1 (MQ=255)                                                     cGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAAtgt                                                                                               <  6:272870/148‑1 (MQ=255)                                                        gCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGaaa                                                                                           <  1:85671/149‑1 (MQ=255)                                                              gCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCgg                                                                                     <  3:168831/149‑1 (MQ=255)                                                               cTGGCAAAATCGCTGGCACATGACGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGt                                                                                    <  1:224350/149‑1 (MQ=255)                                                                  gCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGTGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGACGGTGCg                                                                                 >  2:219812/1‑149 (MQ=255)                                                                            tGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAg                                                                       >  7:282809/1‑149 (MQ=255)                                                                                  aTGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTc                                                                                                           >  6:84978/1‑107 (MQ=255)                                                                                  aTGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTc                                                                                                           <  5:84978/107‑1 (MQ=255)                                                                                    ggcgggAATATTCTCGATCCGGATGTGATTGCCCTGGGGGGCGGGAGGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACg                                                                   <  3:155412/140‑1 (MQ=255)                                                                                    gTCGTGAATATTCTCGATCCGGATTTCATTGTACTCGGGGGAGGCATGAGCAATATACAACCATT‑ATATCAAACGTTATCGCACTAGATTAAACAATTTGTCTTCTGCTCCGAATCTTAAAAGCAGGTCCTTAAGACGAAGCACg                                                                   >  4:155412/1‑145 (MQ=255)                                                                                       gTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAAt                                                                                                 >  7:287609/1‑112 (MQ=255)                                                                                       gTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAAt                                                                                                 <  8:287609/112‑1 (MQ=255)                                                                                         gAATATTCTCGCTCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGaaa                                                                                           >  5:148436/1‑116 (MQ=255)                                                                                         gAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGaaa                                                                                           <  6:148436/116‑1 (MQ=255)                                                                                          aaTATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAATCGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAg                                                         >  8:301401/1‑149 (MQ=255)                                                                                            tatTCTCGATCCGGATGTGATTGTCGTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAgcg                                                       >  5:124074/1‑149 (MQ=255)                                                                                                 tcGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTa                                                  >  3:65914/1‑149 (MQ=255)                                                                                                 tcGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTa                                                  >  4:196743/1‑149 (MQ=255)                                                                                                 tcGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTa                                                  >  1:153230/1‑149 (MQ=255)                                                                                                 tcGATCCGGATGTGATTGTCCTGGGGGGCGGCATGAGCAATGTAGACCGTTT‑ATATCAAACGGTAGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCCGTGCGTAAGGCGAAGCACGTTGATTCCAGCGTCGTa                                                  >  7:299719/1‑149 (MQ=255)                                                                                                    aTCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTAtgc                                               >  3:100079/1‑149 (MQ=255)                                                                                                    aTCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTAcgc                                               <  8:282809/149‑1 (MQ=255)                                                                                                       cGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACgcggc                                            >  1:216195/1‑149 (MQ=255)                                                                                                         gATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGgcgc                                          >  7:31558/1‑149 (MQ=255)                                                                                                               aTTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAATCGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGt                                     <  7:301401/148‑1 (MQ=255)                                                                                                               aTTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTg                                    >  6:111329/1‑149 (MQ=255)                                                                                                               aTTGGCCTGGGGGGCGGGATTAGCAATGTAGACCGTTT‑ATATCAATCGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGt                                     <  8:31558/148‑1 (MQ=255)                                                                                                                   tCCTGGGGGGCGGGATGAGCAATGAAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAAtgt                                                                                               >  7:177454/1‑86 (MQ=255)                                                                                                                   tCCTGGGGGGCGGGATGAGCAATGAAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAAtgt                                                                                               <  8:177454/86‑1 (MQ=255)                                                                                                                       ggggggCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGc                            >  5:354676/1‑149 (MQ=255)                                                                                                                            gCGGGATGTGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACaa                       <  1:219812/149‑1 (MQ=255)                                                                                                                            gCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCcac                         >  7:117727/1‑147 (MQ=255)                                                                                                                            gCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCcac                         >  4:28082/1‑147 (MQ=255)                                                                                                                            gCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACaa                       <  1:363258/149‑1 (MQ=255)                                                                                                                                gATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAagag                    >  7:105303/1‑148 (MQ=255)                                                                                                                                     gCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTaaaa               >  5:302743/1‑148 (MQ=255)                                                                                                                                     gCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTaaaa               >  5:297582/1‑148 (MQ=255)                                                                                                                                      cAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTaaaaaa             >  5:156184/1‑149 (MQ=255)                                                                                                                                         gcccaaCCGTTTGATAT‑AAACGGTCGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACg           >  2:354535/6‑148 (MQ=255)                                                                                                                                           tAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTa         <  5:111329/148‑1 (MQ=255)                                                                                                                                           tAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAg        <  6:124074/149‑1 (MQ=255)                                                                                                                                           ccaaCCGTTTGATAT‑AAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAg        >  1:187997/4‑149 (MQ=255)                                                                                                                                           ccaaCCGTTTGATAT‑AAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAg        <  2:187997/146‑1 (MQ=255)                                                                                                                                             gACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAgg       >  2:328516/1‑148 (MQ=255)                                                                                                                                             gACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAgg       >  7:263393/1‑148 (MQ=255)                                                                                                                                                cGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAAt   <  3:196743/149‑1 (MQ=255)                                                                                                                                                 gTTT‑ATAGCAAACGGTTGGGCAGTTGATTAAACAATTTTTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCAGGGTGATTGCAGCGGCGTACGCGGCGCTGGGTGGTTATGGCCACAAGAGTAAAATACGTAGGCAAtt  <  8:299719/149‑1 (MQ=255)                                                                                                                                                  ttt‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAg        <  3:154131/142‑1 (MQ=255)                                                                                                                                                  ttt‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAg        >  4:154131/1‑142 (MQ=255)                                                                                                                                                     |                                                                                                                                                  TCCGCCTGGTTGAAGAAAGCGATCCGGTAGCGGAACTGGCATTGCGTCGCTACGAGCTGCGGCTGGCAAAATCGCTGGCACATGTCGTGAATATTCTCGATCCGGATGTGATTGTCCTGGGGGGCGGGATGAGCAATGTAGACCGTTT‑ATATCAAACGGTTGGGCAGTTGATTAAACAATTTGTCTTCGGCGGCGAATGTGAAACGCCGGTGCGTAAGGCGAAGCACGGTGATTCCAGCGGCGTACGCGGCGCTGCGTGGTTATGGCCACAAGAGTAAAAAACGTAGGCAATT  >  NC_000913/410775‑411067

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