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	94.3%	138.0 / ‑1.2	54	V357V (GTC→GTT)	queG	epoxyqueuosine reductase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (25/25);  total (27/27)
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

ACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCAGGCGAATAACTCAATGAGTTCC  >  NC_000913/4392750‑4393240                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    acaGCCTCAGGGTGTAATCTGAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTGTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCTCTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCATTTCAACt                                                                                                                                                                                                                                                  <  2:110385/251‑1 (MQ=255)         aGGCTGTAATCTTAATTACAAACTAACTTTTCAACGTGTATAGTAATTTTTTAACCTACTGATAGATAAGTTAAAAATTAAGATACCGTTTTGTTGAGGCTAACACATGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTAGTTTTATTCACAAGCTGGGGACGAATCATGCGTCACGCGGTATGCCTTTTTCAATCTACCGCACCAGACTCTGTTTTTTCGTTAATTGCACTTCAACTAGGCACg                                                                                                                                                                                                                                           <  2:290535/250‑1 (MQ=255)           gCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACCCAt                                                                                                                                                                                                                                        >  2:458873/1‑251 (MQ=255)                 aTCTTAATTTCAAGGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGATGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGtc                                                                                                                                                                                                                                  >  2:527611/1‑251 (MQ=255)                     tAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGGCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGtct                                                                                                                                                                                                                                 <  2:314404/248‑1 (MQ=255)                           cAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTAttt                                                                                                                                                                                                                         <  1:194259/250‑1 (MQ=255)                                       gCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTAGTCTCTCTATTTGCTGCGCAATCGc                                                                                                                                                                                                            <  1:453449/251‑1 (MQ=255)                                            gTGAATAGGAGTTTTTTAACCTAGTGAAAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATATCTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTCCAATCCCCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCCACTATGCGCGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACg                                                                                                                                                                                                       <  2:264716/251‑1 (MQ=255)                                                    tATTTTTTTAACCTATTGATAGATAAGTTAAAACTTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGc                                                                                                                                                                                               >  2:163083/1‑251 (MQ=255)                                                          tttAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCg                                                                                                                                                                                         >  1:315084/1‑251 (MQ=255)                                                                      atagataAGTTAAACATTCAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGCGAAGTGGGTGc                                                                                                                                                                             >  1:433646/1‑251 (MQ=255)                                                                           taAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTTCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCAcc                                                                                                                                                                        >  2:9011/1‑251 (MQ=255)                                                                            aaGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCt                                                                                                                                                                       >  1:63354/1‑251 (MQ=255)                                                                                        aaGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGCGAAGTGGGTGCTCACCTTTACGACTTTc                                                                                                                                                            >  1:247310/1‑250 (MQ=255)                                                                                        aaGAGTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACGTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTc                                                                                                                                                            <  2:247310/250‑1 (MQ=255)                                                                                          gATTCCGTTTTGTCGAGTCGATCACTTGACGCTTTTATTCGACTTGAATTGCAATGCGTTTTTATTTTTATTCACATGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCGATTGCACTTCAACTACGCACGCATTTCGTCTGTCTATTTGCTGCGCCATCGCCCACGCTATGTGCTCCTCGAGAAGTGGATGCTCACCTTTCCGACTTTcccg                                                                                                                                                         >  2:341780/1‑249 (MQ=255)                                                                                          gATTCCGTGTTGTCGAGTCGAACGTTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACTCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTGCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCGATGGGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAg                                                                                                                                                         <  2:24631/251‑1 (MQ=255)                                                                                                       cgagtcgaTCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATc                                                                                                                                            <  1:391946/251‑1 (MQ=255)                                                                                                            cgaTCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTCCGACTTTCCAGCGCTGTCAAAATCGTTTc                                                                                                                                       <  1:432208/251‑1 (MQ=255)                                                                                                              aTCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCAt                                                                                                                                     <  2:315084/251‑1 (MQ=255)                                                                                                                     gACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAg                                                                                                                               <  1:458873/250‑1 (MQ=255)                                                                                                                           tttATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCAtt                                                                                                                        >  2:273674/1‑251 (MQ=255)                                                                                                                                       gAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTAc                                                                                                            >  1:117343/1‑251 (MQ=255)                                                                                                                                            gCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCCTTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTGACCTTTACGACTTTCCAGCGCTGTCACAATCGTTTCATCCCCAGGGGCATTGCCTAATGCTACGGCg                                                                                                        >  2:137129/1‑250 (MQ=255)                                                                                                                                            gCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTGGCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCg                                                                                                        <  1:231491/250‑1 (MQ=255)                                                                                                                                                   gTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACg                                                                                                <  1:341780/251‑1 (MQ=255)                                                                                                                                                        tatttttattCACAAGCTGTGGATGAATCATGCGTCACGCGGTTACCCTTTTTCAATCCCCCGAACCATACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCAAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGc                                                                                            <  2:435677/250‑1 (MQ=255)                                                                                                                                                               attCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGc                                                                                      >  1:183649/1‑249 (MQ=255)                                                                                                                                                               attCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGATCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGc                                                                                      >  1:530771/1‑249 (MQ=255)                                                                                                                                                                    caAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTATGCAGCCAACGCAGGTGAc                                                                               >  1:232799/1‑251 (MQ=255)                                                                                                                                                                     aaGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGAcc                                                                              >  1:45412/1‑251 (MQ=255)                                                                                                                                                                          gtgGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATa                                                                          <  1:480338/250‑1 (MQ=255)                                                                                                                                                                            ggATGAATCAGGCGTAACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTAGGTGCTCATCGAGAAGTGGGTGCTCGCCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCATCCAACGCAGGTGACCAATacg                                                                        <  2:527512/250‑1 (MQ=255)                                                                                                                                                                             gATGAATCAGTCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGACGTGGGTGCTCACCTTTACGACTTTCCAGCGCAGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATacgac                                                                      >  1:192042/1‑251 (MQ=255)                                                                                                                                                                                gAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGAAGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGCAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTACTGCTAAGGCGATATTACGCAGCCAACGCAGGTGACCAATacgacgc                                                                    >  2:216694/1‑249 (MQ=255)                                                                                                                                                                                        cGTCACGTGGTAACCCTTTTTCAATCTCCCGAACCAGACGCTGTTTTTTAGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGTGGCAGTGCTTAATGCTACGGCGATATTACGCCGCCAACGCAGGTGACCAATACGACGAATCGCCGAt                                                           <  2:140309/251‑1 (MQ=255)                                                                                                                                                                                         gTCACGCGGCAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCAAACGCAGGTGACCAATACGACGAATCGCCGATc                                                          >  1:500642/1‑251 (MQ=255)                                                                                                                                                                                                 gTAACACTTTTTCAATCACACGAACCAGACGCTGTTTTTTCGACAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGACCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGCTATTACGCAGCCAACGCAGGGGACCAATACGACGAATCGCCGATCCTTCCgcg                                                  >  2:442859/1‑249 (MQ=255)                                                                                                                                                                                                   aaaccTTGTTCAATCACCCGAACCAGACGTGGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATGTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGa                                                 <  2:207824/247‑1 (MQ=255)                                                                                                                                                                                                                tCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCCTTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGAAGAATCGCCGATCCTTCCGTGACTTTTAAAACCTtc                                   >  2:43520/1‑251 (MQ=255)                                                                                                                                                                                                                 cACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTtc                                   <  1:9011/250‑1 (MQ=255)                                                                                                                                                                                                                  aCCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGCTATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTtctt                                 >  2:183955/1‑250 (MQ=255)                                                                                                                                                                                                                     cGAACCAGTCGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTc                              <  2:63354/251‑1 (MQ=255)                                                                                                                                                                                                                      gAACCAGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACACATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATAGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGACGATCCTTCCGTGACTTTTAAAAACTTCTCTTCg                             >  1:476544/1‑251 (MQ=255)                                                                                                                                                                                                                           aGACGCTGTTTTTTCGGCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCa                        <  1:311955/251‑1 (MQ=255)                                                                                                                                                                                                                           aGACGCTGTTTTTTCGGCAATTCCACTTCAACTATGCACGCATTTCGTATATCTATTTGCTGAGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTACAGCGCTGTCAAAATCGTTTCATCCCACGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGCCCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTccc                        >  2:53665/1‑250 (MQ=255)                                                                                                                                                                                                                             ccacTGTTTTTTCGGCAAATGCTCTACAACTATGCACGCATTTCGTCTCTCTAGTTGCTGCCCAATCGCCCATGCTATGTGCTCATCGAGAAGTGGGTGCTCGCCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCAgg                      <  2:192042/248‑1 (MQ=255)                                                                                                                                                                                                                                    tttttCGGCAATTGCGCTTCGACTATGCACGCATTTCGTCTATCTATTTGCTGTGCAATCGCCCACGCTATCTGCTCATCGCGAAGAGGGTGCTCACCTTTACGACTTTCCAGCGGTGTCAAAATCGTTTCATCCCAAAGGGCATTGCATAATGCTACGGCGATCTTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTACGTGACTTTTAAAAAATTCTATTCGCACCAGGCgtctc                >  2:3718/1‑246 (MQ=255)                                                                                                                                                                                                                                    tttttCGGCAATTGCAAATCAACTATGCACGCATTTCGTCTCTCTCTTTGCTGCGCAATCGCCCACGCTATGTGATCATCGAGAAGTGGGTGCTCACCTTTACGACTTTACAGCGCTGTCAAAATCGTTTCATACCAAGGGGCATTGCCTAATGCTCCGGCGATATTACGCAGCCAAAGCAGGTGACCCATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAATCTTCTCTTCGATCAAGGCGAATa                >  2:257581/1‑250 (MQ=255)                                                                                                                                                                                                                                           gCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCAGGCGAATAACTCAATg        >  1:275982/1‑251 (MQ=255)                                                                                                                                                                                                                                           gCAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCAGGCGAATAACTCAATg        <  1:452451/251‑1 (MQ=255)                                                                                                                                                                                                                                            cAATTGCCCTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCCCGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCAGGCGAATAACTCAATGa       <  2:500642/251‑1 (MQ=255)                                                                                                                                                                                                                                            cAATTGCACTTCAACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCAGGCGAATAACTCAATGa       <  1:182117/251‑1 (MQ=255)                                                                                                                                                                                                                                                 gCACTTCAACTATGCTCGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCAGGCGAATAACTCAATGAGTTcc  <  1:398243/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                    ACAGCCTCAGGCTGTAATCTTAATTTCAAAGAAACTTCGCACGGTGAATAGTATTTTTTTAACCTATTGATAGATAAGTTAAAAATTAAGATTCCGTTTTGTCGAGTCGATCACTTGACGATTTTATTCGTCTTGAATTGCAATGCGTTTTTATTTTTATTCACAAGCTGTGGATGAATCAGGCGTCACGCGGTAACCCTTTTTCAATCACCCGAACCAGACGCTGTTTTTTCGGCAATTGCACTTCGACTATGCACGCATTTCGTCTCTCTATTTGCTGCGCAATCGCCCACGCTATGTGCTCATCGAGAAGTGGGTGCTCACCTTTACGACTTTCCAGCGCTGTCAAAATCGTTTCATCCCAAGGGGCATTGCCTAATGCTACGGCGATATTACGCAGCCAACGCAGGTGACCAATACGACGAATCGCCGATCCTTCCGTGACTTTTAAAAACTTCTCTTCGCTCCAGGCGAATAACTCAATGAGTTCC  >  NC_000913/4392750‑4393240

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