Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,393,552	G→A	F172F (TTC→TTT)	queG ←	epoxyqueuosine reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,393,552	0	G	A	96.5%	178.8 / ‑1.4	60	F172F (TTC→TTT)	queG	epoxyqueuosine reductase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (25/31);  total (28/32)
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.41e-01

GGCAGATAAGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGTTTTTCAGC  >  NC_000913/4393308‑4393792                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 ggCAGATAAGCTGGCAGTCATAGCAACCGTAATTACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGGCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATAGACCAACAATTCGCCTAAAAaaaagacc                                                                                                                                                                                                                                            >  2:289705/1‑249 (MQ=255)             ggCAGCCACCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACGCTTCCGGGATCGCCCCTTCAAGTTCGATGGGGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCCCGTCGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTcg                                                                                                                                                                                                                                <  2:252461/251‑1 (MQ=255)                aGTCATCGCAACCGTAAATACGTTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGa                                                                                                                                                                                                                             <  1:128553/251‑1 (MQ=255)                   cATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTg                                                                                                                                                                                                                          >  1:91165/1‑251 (MQ=255)                   cATCACAACCGTAAATCCGGTTTCCCAGGAAAGGTCGCAACTCTTCCGGTAGCGCCCATTCAAGTGCGAGGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGAtt                                                                                                                                                                                                                           <  2:118721/250‑1 (MQ=255)                     tCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGGGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTAATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTg                                                                                                                                                                                                                          >  1:415861/1‑249 (MQ=255)                        cAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGAt                                                                                                                                                                                                                     >  1:280349/1‑251 (MQ=255)                         aaCCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATa                                                                                                                                                                                                                    <  1:249985/251‑1 (MQ=255)                               aaaTACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAagtg                                                                                                                                                                                                                >  1:18304/1‑249 (MQ=255)                               aaaTACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGGTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGGCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAagtgag                                                                                                                                                                                                              <  1:193496/251‑1 (MQ=255)                                     ggTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCCCGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAgg                                                                                                                                                                                                                       >  2:189856/1‑236 (MQ=255)                                     ggTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAgg                                                                                                                                                                                                                       <  1:189856/236‑1 (MQ=255)                                            cATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTcc                                                                                                                                                                                                 >  1:89703/1‑251 (MQ=255)                                             aTTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGGTCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGGGCTTACCTGTcc                                                                                                                                                                                                 >  2:123090/1‑250 (MQ=255)                                               tAACGGTCGCAACTCTGCCGGGCTCGCCACTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGCAGCTGGAACCGGTCGGGAAAATCATCATACAGGCCACGCATTTGCCGCATCCTTCCTCGCCTGGTTGATCCCCGGGCAGCGGAATGTCGAACAGCAATTCGCCTAAAAACACGAACGAACCGGCCTCGCGATTGAGGCAAAGTGACAGCTGACCTGTccggc                                                                                                                                                                                              >  2:46518/1‑248 (MQ=255)                                                       gCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGCCCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCa                                                                                                                                                                                      <  1:361112/251‑1 (MQ=255)                                                        cAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCCATTCGCCTAAAAAAAAGAACGAACCGGCCTcgcgctt                                                                                                                                                                                                                           >  2:336784/1‑210 (MQ=255)                                                        cAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGAtt                                                                                                                                                                                                                           <  1:336784/213‑1 (MQ=255)                                                          aCTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGCCCAGCCGAGCCCAGCt                                                                                                                                                                                   >  1:164576/1‑251 (MQ=255)                                                               tCCGGGATCGCCCCTTCAATTTCGATGGGGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCACGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCttt                                                                                                                                                                                 <  1:271994/248‑1 (MQ=255)                                                               tCCGGGAGCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCGTCGACGGTAGATGGCTCGACGATGGCACCGGGCGGGCAAATCGTCATACAGGCCACGCAGTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCttttt                                                                                                                                                                               <  2:31563/250‑1 (MQ=255)                                                                  gggATCGCCCCTTCAAGTTCGATGGGGAGATAAGTGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTGCAGc                                                                                                                                                                           <  2:21321/251‑1 (MQ=255)                                                                   ggATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCt                                                                                                                                                                          <  1:17696/251‑1 (MQ=255)                                                                        gCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTCTATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGACCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGCATATAGACCAGCAATTCGCCTAAAAAACAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGCGCCCCGCTTTTTACGCTAACgg                                                                                                                                                                     >  2:242533/1‑251 (MQ=255)                                                                               caatttcgagggggataTAAGAGATACAGCGGCGGGAATCGACGGTAAAGGGCTCGACTAGGGACCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCGTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTGGCGATTGAGGATAAGTGAGTGCTTCCCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCg                                                                                                                                                              <  2:174964/235‑1 (MQ=255)                                                                               cAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAAGATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGAGCTc                                                                                                                                                               <  2:294840/250‑1 (MQ=255)                                                                                aaGTTTGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCCCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGTCCAGCTTTTTCAGCTAACGGGCGCTCga                                                                                                                                                             <  2:18304/251‑1 (MQ=255)                                                                                       aTGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAgg                                                                                                                                                      <  1:242533/251‑1 (MQ=255)                                                                                         ggTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGGCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGgcg                                                                                                                                                    >  2:137666/1‑251 (MQ=255)                                                                                                 aaGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAA‑AAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGATCGAGAATAGGCGCAGAATCg                                                                                                                                            >  1:260241/1‑250 (MQ=255)                                                                                                         ccgcggcgAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACgg                                                                                                                                    >  2:12895/3‑251 (MQ=255)                                                                                                                       gACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGAtt                                                                                                                                                                                                                           >  1:148780/1‑150 (MQ=255)                                                                                                                       gACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGAtt                                                                                                                                                                                                                           <  2:148780/150‑1 (MQ=255)                                                                                                                         cGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGCGAATAGGCGCAGAATCGACAAACGGACTAAAATTCAGCGa                                                                                                                     >  2:354223/1‑250 (MQ=255)                                                                                                                          ggTATATGACTCGACGATGGCACCGGTCGGCCAACTCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGCATATCGACCAGCACTTCGCCTAAAAAAAAGAACGACCCGGCCCCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCCAGCCCCGCTTTTTCAGCTAACGGGCCCTCGGGAATAGGCGCAGCATCGACAACGGGTCTAAAATTCaaccaa                                                                                                                    >  2:194041/1‑245 (MQ=255)                                                                                                                                    tCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGCGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTg                                                                                                         >  2:18504/1‑251 (MQ=255)                                                                                                                                     cgacgaTGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCGCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGtt                                                                                                          <  2:415861/249‑1 (MQ=255)                                                                                                                                      gacgaTGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCt                                                                                                       <  1:12895/251‑1 (MQ=255)                                                                                                                                            ggCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAatca                                                                                                 <  2:91165/251‑1 (MQ=255)                                                                                                                                               aCCGGGCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGTTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCAtct                                                                                              <  1:377205/251‑1 (MQ=255)                                                                                                                                                           aaTCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCtttt                                                                                  >  2:97691/1‑251 (MQ=255)                                                                                                                                                            aTCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCttttt                                                                                 >  1:71914/1‑251 (MQ=255)                                                                                                                                                              cGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGa                                                                               <  1:97691/251‑1 (MQ=255)                                                                                                                                                                tCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGt                                                                             >  1:218699/1‑251 (MQ=255)                                                                                                                                                                      aGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGCCAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTg                                                                       >  1:178959/1‑251 (MQ=255)                                                                                                                                                                          cACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCa                                                                   <  1:123090/251‑1 (MQ=255)                                                                                                                                                                           aCGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTgcgc                                                                    >  2:76643/1‑249 (MQ=255)                                                                                                                                                                              cATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAg                                                               <  1:367295/251‑1 (MQ=255)                                                                                                                                                                              cATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATAATCTCGCCCAGATTTTTGAGTCGGTTGCGCAGa                                                                 >  2:209414/1‑249 (MQ=255)                                                                                                                                                                                        aTCCTTCCCCGACTGGTTGAGCCACGGGCAGCGGAATATCGACCAGCAATCCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAg                                                     <  2:128979/251‑1 (MQ=255)                                                                                                                                                                                            ttcctCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCATCGCGATTGAGGATAAGTGAATGCTTACCTGTCCAGCCGAGCACAGCTTTTTCAGCTAACGGGAGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTACGCGAAACACAATGGTGCTGAATCATCCCGCCCAGCTTTTTGAGTCGGTTGCCCGGAAGTTTGTGAtcatca                                                  >  2:314880/1‑245 (MQ=255)                                                                                                                                                                                             tcctCGACTGGTTGATCCACGGGCAGCGGAATACCGCCCAGCAATTCGCCTAAAAAAAAGAACGAACCGTCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCAc                                                 <  2:280349/250‑1 (MQ=255)                                                                                                                                                                                                 cGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGccc                                            <  2:178959/251‑1 (MQ=255)                                                                                                                                                                                                          gaaccACGGGCAGCTGAATATCGACCAGCAAGTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATaa                                   <  1:46518/248‑1 (MQ=255)                                                                                                                                                                                                             ccACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGCTTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCAc                                                 >  1:213372/1‑234 (MQ=255)                                                                                                                                                                                                             ccACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGATGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCTGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCAc                                                 <  2:213372/234‑1 (MQ=255)                                                                                                                                                                                                               aCGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCt                              >  1:158171/1‑251 (MQ=255)                                                                                                                                                                                                               aCGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTGGTGATAGTCACGGCCCAGCGCATAACGGCt                              <  1:196651/251‑1 (MQ=255)                                                                                                                                                                                                                   gCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACa                          <  1:76643/251‑1 (MQ=255)                                                                                                                                                                                                                                          tCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGTTTTTCa    <  2:71914/250‑1 (MQ=255)                                                                                                                                                                                                                                           cGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGCCAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCTGGTTTGGGGTTTTTCAGc  >  1:56063/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 GGCAGATAAGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGTTTTTCAGC  >  NC_000913/4393308‑4393792

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