Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,392,997	G→A	V357V (GTC→GTT)	queG ←	epoxyqueuosine reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,392,997	0	G	A	98.1%	181.8 / NA	52	V357V (GTC→GTT)	queG	epoxyqueuosine reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (27/24);  minor base C (1/0);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCAGGCGAATAACTCAATGAGTTCCGGT  >  NC_000913/4392750‑4393243                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       acaGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACt                                                                                                                                                                                                                                                     <  1:100138/251‑1 (MQ=255)                        tttCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTAt                                                                                                                                                                                                                              >  1:333381/1‑251 (MQ=255)                         ttCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTAt                                                                                                                                                                                                                              <  2:39287/250‑1 (MQ=255)                           cAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTAttt                                                                                                                                                                                                                            <  1:107232/250‑1 (MQ=255)                                aaaCTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTgcg                                                                                                                                                                                                                      <  1:271927/251‑1 (MQ=255)                                 aaCTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTgcgc                                                                                                                                                                                                                     >  1:97160/1‑251 (MQ=255)                                              gAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAATCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTGCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCt                                                                                                                                                                                                        <  1:305792/251‑1 (MQ=255)                                              gAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCt                                                                                                                                                                                                        >  1:385710/1‑251 (MQ=255)                                               aaTAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCtt                                                                                                                                                                                                       >  2:226415/1‑250 (MQ=255)                                                       ttttttAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTc                                                                                                                                                                                                <  1:103954/250‑1 (MQ=255)                                                             aaCCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCgag                                                                                                                                                                                          <  1:226415/250‑1 (MQ=255)                                                                         gataAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTGTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCa                                                                                                                                                                             <  2:97009/251‑1 (MQ=255)                                                                          ataAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCAc                                                                                                                                                                            <  2:334639/251‑1 (MQ=255)                                                                                   aaaTTAAGATTCCGTTTTGTCGCGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGACGTGGGTGCTCACCTTTACgcc                                                                                                                                                                   >  2:360023/1‑249 (MQ=255)                                                                                        aaGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTcc                                                                                                                                                              >  2:304335/1‑251 (MQ=255)                                                                                            ttCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAgcg                                                                                                                                                          >  1:453111/1‑251 (MQ=255)                                                                                               cGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTg                                                                                                                                                       <  1:387484/251‑1 (MQ=255)                                                                                                  tttGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTc                                                                                                                                                     <  1:8333/250‑1 (MQ=255)                                                                                                             gaTCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTc                                                                                                                                          >  1:73186/1‑250 (MQ=255)                                                                                                                cACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATcc                                                                                                                                      <  2:385710/251‑1 (MQ=255)                                                                                                                    tGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCa                                                                                                                                    >  2:59756/1‑249 (MQ=255)                                                                                                                                   tCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCACTCACCCGAACCAGACGCTGTTTTTTCGGCCATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGCGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATg                                                                                                                   >  1:45824/1‑251 (MQ=255)                                                                                                                                      tGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGGGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCCCCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATg                                                                                                                   <  2:333381/248‑1 (MQ=255)                                                                                                                                                 gCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTa                                                                                                     <  1:304335/251‑1 (MQ=255)                                                                                                                                                      tttatttttattCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACGTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCa                                                                                                 <  1:2913/250‑1 (MQ=255)                                                                                                                                                       ttatttttattCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTAc                                                                                                               <  1:286716/235‑1 (MQ=255)                                                                                                                                                       ttatttttattCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTAc                                                                                                               >  2:286716/1‑235 (MQ=255)                                                                                                                                                          tttttattAACAAGATGTGGATGAATCAGGCGTCACGCGGATACCCTTTTTCCATCACCCGACCCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCACTCGCCCACGCTATGTGCTCCTCGTGATGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGCGCCTTGCCTAATGCTACGGCGATCTTACGCAGCCac                                                                                            >  2:321884/1‑250 (MQ=255)                                                                                                                                                            tttattCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAAc                                                                                           <  1:397071/250‑1 (MQ=255)                                                                                                                                                                     aaGCTGTGGCTGCATCAGGCGTCACGCGGCAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCACCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCACCGCAGGTGAcc                                                                                 >  1:391169/1‑251 (MQ=255)                                                                                                                                                                     aaGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGAcc                                                                                 >  1:61308/1‑251 (MQ=255)                                                                                                                                                                          gtgGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATac                                                                            >  1:445660/1‑251 (MQ=255)                                                                                                                                                                          gtgGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATa                                                                             <  2:328075/250‑1 (MQ=255)                                                                                                                                                                                        cGTCACGCGGTAACCCTTTTTCAATCCCCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGAt                                                              <  2:341206/251‑1 (MQ=255)                                                                                                                                                                                        cGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGAt                                                              >  2:370263/1‑251 (MQ=255)                                                                                                                                                                                             cgcgGTAACCCTTGTTCAATCCCCCGAACCAGACGCTGTTTTTGCGGCAATTGCACTTCAACTATGCCCGCATTTCGTCGCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTGCCAGCGCTGTCAAAATCGTGTCATCCCAAGGGGCATTGCCGAATGCTACGGCGATATTACGCAGCCAAGGCAGGTGGCCAATACCACGAATCGCCTTTCCTTc                                                         <  2:198298/251‑1 (MQ=255)                                                                                                                                                                                               cgGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCg                                                       >  1:78822/1‑251 (MQ=255)                                                                                                                                                                                                   aaCCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGAc                                                   >  1:167797/1‑251 (MQ=255)                                                                                                                                                                                                    aCCCTTTTTCAATCACCCGAACCAGACGCTGTTTTT‑CGGCAATTGCACTTTAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCt                                                                                                                                                                                                        >  2:33419/1‑100 (MQ=255)                                                                                                                                                                                                      ccTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACttt                                                <  2:97160/251‑1 (MQ=255)                                                                                                                                                                                                                  aCCCGAACGAGACGCTGTTTTTTCGGCAATTGCACTTCCACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGGGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTtctc                                    >  2:326223/1‑251 (MQ=255)                                                                                                                                                                                                                              cGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGTCTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAACAACTTCTCTTCGCTCCAGGc                        >  2:366211/1‑251 (MQ=255)                                                                                                                                                                                                                               gCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCAGGCg                       >  1:374874/1‑251 (MQ=255)                                                                                                                                                                                                                               gCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTACGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCTACGCAGGGGACCACTACGACGAGTCGCCGATCCTTCCGTGAATTTTAAAAACTTCTCTTCGCTCCAGGCg                       >  2:180360/1‑251 (MQ=255)                                                                                                                                                                                                                                      tttCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCAGGCGAATAAc                 <  2:4319/250‑1 (MQ=255)                                                                                                                                                                                                                                            cAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCACCGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCAGGCGAATACATCAATg           >  2:456881/1‑250 (MQ=255)                                                                                                                                                                                                                                            cAATTGCACTTCAACAAGGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCAGTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACATGTAAAAACTTCTCTTCGCTCCAGTCGAATAACTCAATGa          <  2:45824/251‑1 (MQ=255)                                                                                                                                                                                                                                             aaTTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGCGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCAGGCGAATAACTCAATg           >  2:77381/1‑249 (MQ=255)                                                                                                                                                                                                                                             aaTTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCAGGCGAATAACTCAATg           >  2:130774/1‑249 (MQ=255)                                                                                                                                                                                                                                                 gCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTGCCGTGACTTTTAAAAACTTCTCTTCGCTCCAGGCGAATAACTCAATGAGTTcc     <  2:61308/251‑1 (MQ=255)                                                                                                                                                                                                                                                    cTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAACATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCTCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCAGGCGAATAACTCAATGAGTTCCGGt  >  2:406456/1‑251 (MQ=255)                                                                                                                                                                                                                                                    cTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGACTCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCAGGCGAATAACTCAATGAGTTCCGGt  >  2:297001/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       ACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCAGGCGAATAACTCAATGAGTTCCGGT  >  NC_000913/4392750‑4393243

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