Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	65,043	A→C	L246L (CTT→CTG)	polB ←	DNA polymerase II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	65,043	0	A	C	87.9%	134.4 / 5.7	59	L246L (CTT→CTG)	polB	DNA polymerase II
Reads supporting (aligned to +/- strand):  ref base A (4/3);  new base C (32/19);  total (36/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.57e-01

GCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCCTGAGCGACAGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCGCTATTATCGCGCCCAAGACGCAGCGGAAGACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGAT  >  NC_000913/64801‑65279                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             gcAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCCTGAGCGACAGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCGCTATTATCGCGCCCAAGACGCAg                                                                                                                                                                                                                                      <  1:43893/251‑1 (MQ=255)       ttATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCCCACTGTTATCGCGCCCCAGCCGCAGCGGaat                                                                                                                                                                                                                                >  1:64654/1‑250 (MQ=255)             tCGGCGAAACGGCGGTCAATTTCGTCCATTCGAACCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTGTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGt                                                                                                                                                                                                                           <  1:315894/250‑1 (MQ=255)              cGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCCCCCCCGCCGCAGCGGAATACGGTa                                                                                                                                                                                                                          >  1:229418/1‑250 (MQ=255)                  gAAACGGCGGTCAATTTCCTCCCTTCGTTCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATTGAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAACCGACCTTTAGCC‑TGGGCAAAAAAAACGCCGTTTTTAAAGCCGTGCTCGCGCCCCTCCCGCTCACTGTTATCGCGCCCCCGCCGCCGCGGAATACGGTAACGcc                                                                                                                                                                                                                     >  2:58536/1‑250 (MQ=255)                   aaaCGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTc                                                                                                                                                                                                                    <  1:248861/251‑1 (MQ=255)                       ggcggTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGc                                                                                                                                                                                                                 >  1:384918/1‑250 (MQ=255)                              aaTTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGttttt                                                                                                                                                                                                         >  2:375269/1‑251 (MQ=255)                                tttCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATCCGGTAACGCTCGGCATGTTTTTgc                                                                                                                                                                                                       >  1:167010/1‑251 (MQ=255)                                 ttCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTgc                                                                                                                                                                                                       <  1:159008/250‑1 (MQ=255)                                           cGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGCCCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGc                                                                                                                                                                                              >  1:68758/1‑249 (MQ=255)                                                  aCGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCCGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAAc                                                                                                                                                                                     >  1:224245/1‑251 (MQ=255)                                                   cGGGTTATCGCTAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGCAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTACCGCTCGGCCTGTTTTTgaagaat                                                                                                                                                                                                  >  2:436815/1‑231 (MQ=255)                                                   cGGGTTATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTCCCAGCGAGAATGAAGAGAAATTCCAGAACGCTGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCAt                                                                                                                                                                                                  <  1:436815/237‑1 (MQ=255)                                                    agattATAGATAGATTTTCCTTCGCATAATAACTCCGGGGCGACGGTTTACAGAGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGGGCCTCGAGACCGTCGATAATTAGACGACCTTGAGAC‑TGGGCAACAAAGACGGCGTTTTTAAAGCCGTGCTCGCGCTACTCCAGCACACTGATATCGCGCACCAGCCGCAGCGGAATCCGGTAACGCTAGGCATGTTATTGCAGACTTCGCAGCTAGAACtc                                                                                                                                                                                   >  2:276648/4‑250 (MQ=255)                                                       ttATCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGcac                                                                                                                                                                                >  2:86426/1‑251 (MQ=255)                                                         aTCGATAGATTTTCCTTCGCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTGGCAGCATGCGCAGATCGAACGGcacc                                                                                                                                                                               <  1:129183/250‑1 (MQ=255)                                                                        ttCGCCTAATAGCTCCTGAGCGACAGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCGCTATTATCGCGCCCAAGACGCAGCGGAAGACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCAAACTGCACCACGTTCCACCCgatg                                                                                                                                                                >  2:31334/1‑250 (MQ=255)                                                                         tCGCCTAATAACTCCTGAGCTACTGTTGCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATCCCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACTCCGTTTTTAAAGCCGTGCTCGCGCCACTCTAGCTCACTTTTAGCGCGCCCCAGCCGCAGCGGAATACGGTAATGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCgatga                                                                                                                                                               <  1:373119/250‑1 (MQ=255)                                                                           gCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATcc                                                                                                                                                            >  1:27931/1‑250 (MQ=255)                                                                           gCCTAATAACTCCTGAGCGACGGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATca                                                                                                                                                            <  1:246427/251‑1 (MQ=255)                                                                                                     ccAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCa                                                                                                                                   <  1:58536/250‑1 (MQ=255)                                                                                                       aGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGc                                                                                                                                >  2:140129/1‑251 (MQ=255)                                                                                                         cGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTACCGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGCTGATCACATCAGGATCGTAGTTGGCAAACCAGGCGt                                                                                                                              >  2:395723/1‑251 (MQ=255)                                                                                                                    agaAATTCCAGAACGCGGATTTCAGCGCCGCGATACCGCCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTcc                                                                                                                   <  1:433952/251‑1 (MQ=255)                                                                                                                                aCGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTGTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTgcggg                                                                                                       <  1:113202/251‑1 (MQ=255)                                                                                                                                           aGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGAc                                                                                            <  1:140129/251‑1 (MQ=255)                                                                                                                                               ccTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAACGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGCCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGCCGTAt                                                                                        >  1:132493/1‑251 (MQ=255)                                                                                                                                                        cGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTca                                                                               <  1:195725/251‑2 (MQ=255)                                                                                                                                                           cGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCGCTATTATCGCGCCCAAGACGCAGCGGAAGACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACt                                                                                                                                                                                    <  1:234706/147‑1 (MQ=255)                                                                                                                                                           cGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCGCTATTATCGCGCCCAAGACGCAGCGGAAGACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACt                                                                                                                                                                                    >  2:234706/1‑147 (MQ=255)                                                                                                                                                            gATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTCTTTAAAGCCGTGCACGCGCCACTCCACCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCACCCGATGATCACCTCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTCCCAGTAGCTGCGGGCGGCTGCCGACggcttcacgctcaaaat                                                                           >  2:311905/1‑235 (MQ=255)                                                                                                                                                              tAATTAGCCGCCCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGCCGTATTCCAGTTCAAAATcg                                                                         >  2:361300/1‑250 (MQ=255)                                                                                                                                                               aaTTAGCCGCCCTTTAGCC‑TGGGCAAAAAAGTCGCCGTTTTAAAACCCGTGCTCGCGCCACGCCAGCTCACTGTTAGCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATcg                                                                         <  2:224245/250‑2 (MQ=255)                                                                                                                                                                aatagtctaCCTGTAGCCGTGGCCAAAAAAGCTGCGGTTTTTAA‑GCCGTGCCCGCGCCACCCCAGCTCTCTGTTAGCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGCTTTGGCAGCAGTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACAGCAGGATCGTATTTGGCAAACCAGGCGTTGATTTCTGCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAggtcacaatcggg                                                                       <  2:426816/243‑13 (MQ=255)                                                                                                                                                                       gcccTTTAGCC‑TGGGCAAAAAAGCCGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCGAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCgcggat                                                                <  2:64654/249‑2 (MQ=255)                                                                                                                                                                       gACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCgcgg                                                                  >  2:92475/1‑246 (MQ=255)                                                                                                                                                                       gACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACGGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCgcggag                                                                <  1:243382/251‑1 (MQ=255)                                                                                                                                                                         ccTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATcac                                                                                                                                                           >  1:425008/1‑158 (MQ=255)                                                                                                                                                                         ccTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGGGCTCGCGCCACGCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATcac                                                                                                                                                           <  2:425008/158‑1 (MQ=255)                                                                                                                                                                             aaGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGAGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAAAGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCAACACGATCCAACCGATGATCACATCAGGATCGTAGATGGCAAAACAGGAGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAATTCGAGCGCGGAGGCGTCg                                                          >  1:455459/2‑251 (MQ=255)                                                                                                                                                                                 c‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCCTTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACAGATTCCAGTTCAAAATCGAGCGCGGAGGCGTAGCCAt                                                      >  2:241127/1‑251 (MQ=255)                                                                                                                                                                                   tGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTAGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACAAGGCGTTGAGTTTTTCCAGGAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATAGAGCGCGGAGGCGTCGCCAtt                                                     >  1:57941/1‑251 (MQ=255)                                                                                                                                                                                     ggCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATtct                                                   >  1:304728/1‑251 (MQ=255)                                                                                                                                                                                      gCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCGCTATTATCGCGCCCAAGACGCAGCGGAAGACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCAAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATtctc                                                  <  1:31334/251‑1 (MQ=255)                                                                                                                                                                                       cAAAAAAGACGCCGTTTTTAAAGCAGTGCTCGCGCCACTCAAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATtctct                                                 >  1:392228/1‑250 (MQ=255)                                                                                                                                                                                              gACGCCGTTTTTAAAGCAGTGCTCGCGCAACTCAAGCTGGCTATTATCGCGCCAAAGACGCAGCGGAAGACGGTAACGCTCGGCATGTTTTTGCAGCATTAGCAGATCGAACTGCACCACGTTCCAACAGATGATCACATCAGGATCGTAGTTGGCAAACAAGGCGTTGAGTTTTTCCAGCAATTGCGGGCGGCTGGCGACGTATTACAGTTCGAAATCAAGCGAGGAGGCGTCGCAATTCTCCGGcacc                                           >  1:150908/1‑247 (MQ=255)                                                                                                                                                                                                 gCCGTTTTTAAAGCCGTGCTCGCGCCAATCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGc                                        >  1:359698/1‑250 (MQ=255)                                                                                                                                                                                                       tttAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCATTAGCTGCGGGCGGCTGGCGACTTATTCCAGTTCGAAATCGAGCGAGGAGGCGTCGCCATTCTCTGGCCCCAGAATATaa                                  >  2:340101/1‑250 (MQ=255)                                                                                                                                                                                                       tttAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATaa                                  >  1:294193/1‑250 (MQ=255)                                                                                                                                                                                                        ttAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGCGTTTTTCCAGTAGCTGCGGGCGGCTGGCGTCGTATTCCAGTTCAAACTCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCGGCATATAAAc                                >  2:386700/1‑251 (MQ=255)                                                                                                                                                                                                         tAAAGCCGTGCTCGCGCCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAAc                                <  1:371396/250‑1 (MQ=255)                                                                                                                                                                                                                        gcCACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGc                  >  1:67694/1‑249 (MQ=255)                                                                                                                                                                                                                          cACTCCAGCTCACTGTTATCGCGCCCCAGCCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCCCCATTCTCTGGCCCCAGCATATAAACCATCCGCTGCCCGCAGcc              >  2:237215/1‑251 (MQ=255)                                                                                                                                                                                                                             tCCAGCTCGCTATTATCGCGCCCAAGACGCCGCGGAAGACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATCTAAACAATGCGCTGCCCGCAGCCTTc           >  2:171719/1‑251 (MQ=255)                                                                                                                                                                                                                             tCCAGCTCACTGTTATCGCGCCCCAGCCGCCGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAACATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAACCGATGCGCTGCCCGCAGCCTTc           >  2:366417/1‑251 (MQ=255)                                                                                                                                                                                                                             tCCAGCTCACTGTTATCGCGCCCCAACCGCCGCGGAATACGGTCACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATCTAAACGATGCGCTGCCCGCAGCCTTc           >  1:362278/1‑251 (MQ=255)                                                                                                                                                                                                                             cccAGCTCCCTGTTATCGCGCCCCAGTCGCAGCGGAATACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCGATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTc           <  1:126160/250‑1 (MQ=255)                                                                                                                                                                                                                                       tgtctccgcgCCCCAGCCGCAGCGGAATACGGTAACGCCCGGCATGGGTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGAGCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGTAGCTGCGGGCGGCTGGCGACGTATTCCAGTTCAAAATCGAGCGCGGAGGCGTCGCCATTCTCTGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGAt  <  2:386715/244‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             GCAGGTTTATCTTCGGCGAAACGGCGGTCAATTTCGTCCATTCGATCCCACGGGTTATCGATAGATTTTCCTTCGCCTAATAGCTCCTGAGCGACAGTTTCCAGCGAGAATGAAGAGAAATTCCAGAACGCGGATTTCAGCGCCTCGATACCGTCGATAATTAGCCGACCTTTAGCC‑TGGGCAAAAAAGACGCCGTTTTTAAAGCCGTGCTCGCGCCACTCCAGCTCGCTATTATCGCGCCCAAGACGCAGCGGAAGACGGTAACGCTCGGCATGTTTTTGCAGCATTCGCAGATCGAACTGCACCACGTTCCAACCGATGATCACATCAGGATCGTAGTTGGCAAACCAGGCGTTGAGTTTTTCCAGCAACTGCGGGCGGCTGGCGACGTATTCCAGTTCGAAATCAAGCGAGGAGGCGTCGCCATTCTCCGGCCCCAGCATATAAACGATGCGCTGCCCGCAGCCTTCCAGGCCGAT  >  NC_000913/64801‑65279

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: ‑