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	90.5%	104.4 / 4.8	42	F172F (TTC→TTT)	queG	epoxyqueuosine reductase
Reads supporting (aligned to +/- strand):  ref base G (0/4);  new base A (17/21);  total (17/25)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01

ACGGGCAGATAAGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGTTTTTCAGCGT  >  NC_000913/4393305‑4393794                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   aaggtcaaaTAAGTTGCCAGTCAACTCAAACGTAAATTCGTGTTCCTAATAACGGTCGCAACTCCTCCGGGATCGCCCCTTCAAGTTCGAGGGTGAGATAAGGAATACAGCGGCGCGCATCGACGGTATATGGCTCGACGATGGCCCCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCCTCCTTCCACGACTGTTTGAAACACTGGCAGCGGAATATCGACCAGCAATTCGCGTAAAaaaaag                                                                                                                                                                                                                                                 <  2:2193/243‑1 (MQ=255)  cGGGCAGATAAGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGAGCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGGATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCGCGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAaaaaga                                                                                                                                                                                                                                                <  1:348573/251‑1 (MQ=255)       aGATAAGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGTCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAgaacga                                                                                                                                                                                                                                            >  1:442272/1‑250 (MQ=255)         aTAAGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAAAACGAAcc                                                                                                                                                                                                                                         >  2:386198/1‑251 (MQ=255)           aaGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCg                                                                                                                                                                                                                                        <  1:386198/250‑1 (MQ=255)           aaGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGGCGCAACCCTTCCGTGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGAGACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGATGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCg                                                                                                                                                                                                                                        <  2:117343/250‑1 (MQ=255)                   aGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTcgcg                                                                                                                                                                                                                                >  1:204208/1‑250 (MQ=255)                    gTCATCGCAACCGTAAATACGTTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGAt                                                                                                                                                                                                                              <  2:265303/251‑1 (MQ=255)                     ccATCGCAACCGTAAATACGGTTTCCCATTATCGGTCGCAACGCTTCCGGGATCGTCCCTTCAGGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCGCGACTGGTTGATCCCCGGGCAGCGGAATATCGAACAGCAATGCGCCTAAAAAAAAGAACGAACCGGCCTCGCGAtt                                                                                                                                                                                                                             <  2:3445/250‑1 (MQ=255)                      cATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACAATGGCACCGGTCGGGCAAATCGTCATACAAGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCAACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTg                                                                                                                                                                                                                            <  1:133712/251‑1 (MQ=255)                                gTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAACAAAGAACGAACCGGCCTCGCGCTTGAGGATAagtg                                                                                                                                                                                                                  >  2:262790/1‑251 (MQ=255)                                gTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGACTCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCATCTCGATTGAGGATAaggg                                                                                                                                                                                                                  >  2:66417/1‑249 (MQ=255)                                  aaaTACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCACTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAATAACGAACCGGACTCGCGATTGAGGATAagtgag                                                                                                                                                                                                                <  1:146891/251‑1 (MQ=255)                                       cGGTTGCCCATTAATGGTCGCAACGCGTCCGTGACAGCCTTGTCAAGTTCGATGGGGAGATAAGATATACAGCGGCGAGGATCGACGGTATATGGGGCGACGATGGAACCGGTCGGGCAAATCGTAATACAGGCCACGCATTTGCCGCCTCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCTGCCTCGCGATTGAGGATAATTGAGTGcat                                                                                                                                                                                                           <  2:530210/251‑3 (MQ=255)                                                            aaCTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCGGTCCAGCCGAGCCCAGc                                                                                                                                                                                      <  2:114653/251‑1 (MQ=255)                                                                    cGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTc                                                                                                                                                                                >  2:292827/1‑249 (MQ=255)                                                                          cGCCCCTTCAAGTTCGATGGTGAGAAAAGAGATACAGCGGCGAGCATCTACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACg                                                                                                                                                                        <  2:222617/251‑1 (MQ=255)                                                                                    aGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCCCGCATTTGCCGCATCCTTCCTCGACTGGTTGAGCCACGGGCACCGGAATCTCGACCAGCACTCCGCCTAAAAAAAAGAACGAACCGGCCCCGCGATTGCGGATCAGTGAGTGCTTACATGTCCAGCCGACCCCAGCTTTTTCAGTTAACGAgcgaccgcg                                                                                                                                                              >  2:172195/1‑245 (MQ=255)                                                                                          aTGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAgg                                                                                                                                                        <  1:399742/251‑1 (MQ=255)                                                                                               gagaTAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAg                                                                                                                                                   <  2:204208/251‑1 (MQ=255)                                                                                                    aaGAGGTACATCGGCGGGCAGCGACGGTATAGGGCGCGACGATGGCACCGGTCCGGCAAAGCGTTATCCAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCGTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGaagcg                                                                                                                                              <  2:176345/251‑4 (MQ=255)                                                                                                                     gCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAAtt                                                                                                                             <  1:216694/251‑1 (MQ=255)                                                                                                                           aCGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGa                                                                                                                       <  1:292827/251‑1 (MQ=255)                                                                                                                                    ggCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATg                                                                                                              <  1:233429/251‑1 (MQ=255)                                                                                                                                      ttCGACGACGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCAGTTGCCGCATTCTTCCTCGACTGGGTGATCCCCGGGCAGCGGAATATCGACCATCAAGCCGCCTAAAAAAAAGAACTAACCGGCCTCGCTATCGAGGAGAAGTGAGTGCTTACCGGTCCAGCCGAGCCCAGCTTTTGCAGCTAACGTGTGCCCGAGAATAGGCGCAGAAGCGACAAATGGTCGAAATTTCAGCGCAACACAATGt                                                                                                             <  2:65726/249‑1 (MQ=255)                                                                                                                                         gacgaTGGCACCGGACGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGACGCTCGAGAATAGGCCCAGAATCGACAAACGGTCTAAAATTCAGCGAAAAACAAtttttct                                                                                                         >  2:415776/1‑245 (MQ=255)                                                                                                                                               ggCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAatc                                                                                                    >  1:184053/1‑250 (MQ=255)                                                                                                                                               ggCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGAGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAatc                                                                                                    <  1:43520/250‑1 (MQ=255)                                                                                                                                                        cggtcAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGAGCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCa                                                                                          <  2:209277/247‑1 (MQ=255)                                                                                                                                                                     gtaccggCCACGCATTTGCCGCCTCCTTCCTCGACGGGTGGATCCACGGGCAGCGTAAGATCGCCCAGCAATTCGCATAAAAAAAAGAACTATCCGGCCTCGCGATGGAGAATAAGAGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGGCTAAAATTCAGCGAAAGACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCg                                                                             <  2:330239/246‑1 (MQ=255)                                                                                                                                                                     aTACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCg                                                                             >  1:434033/1‑251 (MQ=255)                                                                                                                                                                        cAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGACCGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACACACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGtt                                                                          >  2:519999/1‑251 (MQ=255)                                                                                                                                                                         aGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTg                                                                         >  2:111857/1‑251 (MQ=255)                                                                                                                                                                                gCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTTAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGAGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTAGCCCAGCTTTTTGAGTCGGTTGCGCAGa                                                                   >  1:75477/1‑250 (MQ=255)                                                                                                                                                                                    ttgccggtgccttcctGGAATGGGTGAGCCACGGGCAGCGGAATATCTACCAGCAATTCGCAGAAAAAAAAGAACGACCCGGCCTCGCGAGTGAGGATAAGTGAGGGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGtt                                                               <  2:542070/241‑1 (MQ=255)                                                                                                                                                                                          cATCCTTCCTCGACTGGTTGCTCCACGGGCAGCGGAATATCGACCAGCAATTCGCATAATAAATAGAACGAACCGGACTCGCGCTAGAGGATAAGTGATAGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCCGCTAACGGGCGCTCGAGAATAGGCGCAGAAACGACAAACGGTCTAAAATTCAGCGCAACACAATGTTGCTGAATCATCTAGCACAGCATTTTGAgacgaatgtgcaaaaggttgta                                                           >  2:387238/1‑227 (MQ=255)                                                                                                                                                                                                 cctCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACgg                                                 <  1:491026/251‑1 (MQ=255)                                                                                                                                                                                                  ctCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGc                                                >  1:515906/1‑251 (MQ=255)                                                                                                                                                                                                    cGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGccc                                              <  1:262790/251‑1 (MQ=255)                                                                                                                                                                                                           ttGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTTATAGTCACGGCCCAGCGCAt                                       >  1:301962/1‑251 (MQ=255)                                                                                                                                                                                                            tGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCAt                                       <  1:3718/250‑1 (MQ=255)                                                                                                                                                                                                                     ggCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTACAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGGTTGTGATaagaacggacaagcgcataaaagct                                >  2:253904/1‑224 (MQ=255)                                                                                                                                                                                                                          cGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGc                        <  1:247330/251‑1 (MQ=255)                                                                                                                                                                                                                                                ccTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGATCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGTTTTTCAGCGt  <  2:442272/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   ACGGGCAGATAAGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGTTTTTCAGCGT  >  NC_000913/4393305‑4393794

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