Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,790,949	C→T	G99D (GGC→GAC)	selO ←	protein adenylyltransferase SelO

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,790,949	0	C	T	97.7%	126.0 / NA	45	G99D (GGC→GAC)	selO	protein adenylyltransferase SelO
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (17/26);  total (18/27)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGC‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATACCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAAATCAGCCGGGCATTATTTAAAGG  >  NC_000913/1790703‑1791175                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                    aaaTGACCAAAGCGCAGATGACTTGGTGCCACACGCCGCAGCATCGCGCCTGCCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCGCTTATCGTCGACCGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATc                                                                                                                                                                                                                                <  1:320631/251‑1 (MQ=255)     gACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCagc                                                                                                                                                                                                                             >  1:570028/1‑250 (MQ=255)     gACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTAGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCagc                                                                                                                                                                                                                             >  1:429772/1‑250 (MQ=255)                  aTGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGc                                                                                                                                                                                                               >  1:138574/1‑251 (MQ=255)                      cTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCCGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCg                                                                                                                                                                                                             >  1:488334/1‑249 (MQ=255)                        tGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGt                                                                                                                                                                                                          >  1:207253/1‑250 (MQ=255)                                        aGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATc                                                                                                                                                                                         >  1:54197/1‑251 (MQ=255)                                                gcCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTCACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCCGCCAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGtt                                                                                                                                                                                 >  1:39212/1‑251 (MQ=255)                                                    ggCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATGCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTgg                                                                                                                                                                               <  1:238906/249‑1 (MQ=255)                                                      cTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGgcccg                                                                                                                                                                           >  1:199157/1‑251 (MQ=255)                                                        ccgccgTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGgcccgcc                                                                                                                                                                         >  1:69696/1‑251 (MQ=255)                                                           ccgTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGTACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAg                                                                                                                                                                      <  1:144992/251‑1 (MQ=255)                                                                     aTACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACt                                                                                                                                                            <  1:2411/251‑1 (MQ=255)                                                                        cacTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACtgat                                                                                                                                                         <  1:462308/251‑1 (MQ=255)                                                                              cgaatgccTGGGGACGAAACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCAAGACTTTCTCGTATCGTCGAACGTAAAACGGCCCGTCCTTCCCCCCATCGCGAATAAAGCATCAGGCCAGCACCTTTCAGATCCCAGCACATTGTA‑TACACATCAACAAGCAGATGTTCGACACGTAGAACGCCCCGCCCATCCCCCAGTCAGCCAGCCCCGACGCCGAGCTGATGACCAc                                                                                                                                                   >  2:196421/8‑251 (MQ=255)                                                                              cGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGACGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGATGGCCCGCCCAGACGCCGAACTGATGAcacc                                                                                                                                                   >  2:196423/1‑248 (MQ=255)                                                                               ggaTCGCTGGTGACGATACTGAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACt                                                                                                                                                  <  1:387786/249‑1 (MQ=255)                                                                                 atcgcggttgacgttccTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCACCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGGTCGCCGAGTAAAATGCCGCGCCCATCACCCAGGTGGCCCGCCCAGACGCCGAACTGATGACCACTGt                                                                                                                                                <  2:492753/235‑1 (MQ=255)                                                                                          gccgatatgtAACGCGCGGGCCTGCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTGCTCGTAGCGTCCAACGTAAAACCGCCCGTCCATCACCCAGTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTgcg                                                                                                                                       <  2:92809/242‑1 (MQ=255)                                                                                            cGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGcc                                                                                                                                     <  1:501718/251‑1 (MQ=255)                                                                                                   tAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAgtgg                                                                                                                                <  1:14719/249‑1 (MQ=255)                                                                                                    aaCGCGCGTGCCGCCGGAATGCCCATATAATGCATCTCCTCACTGGCGATCCTTGCACGTATCGTCGAACGTAAATCCGCCCGTCCAGCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAACCCTGCGCCAgtggtg                                                                                                                              <  2:207253/250‑1 (MQ=255)                                                                                                    aaCGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCCTCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCTCCCAGCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGCCCACCGTAAACCTGCACCAggggtac                                                                                                                             >  2:380934/1‑249 (MQ=255)                                                                                                    aaCGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCCCCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGa                                                                                                                             <  1:177955/251‑1 (MQ=255)                                                                                                     ccgcgcgGGTCGTCGGAATGCCCAGATAATGCAACGCCTCACTGGCGAGTCTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTTATGACCACTGTAACCCTGCGCCAGTGGTGAc                                                                                                                            <  2:391094/250‑1 (MQ=255)                                                                                                               cgtaggAATGCCCAGATAATGCACCGCCTCACTGGCGAAACTTTCTCGGAACGTCGAACGTAAAACCGCCCGTCCATCACCCATGCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCa                                                                                                                  <  2:174653/247‑1 (MQ=255)                                                                                                                          ccAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTcg                                                                                                       >  1:142675/1‑251 (MQ=255)                                                                                                                                            cACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTgc                                                                                      <  1:101202/250‑1 (MQ=255)                                                                                                                                                   gagaCTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCAtttt                                                                              >  1:297314/1‑251 (MQ=255)                                                                                                                                                        tttCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAAc                                                                         >  1:99349/1‑251 (MQ=255)                                                                                                                                                            ccgtaccgccgAAGGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGcg                                                                     <  2:259022/242‑1 (MQ=255)                                                                                                                                                             cGTATTGTCGAACGTCAAACCGCCCGTCCAACACCCATGCTCGAATAAGGCGGCAGGTCAGCACCGTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTTGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTTCGCTAGTGGTGACATATCAGGCAGTATGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGcga                                                                    <  2:519848/251‑1 (MQ=255)                                                                                                                                                                        aCGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATAc                                                         <  1:453536/251‑1 (MQ=255)                                                                                                                                                                            aaaaCCGCCCGTCCATCACCCATTCGCGAGTAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAg                                                     <  1:36170/251‑1 (MQ=255)                                                                                                                                                                                          aTCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCATTTGTACGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGATCTGATGACCACTGTAAACCTGCGCCAGTGGTTACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGtt                                       <  2:497300/251‑1 (MQ=255)                                                                                                                                                                                                 aTTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATAGCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGtatt                                <  1:355306/251‑1 (MQ=255)                                                                                                                                                                                                       gAATACGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCAAGTTGGCCCGCCCAGACGACGAACTGATGACCACTGTCAACCTGCGCCAGTGGTGAAATATCCGGCTGTAGGGCTTCGCCGCCCCAGACACCTACGCCATTTTTAAACAGCGGCGCTGGa                                                             >  2:455814/1‑216 (MQ=255)                                                                                                                                                                                                       gAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGACTAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGCCACCTGCGCCATTTTTAAACAGCGACGATGGa                                                             <  1:455814/216‑1 (MQ=255)                                                                                                                                                                                                         aTAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCaa                         >  1:10312/1‑250 (MQ=255)                                                                                                                                                                                                                   aGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCa                                                      <  2:441438/211‑1 (MQ=255)                                                                                                                                                                                                                   aGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCa                                                      >  1:441438/1‑211 (MQ=255)                                                                                                                                                                                                                      ccagcatcgtacAGATGCCATGCCAGTGTAGAGA‑CCTCAGCAAGCAGTTGTTCGCCGAGTAAAATTCCGCGGCCATCCCCCAGTTGGCCCGCCCAGACGCCGAACTGAGGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGGAGGTCTTCGCCGCCCCAGACGCCTGTGTCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCAAAATCCGCCGGGCatt           <  2:437715/240‑1 (MQ=255)                                                                                                                                                                                                                         gCACCTTTCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAAATCAGCCGGGCatgatt        <  1:339308/251‑1 (MQ=255)                                                                                                                                                                                                                           aCCTTTCAGATGCCATTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAAATCAGCCGGGCATTATTTa      <  1:27955/251‑1 (MQ=255)                                                                                                                                                                                                                                tCAGATGCCAGTCCATTGTAGTGT‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATATCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAAATCAGCCGGGCATTATTTAAAgg  <  1:78347/250‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                    AAATGACCAAAGCGCAGATGACTTGGTGCCACACGCATCAGCATCGCGCCTGGCTCCGCCGTTTCCCGATACACTGGCGAATCGCTGGTGACGATACTTAACGCGCGGGTCGTCGGAATGCCCAGATAATGCATCGCCTCACTGGCGAGACTTTCTCGTATCGTCGAACGTAAAACCGCCCGTCCATCACCCATTCGCGAATAAGGCGTCAGGCCAGCACCTTTCAGATGCCAGTCCATTGTAGTGC‑CATCAGCAAGCAGTTGTTCGCCGAGTAAAATGCCGCGCCCATCACCCAGTTGGCCCGCCCAGACGCCGAACTGATGACCACTGTAAACCTGCGCCAGTGGTGACATACCAGGCAGTAGGGCTTCGCCGCCCCAGACACCTGCGCCATTTTTAAACAGCGACGATGGAATACTCAGCGTGTTAGCCAGTTCGGTATTATGCCAAATCAGCCGGGCATTATTTAAAGG  >  NC_000913/1790703‑1791175

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