Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	790,832	Δ1 bp	coding (198/1047 nt)	galT ←	galactose‑1‑phosphate uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	790,832	0	G	.	95.1%	169.7 / 1.7	41	coding (198/1047 nt)	galT	galactose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base . (22/17);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01

CATGGGTACGTTTTCCCCAGTTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCC  >  NC_000913/790590‑791073                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  cATGGGTTCGTTTTCCCCAGTTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCTGG‑TTTTTATCg                                                                                                                                                                                                                                          <  2:417606/251‑1 (MQ=255)     ggTACGTTTTCCCCAGTTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTg                                                                                                                                                                                                                                      <  1:49172/251‑1 (MQ=255)                            gTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAACGACATTAGGGAAAACGTAAGTCCCGGCGTAATCGGG‑TTTTTATCGCCTGTAACCCGCAATTTACCttctca                                                                                                                                                                                                                >  2:203905/1‑245 (MQ=255)                                  tCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGa                                                                                                                                                                                                           >  1:524040/1‑249 (MQ=255)                                   ccTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAAT‑AGGGTTTTTATCGCCTGTCACCCGCCCATTACCTGCGCAGGGGAAGc                                                                                                                                                                                                        >  2:470197/1‑251 (MQ=255)                                         aGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTg                                                                                                                                                                                                  <  1:70421/251‑1 (MQ=255)                                          ggTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCGGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTgg                                                                                                                                                                                                 <  2:29696/251‑1 (MQ=255)                                               tGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCg                                                                                                                                                                                             <  1:184199/250‑1 (MQ=255)                                                       tCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCCCCCGCACATTACCTGCGCCGAGGAAGCAATCTGGATCGTGCGCAgg                                                                                                                                                                                     <  2:160082/250‑1 (MQ=255)                                                             aaTGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACAcc                                                                                                                                                                              <  1:213648/251‑1 (MQ=255)                                                                   gcAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTtg                                                                                                                                                                        >  1:97304/1‑251 (MQ=255)                                                                   gcAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTtg                                                                                                                                                                        <  2:229590/251‑1 (MQ=255)                                                                        gCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTtggctgg                                                                                                                                                                   >  1:375466/1‑251 (MQ=255)                                                                                ctGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑ATTTTATCGCCTGTCACCCGCACATTACCTGCGCCGAGGAAGCAATCTAGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCt                                                                                                                                                           >  1:484205/1‑251 (MQ=255)                                                                                           ttACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTAACTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTg                                                                                                                                                 <  1:162622/250‑1 (MQ=255)                                                                                            tACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTTCGCCCCCTGcc                                                                                                                                               >  2:83973/1‑251 (MQ=255)                                                                                             aCTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGGAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGcc                                                                                                                                               >  2:295439/1‑250 (MQ=255)                                                                                                    tCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGAACGTGACTTTCTGGCGCATCTGGCGAGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCCGGGgctc                                                                                                                                       >  1:192504/1‑249 (MQ=255)                                                                                                      cGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCtt                                                                                                                                     <  1:522791/251‑1 (MQ=255)                                                                                                        gTGAAAATCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCATCGCATCAGCGGATCGTGTCTTGCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCTGGTCACCCGCACATTTCCCGCTCTGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGTCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAg                                                                                                                                   <  1:491908/251‑1 (MQ=255)                                                                                                                   aTCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGCCGCATCTGGCGAGTCAGACATCCAAGCCGCAAAGTCCTTCGTGAAAAAGTAAGTCCCGGAGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCgcggtgcggg                                                                                                                        >  1:351612/1‑250 (MQ=255)                                                                                                                    tCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGAACGTGACTTTCTGGCGCCTCTGGCGTGTCAGACATCAAAGCCGCACAGTCATTAGTGAAAACGTACGTCACGGTGTCATCGGG‑TTTTTATCGCATGTCACCCGCACCTTCACTGCGCAGAGGCAGCAATCTGGATCGGGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCACTGCCAGGGGCGCTTAGCGcggtgcggtg                                                                                                                       >  1:503145/1‑251 (MQ=255)                                                                                                                         cGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAAcc                                                                                                                  >  1:104577/1‑251 (MQ=255)                                                                                                                                  ccGCTCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCTCGGTGCGGTGAAACCAGAATCCAt                                                                                                         <  1:195267/251‑1 (MQ=255)                                                                                                                                                 ccactcATCAGCGGTTCGAGACTTTCTGGCGCAACTGGCGTGTCAGACATCAAATCCGCAAAGTCATTAGTGAAAACTTATGTCCCGGTGTAATCGTG‑TTTTTATCGCCTGTCACCCGCACATTTCCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTGGCGCGGTGCGGCGAAACCAGAATCCATTGCCCGGTGAGCggg                                                                                          <  2:90747/246‑1 (MQ=255)                                                                                                                                                       aTCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAg                                                                                    <  1:71153/251‑1 (MQ=255)                                                                                                                                                              gATCGTGACTTTCTGGCGCATCAGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTACCGGTGTAATCGGG‑TTTTTATCGCCTGTCAACCGAACATTACCTGCGCAGAGGAAGCCATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGGGAAACCAGAACCCATTGCCCGGTGGACGGGTGGTAgcggagcg                                                                             >  2:360101/1‑247 (MQ=255)                                                                                                                                                               aTCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGAt                                                                             >  1:110105/1‑251 (MQ=255)                                                                                                                                                               aTCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGAtg                                                                            >  1:223313/1‑251 (MQ=255)                                                                                                                                                                                cATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGAt                                                           >  1:130742/1‑251 (MQ=255)                                                                                                                                                                                   cTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTa                                                         >  1:119489/1‑250 (MQ=255)                                                                                                                                                                                   cTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGAGATGGGGATGCTCAAAGGGATTa                                                         >  2:463180/1‑250 (MQ=255)                                                                                                                                                                                    tGGCGTGTCAGACATCAAAGCCGCAAAGTCCTTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTaa                                                        >  1:504121/1‑250 (MQ=255)                                                                                                                                                                                                 aTCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCa                                                                                                                 >  2:17664/1‑180 (MQ=255)                                                                                                                                                                                                 aTCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCa                                                                                                                 <  1:17664/180‑1 (MQ=255)                                                                                                                                                                                                 aTCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTc                                          <  1:470197/251‑1 (MQ=255)                                                                                                                                                                                                   cAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGt                                        >  1:414408/1‑251 (MQ=255)                                                                                                                                                                                                   cAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGt                                        >  1:359042/1‑251 (MQ=255)                                                                                                                                                                                                                                     ccGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCggtac      >  2:107660/1‑247 (MQ=255)                                                                                                                                                                                                                                     ccGGTGTAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGAc      <  1:413047/251‑1 (MQ=255)                                                                                                                                                                                                                                          gtAATCGGG‑TTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGGTCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACtgcc  >  2:418214/1‑250 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  CATGGGTACGTTTTCCCCAGTTCTGCGGTTTGCTCCTGCCAGGTTTTGACGATTTCCGTCAATGCTGCAACGCTGAGCTCTGGCAGCGTTTTACTGTGATCCGGTGAAAAGCAGATCACCCGGCTGGTGCCGCGCGCGCTCTGGCAACGCATCAGCGGATCGTGACTTTCTGGCGCATCTGGCGTGTCAGACATCAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCC  >  NC_000913/790590‑791073

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