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	94.8%	173.3 / 1.0	58	F172F (TTC→TTT)	queG	epoxyqueuosine reductase
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (25/30);  total (26/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 = 2.29e-01

GGCAGATAAGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGTTTTTCAGCGTGCTG  >  NC_000913/4393308‑4393798                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        ggCAGTTAAGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATGAACGGTCGCAACTCTTCCGGGATCGCCCCGTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCTCCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑gaac                                                                                                                                                                                                                                                  <  2:643048/251‑1 (MQ=255)      aTAAGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAAcc                                                                                                                                                                                                                                             >  1:131858/1‑251 (MQ=255)        aaGCTGGCAGTCATCGCAACCGTAAATACGGGTTCCCAATAACGGTCGCAACGCTGCCGGGATCGCCCCGTCAAGGTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCAGTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCg                                                                                                                                                                                                                                            <  2:150772/250‑1 (MQ=255)         aGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGc                                                                                                                                                                                                                                          <  1:678329/251‑1 (MQ=255)          gCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGcc                                                                                                                                                                                                                                         >  1:221987/1‑251 (MQ=255)             ggCAGTCATCTCAACCGTAAAAACGGGTTCTCATTAACGGGCGCAACTCGTCCGGGAACGCCCCTTCAAGTGCGATGGTGATATAAGAGATACAGCGGCGAGCATCGACGGTAGATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATGTGTCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATAACGACCAGCAATCCGCCTAAAAAAAA‑TAACGAACCGGctt                                                                                                                                                                                                                                        <  2:97935/249‑3 (MQ=255)             ggCAGTCATCGCAACCGTAAATACGGTTTCCCATTACCGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCCTCGACGGTATATGGCTCGACGAAGGCACCGGTCGGGCAAATCGTCAAAAAGGCCACGCATTTGCCGCATCCTTACTCGACTGGTTGTTCCACGGGCAGCGGAATATTGACCAGAAATTCGCCTAAACAAAA‑GAACGAAACGGcagcg                                                                                                                                                                                                                                      >  2:529756/1‑247 (MQ=255)                      cGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAgg                                                                                                                                                                                                                             <  2:359731/251‑1 (MQ=255)                        caaccggaaaatcGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGag                                                                                                                                                                                                                           <  2:167589/239‑2 (MQ=255)                        cAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGCAGATGGCACCGGTCGGGCAAATCGTCATCCAGGCCACGCATTTGCCGCCTCCTTCCTCGACTGCTTGATCCACGGGCAGCGGCAGATCGACCAGCACTTCGCCTAAAAAAAA‑GAACGAAACGCCCACGCGCTTGAGGa                                                                                                                                                                                                                            >  2:613267/1‑250 (MQ=255)                           ccGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGTGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGACTCGCGCTTGAGGATa                                                                                                                                                                                                                          >  1:586163/1‑249 (MQ=255)                           ccGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAag                                                                                                                                                                                                                        <  1:232355/251‑1 (MQ=255)                               aaatttgggttCCCAGTAACGGCCGCAACCCTTCCGGGATCGCCCCTGCAATTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGTCGAGGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGGCCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAagtga                                                                                                                                                                                                                     <  2:553360/241‑1 (MQ=255)                                     ggTTGCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGGCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTa                                                                                                                                                                                                              <  2:548848/251‑1 (MQ=255)                                      gTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCAACGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCCTCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCACTTCGCCTAAAAAAAA‑GAACAAACCGGCCTCGCGATTATGGATAAGTGCGTGCTTAc                                                                                                                                                                                                             >  2:98234/1‑251 (MQ=255)                                        ttCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGCACCGGCCTCGCGATTGCGGCTAAGTGAGTGCTTACCt                                                                                                                                                                                                           >  1:31631/1‑251 (MQ=255)                                               tAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAA‑GAACGAACCGGCCTc                                                                                                                                                                                                                                       >  1:477277/1‑216 (MQ=255)                                               tAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGc                                                                                                                                                                                                    >  1:397816/1‑251 (MQ=255)                                                   ggTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCg                                                                                                                                                                                                  >  1:135210/1‑249 (MQ=255)                                                          aCTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCt                                                                                                                                                                                         <  1:157105/251‑1 (MQ=255)                                                                       cGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACg                                                                                                                                                                            >  1:196502/1‑251 (MQ=255)                                                                                   tgtgATGGTGAGATAATGGTAACAGCGTCGAGCATCGGCGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCgaga                                                                                                                                                                 <  2:569610/247‑1 (MQ=255)                                                                                       aTGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAgg                                                                                                                                                            <  1:64170/251‑1 (MQ=255)                                                                                        gggTGAGATAAGTGATTCAGCGGCGAGCATCGCCGGTATATGGCTCGACGATGGCACCGGTCGGGCAAACCGTCATACAGGCCACGCATTTTCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGgc                                                                                                                                                           <  2:454096/250‑1 (MQ=255)                                                                                         ggTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGgc                                                                                                                                                           <  1:385686/250‑1 (MQ=255)                                                                                         ggTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGgcg                                                                                                                                                          >  1:572349/1‑251 (MQ=255)                                                                                          gTGAGATAAGAGATAAAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACAGGCCTCGCGGTTGAGGATAAGTGGGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAACTAACGGGAGCGCGAGGAAAGgcgc                                                                                                                                                         >  2:659865/1‑251 (MQ=255)                                                                                             agaTAATAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGAGGGAACCAGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCCGCGGAATATCGACCAGCAATTAGCCTAAAAAAAA‑GAAAGAACCGGCCTCGAGCTTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCCGCTTTTTCTGGAAACGGGTGCTCGAGaaaaagcaacgg                                                                                                                                                      >  2:480577/1‑241 (MQ=255)                                                                                              gaTAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGaa                                                                                                                                                     <  1:529463/251‑1 (MQ=255)                                                                                                  agagATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGGTGATCCACGGGCAGCGGAATATCGACCAGCAAGTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGa                                                                                                                                                 <  1:211738/251‑1 (MQ=255)                                                                                                    agaTACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACa                                                                                                                                               <  1:209418/251‑1 (MQ=255)                                                                                                     gaTACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACaa                                                                                                                                              <  1:198666/251‑1 (MQ=255)                                                                                                      aTACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCATGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACaaa                                                                                                                                             <  1:353469/251‑1 (MQ=255)                                                                                                         cAGCGGCGAGCATCGGCGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGTCCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACgg                                                                                                                                          <  1:457597/251‑1 (MQ=255)                                                                                                               cgAGCATCGGCGGTATATGGCTCGACGCTGGCACCGGTCGGGCAAATCGTCATACAGGCCCCGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCGGCAATTCGCTTAAGAAAAT‑GAACGAACCGGCCTCGCGATTGAGGGTAGGTGAGGGCTTCCCTGTCCTGCCGTGCCCAGCTTTTTCAGCTAACGTGTGCTAGAGGATTGGCGCAGAATCGACAAACGGTCTaa                                                                                                                                     >  1:401642/1‑250 (MQ=255)                                                                                                                gagtaTCGAGGGTAAATGGGTGGACGATGGCATCGGTCGGGCAAAGCGGCAGACAGGCAACGCAGGTGGCGCATCCTTCCGCGACGGGTTGATCCACGGGAAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTaaaa                                                                                                                                   <  2:235316/247‑1 (MQ=255)                                                                                                                    aTCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGAACCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCa                                                                                                                               <  2:167535/251‑1 (MQ=255)                                                                                                                      cGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAg                                                                                                                              <  1:157347/250‑1 (MQ=255)                                                                                                                       gACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCg                                                                                                                            >  1:247332/1‑251 (MQ=255)                                                                                                                                                          aaaTCGGCATAAAGGCCTAGCATATGCCGAATCATTCCTCGACTGGTTGCTCCACAGGCAGAGGAATATAGACCAGCAATTCGCCTAAAAGAAA‑GAACGCACCGGACTCGCGATTGAGGATAAGTGAGAGCTTACCTGTCCAGCCGAGCCCAGCTATTTAAGCTAAGGGGCGCTCGAGAATAGGCGCAGACTCGCCAAACGGGCTACAATGCAGCGAAACACAATGTTGCTGAATCATCTCGGCTAGCttt                                                                                         >  2:98811/1‑251 (MQ=255)                                                                                                                                                           aaTCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAatca                                                                                                       <  1:464745/236‑1 (MQ=255)                                                                                                                                                           aaTCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCACAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACAAAATGTTGCTGAatcaagatcgaaagaacgt                                                                                        >  2:464745/1‑236 (MQ=255)                                                                                                                                                                   tACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAAtt                                                                                                                                 >  2:586760/1‑202 (MQ=255)                                                                                                                                                                   tACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAAtt                                                                                                                                 <  1:586760/202‑1 (MQ=255)                                                                                                                                                                     caggcaccGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGtt                                                                              <  2:141419/244‑1 (MQ=255)                                                                                                                                                                           aCGCATTTGCCGCATCCTTACTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCCGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCATTTTGAGTCGGTTgcgc                                                                          >  2:412859/1‑249 (MQ=255)                                                                                                                                                                                          ccttcctCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAAACGGCCTCGCGATTGAGGATAAGTGAGTGCTTTCCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCGCGAGAATAGGCGCATAATCGACAAACTGTCTAAAATTCAGCTAAACTCAATGTTGCTGAATCATCGCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGGTAGTc                                                         <  2:654286/251‑1 (MQ=255)                                                                                                                                                                                                     gggTTGATCCCCGGGCAGCGGAATATCGACCAGCAATTCTCCTAAAAAAAA‑GAACGAACCTGACTCGCGATTGAGGATAATTGAGTGCTTCCCTGTCCAGCCTAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATGCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAgcg                                              <  2:586163/250‑1 (MQ=255)                                                                                                                                                                                                              cACGGGCAGCGGAATCTCGAACAGCAATTCGCATAAAAAAAA‑GAAAGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGACCAGCTTTTTCAGCTAACGGGCGATCGAGAATAGGCGCAGAATCGACCAACGGTCTAAAATTAAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGGTGCGCAGAAGTTTGTAATAATCACGGCCCATCGCAtcaaggc                                     >  2:432715/1‑245 (MQ=255)                                                                                                                                                                                                                  ggCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTaa                                  >  1:330608/1‑250 (MQ=255)                                                                                                                                                                                                                     aGCGGAATATCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATa                              <  1:362286/251‑1 (MQ=255)                                                                                                                                                                                                                             atCGACCAGCAATTCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCGCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGt                       <  2:247332/250‑1 (MQ=255)                                                                                                                                                                                                                             agcGACCAGCAATTCGCCTAAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAg                        <  1:585569/248‑1 (MQ=255)                                                                                                                                                                                                                                         ttCGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGATGCGCAGAAGATTGTGATAGACACGGCCCAGCGCATAACGGCTAACATAGCCGGGTTTGGGGTTTTTc           >  2:667424/1‑250 (MQ=255)                                                                                                                                                                                                                                           cGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGACTAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCCCGGCCCa                                                 >  2:684747/1‑210 (MQ=255)                                                                                                                                                                                                                                           cGCCTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCa                                                 <  1:684747/210‑1 (MQ=255)                                                                                                                                                                                                                                             ccTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAATTTGTGATAGTCCCGGCCCAGCGCATCACCGCTAACATAGACGAGATTGGGGGTTTTCAggtg      >  2:380074/1‑248 (MQ=255)                                                                                                                                                                                                                                              cTAAAAAAAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGTTTTTCAGCGTg     >  1:474838/1‑251 (MQ=255)                                                                                                                                                                                                                                                aaaaaaaa‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGTTTTTCAGCGTg     <  1:325715/249‑1 (MQ=255)                                                                                                                                                                                                                                                 aaaaaaA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCGAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGTTTTTCAGCGTGCTg  <  2:114954/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                        GGCAGATAAGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAA‑GAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAGCTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGTTTTTCAGCGTGCTG  >  NC_000913/4393308‑4393798

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