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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,393,552	0	G	A	96.1%	146.7 / ‑3.2	51	F172F (TTC→TTT)	queG	epoxyqueuosine reductase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (27/22);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.47e-01

ATAAGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGTTTTTCAGCGTGCT  >  NC_000913/4393313‑4393797                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       aTAAGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAAcc                                                                                                                                                                                                                                             >  1:271507/1‑251 (MQ=255)   aaGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCg                                                                                                                                                                                                                                            <  1:115035/250‑1 (MQ=255)   aaGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCAGCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCg                                                                                                                                                                                                                                            <  1:271354/250‑1 (MQ=255)     gCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGATCGTCGATGGCACCGGTCGGGCAAATAGTCATAAAAGCCACGCATTTACCGCATCCTTCCTCGACTGGTTGATCCACGTGCAGCGGAATATCGACCCACAATTCGCCTAAAAAAAAGAACGAACCGGcc                                                                                                                                                                                                                                         >  2:292181/1‑251 (MQ=255)               aTCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTg                                                                                                                                                                                                                                >  1:2407/1‑250 (MQ=255)                    aaCCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGAt                                                                                                                                                                                                                           >  2:108571/1‑250 (MQ=255)                    aaCCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGAGATTGAGGAt                                                                                                                                                                                                                           >  1:10750/1‑250 (MQ=255)                          aaaTACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCGCGACGATGGCACCGGCCGGGAACATCGTCATACAGGCCACGCATTTTCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCATCAATTCGACTAAAAAAAAGACCGAACCGGCCTCGCGATTTAGGATCagtg                                                                                                                                                                                                                      >  2:74941/1‑249 (MQ=255)                                       ccttAACGGTCGCAACTCTTCCGGGCTCGCCCCTACAAGTTCGATGGTGAGATAAGAGATACAGCGGAGAGCATCGACGGTATATGGCTCGAGGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGACTTTGCCGAATCCTTCCTCGACTGGTTGATCCACGGGCCGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGCACCGGACTCGAGATTGGGGaggagggggtgctatacagtac                                                                                                                                                                                                       >  2:324265/3‑230 (MQ=255)                                       cATTAACGGTAGCAACTCTTACGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGCTACAGCGGCGAGCATCGACGGTATATGGCTCGCAGATGGCCCCGGTCGGGCAAATCGTCATACACGCCACGCATTTGCCGCATCCTTCATCCACTGGTTGCTCCACGGGCAGCGGCATATCGACCAGAAATACGCCTAAACAAAAGCACGAACCGGCCTCGCAATTgcagatacgggagttattacctatcc                                                                                                                                                                                                       >  2:69870/1‑226 (MQ=255)                                          tAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAg                                                                                                                                                                                                     <  2:184179/250‑1 (MQ=255)                                             cGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGCCGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGCTTGAGGTTAAGTGAGTGCTTACCTGTCCAGCCg                                                                                                                                                                                                  >  2:61980/1‑250 (MQ=255)                                              ggTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCg                                                                                                                                                                                                  >  1:86175/1‑249 (MQ=255)                                                     aCTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCGTACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAACAAAGAACGAACCGGCCTCGCGATTGGGGATAAGGGAGTGCTTACCTGACCaaaggagaccagct                                                                                                                                                                                         >  2:378034/1‑238 (MQ=255)                                                          tCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTc                                                                                                                                                                                   >  1:436749/1‑251 (MQ=255)                                                                         gcAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCAGCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGgcgcgc                                                                                                                                                                    <  2:140439/250‑3 (MQ=255)                                                                                 gATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAg                                                                                                                                                            <  1:394271/251‑1 (MQ=255)                                                                                 gATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCACTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAg                                                                                                                                                            <  1:88264/251‑1 (MQ=255)                                                                                         gaTAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGTAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGa                                                                                                                                                     >  1:97776/1‑250 (MQ=255)                                                                                                gaTACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACa                                                                                                                                              >  1:41970/1‑250 (MQ=255)                                                                                                aaTACAACGGAGACCATCGAAGGAATATGGCTCGACGCTGGAACCGGTCCAGAAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGCCTGGTTGCTCCAAAGGCAGCGGAATATAGACCAGCAATCAGACTAAAAAACAGAACGAACCGGCATCGCGACTGAGGATAAGGGAGTGATTACCTGTACAGCCGAGCC‑AGACTATTTCAGCTAAAGGGCGCTCGAGAACAGGCGCCGaaacggca                                                                                                                                              >  2:209122/2‑244 (MQ=255)                                                                                                  aaCAGTGGCGAGCATCGACGTTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGGCATGCAGTCCACGCATTTGCCGCCTCCTTCCCCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACaagc                                                                                                                                           <  2:24845/250‑3 (MQ=255)                                                                                                          cgAGCAGCGACGGTATAAGGCTCGACGATGGCACCGTTCGGGCAAAACGGCATACAGGCCACGCTTGTGCTGCAACCTTCCTCGACTGGTTGATCCACGGGCAGTGGAATATCGACAAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCACGCGATTGATGATAAGTGAGTGCTTACCTGTCCAGCCTAGCCCAG‑CTTGTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTaaa                                                                                                                                   <  2:374400/251‑1 (MQ=255)                                                                                                              cctCGACGGTACAGGGCTCGACTATGGCACCGGTCGGGCAAATCGTCTGACAGGCCCCGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCATCAATTCGCCTACAAAAAAGAACGACCCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCGGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTc                                                                                                                               <  2:11256/249‑1 (MQ=255)                                                                                                                  gACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTTCCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATGCAGCg                                                                                                                           <  1:324265/251‑1 (MQ=255)                                                                                                                     gggatatGGCTCGAGGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACTCATTTGTTGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGTCCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTTAGCTAACTGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGaaa                                                                                                                        <  1:464499/248‑1 (MQ=255)                                                                                                                     ggTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCCCGCGATTGAGGATAAGTGCGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGATAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAACTTCAGCGaa                                                                                                                         >  1:282410/1‑250 (MQ=255)                                                                                                                      gTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAac                                                                                                                       <  1:61980/251‑1 (MQ=255)                                                                                                                            ggCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATg                                                                                                                 >  1:204826/1‑251 (MQ=255)                                                                                                                                     agggCACCGGTCGGGCAAATCGTTATACAGGCCACGCAATTGCCGCATCCTTCCTCGACGGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAat                                                                                                        <  2:441810/249‑1 (MQ=255)                                                                                                                                       ggCACCGGGCGGGCAAAGCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAatc                                                                                                       <  2:271507/250‑1 (MQ=255)                                                                                                                                                       aTCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCttttt                                                                                      <  1:159638/251‑1 (MQ=255)                                                                                                                                                                    ccACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGCGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTgggc                                                                         >  2:74728/1‑248 (MQ=255)                                                                                                                                                                     cACGCATTTGCCTCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATAGCGACCAGCAATTCGCCTAAAAAAAATAACGAACCGGCCTCTCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCGCGAGAATAGGCGCAGAATCGTCAAACGGTCTAAAATTCATCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCa                                                                        <  1:209122/251‑1 (MQ=255)                                                                                                                                                                                cGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTtgtg                                                              >  2:467430/1‑250 (MQ=255)                                                                                                                                                                                       ttcctCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGCGTCGGGTGCGCAGa                                                                      >  2:38898/1‑235 (MQ=255)                                                                                                                                                                                       ttcctCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGa                                                                      <  1:38898/235‑1 (MQ=255)                                                                                                                                                                                        tcctCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCt                                                                                                            <  1:280128/196‑1 (MQ=255)                                                                                                                                                                                        tcctCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCAAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTAAGAGACACACAATGTTGCt                                                                                                            >  2:280128/1‑196 (MQ=255)                                                                                                                                                                                                   ttGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGCGTCGGTTGCGCAGCAGTTTGTGATAGTCACGGCCCCGCGCAt                                          >  1:214672/1‑251 (MQ=255)                                                                                                                                                                                                      aTCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGACACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGGCGGTTTCGCAGTAGTTTGTtatatgccacgccccgcgcctaac                                       >  2:146538/1‑231 (MQ=255)                                                                                                                                                                                                        ccACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACg                                      <  1:253694/250‑1 (MQ=255)                                                                                                                                                                                                                 gCGGAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCTACAAACGGTCTAAAATTCAGCGAACCACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATCGTCACGGCCCAGCGCATAACGGCTAACAt                              >  1:432978/1‑249 (MQ=255)                                                                                                                                                                                                                   ggAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGa                                                                                                                                                <  2:219019/133‑1 (MQ=255)                                                                                                                                                                                                                   ggAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGa                                                                                                                                                >  1:219019/1‑133 (MQ=255)                                                                                                                                                                                                                   ggAATATCGACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGGTAGTCACGGCCCCGCGCATAACGGATAACATAGcc                          >  2:77748/1‑251 (MQ=255)                                                                                                                                                                                                                           gACCAGCAATTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTgg                  <  1:108571/251‑1 (MQ=255)                                                                                                                                                                                                                                gCAATCCGCCTAAAAAAAAGAACGAACCGGCCTCGAGATTGAGGATATGTGAGGGCTTACCTGTCCATCCGAGCCCAT‑CGTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCGAAAATGCAGCGAAACACAATGTTTCTGAATCATCTCGTCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGACCCGCGCATAACTGCTAACATAGCCGAGTTTGGGGtt              <  2:72471/250‑1 (MQ=255)                                                                                                                                                                                                                                  aaTTCGCCTAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGttttt           <  2:41970/251‑1 (MQ=255)                                                                                                                                                                                                                                     tCGCCTAAAAAGAAGAACGAACCGGCCTCGAGATTGAGGATAAGTGAGTGCTTACATGTCCAGCCGAGACCAG‑CTTTTTCAGATAACGGGCGCTCGCGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCCGCGTAACACAACGATGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGCTTGTGATAGTCACGGCCAAGCGCATAAAGGCTAACATAGCCGAGATTGGGGTTTTTCa         >  1:86481/1‑250 (MQ=255)                                                                                                                                                                                                                                          tAAAAAAAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGAGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTATAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAAAATAGCCGAGTTTGGGGTTTTTCAGCGTGc   >  1:47743/1‑251 (MQ=255)                                                                                                                                                                                                                                           aaaaaaaaGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAACCGAGTTTGGGGGTTTTCAGCGTGCt  >  1:432042/1‑251 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                       ATAAGCTGGCAGTCATCGCAACCGTAAATACGGTTTCCCATTAACGGTCGCAACTCTTCCGGGATCGCCCCTTCAAGTTCGATGGTGAGATAAGAGATACAGCGGCGAGCATCGACGGTATATGGCTCGACGATGGCACCGGTCGGGCAAATCGTCATACAGGCCACGCATTTGCCGCATCCTTCCTCGACTGGTTGATCCACGGGCAGCGGAATATCGACCAGCAATTCGCCTAAAAAGAAGAACGAACCGGCCTCGCGATTGAGGATAAGTGAGTGCTTACCTGTCCAGCCGAGCCCAG‑CTTTTTCAGCTAACGGGCGCTCGAGAATAGGCGCAGAATCGACAAACGGTCTAAAATTCAGCGAAACACAATGTTGCTGAATCATCTCGCCCAGCTTTTTGAGTCGGTTGCGCAGAAGTTTGTGATAGTCACGGCCCAGCGCATAACGGCTAACATAGCCGAGTTTGGGGTTTTTCAGCGTGCT  >  NC_000913/4393313‑4393797

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