Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,434,861	0	G	T	28.3%	61.9 / 12.4	46	P259Q (CCA→CAA)	truA	tRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base G (3/30);  new base T (13/0);  total (16/30)
Fisher's exact test for biased strand distribution p-value = 5.50e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACCAGACCGGTTTCACAGAACAAAATTAAGAACAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCGTTGCT‑GCCGCCAGCGTT‑CTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGTTG  >  NC_000913/2434611‑2435109                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                               accaGACCGGTTTCACAGAACAAAATTAAGAACAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAAGCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTg                                                                                                                                                                                                                                                                <  1:244331/251‑1 (MQ=255)         ggTTTCACAGAACAAAATTAAGAACAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTggggagtg                                                                                                                                                                                                                                                         >  2:596194/1‑244 (MQ=255)            ttCACAGAACAAAATTAAGCACAAAATGGCATACACCCAGACGCCGTACTCCGCGACCAGCTCTGACAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGACTACAGCCTGAACTGTTCCGAGCATTCTTCTCGTTAGTCCGCCAGAAATCGCGGGCCCATTGGCGGTTTTgttcgattatac                                                                                                                                                                                                                                                     >  2:571377/1‑240 (MQ=255)            ttCACAGAACAAAATTAAGAACAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGACATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAACTAGCGGGCCCATTGGCGGGTTTTGTAGATCAttc                                                                                                                                                                                                                                                     >  1:551385/1‑249 (MQ=255)                 aGAACAAAATTAAGAACAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTc                                                                                                                                                                                                                                               <  1:136882/251‑1 (MQ=255)                     cAAAATTAAGAACAAAAAGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGTATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCGGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCAGTATTCTCGTTAGTCCGCCTGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATC‑CCGGTCAggg                                                                                                                                                                                                                                           <  2:510557/251‑1 (MQ=255)                      aaaaTTAAGAACAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGt                                                                                                                                                                                                                                          <  2:231680/251‑1 (MQ=255)                                            aaaCCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCCTA‑CCCGGCCGGGGAAACCACCGCGAACAGAAACa                                                                                                                                                                                                                     >  1:425128/1‑250 (MQ=255)                                            aaaCCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAGGAACATA‑CCGGTCAGGGGAATCCAcccccgcccgaaaac                                                                                                                                                                                                                     >  1:346147/1‑237 (MQ=255)                                                 cAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGCAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTGGGAGATCATAACCGG‑CAGGGGAATCCACCCCGACCCGATACA‑Gc‑ccct                                                                                                                                                                                                             >  1:390002/1‑249 (MQ=255)                                                    aCGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGTAGATCATA‑CCCGGCAGGGGAATTCACCCCGGCCCAATAAA‑GC‑CCTT‑ccc                                                                                                                                                                                                         >  1:222840/1‑250 (MQ=255)                                                      gCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCCATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATTATA‑ACGGTCAGGGGAATCCACCCCGCCCCGAAACA‑GC‑CCCT‑ccccc                                                                                                                                                                                                       >  1:31911/1‑248 (MQ=255)                                                           aCTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATA‑CCGGGCAGGGGAATCCCCCGCCGCCAGATAAA‑AC‑CCCT‑CCGCTTTTgt                                                                                                                                                                                                  >  1:132426/1‑250 (MQ=255)                                                                   aCCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCGTtgct‑g                                                                                                                                                                                         <  2:191148/251‑1 (MQ=255)                                                                    ccAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCTCTTTTGCCGTtgct‑gc                                                                                                                                                                                        <  1:413344/251‑1 (MQ=255)                                                                       gCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATTCTA‑CCCGGCAGGGGAAACCACCCCGGACAGCTACA‑GCACCTT‑CCCCttttcccttgccgcccc                                                                                                                                                                                       >  2:450846/1‑238 (MQ=255)                                                                        cactGCCAGATGCACGTCAATGTGCAGTATAAAATCGATGAGGAAATAAAGCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATGCTCGTTAGTCCGCCAGAAATAGCGGGCCCATGGGCGGTTTTGGAAGATCATA‑CCGTTCAGGGTAATCCACCGCGACCAGAGACA‑GC‑CCTT‑CCGCTTTTGCCGTTGCT‑gccgcc                                                                                                                                                                                    <  2:573936/249‑1 (MQ=255)                                                                                      cGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCGTTGCT‑GCCGCCAGCGTT‑CTGTCCtt                                                                                                                                                                     <  1:412966/251‑1 (MQ=255)                                                                                             gtgCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATAGTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTAGTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCGTTGCT‑GCCGCCAGCGTT‑CTGTCCTTTgccgc                                                                                                                                                               <  2:115778/250‑1 (MQ=255)                                                                                               gCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTCCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCGTTGCT‑GCCGCCAGCGTT‑CTGTCCTTTGCCGCCAg                                                                                                                                                            <  2:256504/251‑1 (MQ=255)                                                                                                  ggATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCtttt                                                                                                                                                                                                    <  1:433907/210‑1 (MQ=255)                                                                                                    aTAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCGtt                                                                                                                                                                                              <  2:136525/214‑1 (MQ=255)                                                                                                      aaaaTCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATA‑CCGGTCCGGGGAATCCACCGCGACCCGATAAA‑GC‑CCTTTCCCCTTTT‑CCCTTGCT‑GCCGCCCGCGTT‑CTGTCCCTTgccgcccgccacccc                                                                                                                                                     >  1:442128/1‑242 (MQ=255)                                                                                                        aaTCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCCTA‑CCGGTCAGGGGAATACCCCGCGGCCAGATACA‑AA‑CCCT‑CCCCTTTTTCCGTTGCTTGCC‑CCACGGGT‑CTGTCCCTTTCCGCCCGCAActcttc                                                                                                                                                   >  1:433474/1‑249 (MQ=255)                                                                                                             atggggAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCTTCATCTTTCAGGCTTGAACCACGGTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGACCGCCAGAAATAGCGGGCCCATTGGCGGTTTGGGAAGATCATG‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCGTTGCT‑GCCGCCAGCGTT‑CTGTCCTTTTCCGCCAGCAACTCGGCTATcc                                                                                                                                              <  2:243633/247‑1 (MQ=255)                                                                                                             aTGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCGTTGCT‑GCCTCCAGCGTT‑CTGTCCTTTGCCGCCAGCAACTCTGCTATcc                                                                                                                                              <  2:489920/251‑1 (MQ=255)                                                                                                                 ggAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATTATA‑CCGGGCAGGGGAATCCACCGCGCCCAGATACA‑GC‑CCTTTCCCCTTTT‑CCGTTGCT‑TCCGCCCGCGTTTCTGCCTTTTGCCCCCA‑CCACTCCGCTATcccgcc                                                                                                                                          >  1:16892/1‑247 (MQ=255)                                                                                                                                  atatTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGTAGATCAAA‑CTGGGTAGGGGAATTCACCACGGCCAGAGACAATC‑CCTT‑CC‑CACTTTACGTTGCA‑CCCGCCCGCGTT‑TTGTCCTCTTCCGCCCCCAACTACGCTATCCCCCTCTCCGGCGGGTTGTggg                                                                                                                         >  1:61122/1‑251 (MQ=255)                                                                                                                                      ttGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCGTTGCT‑GCCGCCAGCGTT‑CTGTCCTTTGCCGCCAGCAACTCTg                                                                                                                                                    <  2:324446/220‑1 (MQ=255)                                                                                                                                       gtccTATTCCTTGTACTCGTCATCGTTCATGCTGTACCTACGTGGGCTACAGCCTGAAATGTTCCGAGCATGATGCTCGTGATTCCGCCAGAAAGAGCGGGCCCATGGTCGGTTTTGGAATATCATA‑CCGGTCAGTGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCTTTGCT‑GCCGCCATCGTT‑CTGTCCGTTGCCGCCAGCAACTCGGCTATCCAGCTCGCCTGCTTGTTGTGGGcccc                                                                                                                     <  2:342265/248‑4 (MQ=255)                                                                                                                                            aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATA‑CCGGTCCGGGGAATCCACCGCGACCCGATACA‑GC‑CCTT‑CCGCCTTTTCCGTTTCT‑GCCGCCCGCGTT‑CTGTCCCTTGCCGCCAGCAACTCTGCTTTCCCGCTCTCCGGCTGGGTGTGGGCCCCCACTTc                                                                                                               >  1:490137/1‑251 (MQ=255)                                                                                                                                                                ttaatgctttAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGGCCGCCAGAAAGAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCGTTGCT‑GCCGCCAGCGTT‑CTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa                                                                                               <  2:551385/238‑1 (MQ=255)                                                                                                                                                                          aaCTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCGTTGCT‑GCCGCCAGCGTT‑CTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCat                                                                                 <  2:31911/251‑1 (MQ=255)                                                                                                                                                                              aCGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCGTTGCT‑GCCGCCAGCGTT‑CTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCatat                                                                               <  2:132426/249‑1 (MQ=255)                                                                                                                                                                              aCGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCGTTGCT‑GCCGCCAGCGTT‑CTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATatg                                                                              <  1:110177/250‑1 (MQ=255)                                                                                                                                                                                       aCAGCCTGAAATGGTCCGAGCATTATTCTCGATAGTCCGCCAGAAATAGCGGGCCCATTGGCGTTTTTGGAAGATCATG‑CCGGTCAGTGGAATCCTCCGCGACCAGATACA‑GC‑CCTT‑CCGCGTTTGCCGTTGCT‑GCCGCCAGCGTT‑CTGTCCTTTGCCGCCAGCAACTCTGCTGTCCAGCTCTCCGTCTGGTTGTTGGCACCGGCTTGCAACAGGCTGCCGACAATATTCCTGACCATATGATGTACaaa                                                                     <  1:72161/250‑1 (MQ=255)                                                                                                                                                                                         aGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCGTTGCT‑GCCGCCAGCGTT‑CTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGc                                                                  <  1:254397/251‑1 (MQ=255)                                                                                                                                                                                              gaaagcttCCGAGCATTATTCTCGGTAATCCGCCAGAAATCGCGGGCCCGTTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCGTTGCT‑GCCGCCAGCTTT‑CGGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTc                                                              <  2:3521/244‑1 (MQ=255)                                                                                                                                                                                              gAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCGTTGCT‑GCCGCCAGCGTT‑CTGTCCTTGGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTc                                                              <  2:390002/250‑1 (MQ=255)                                                                                                                                                                                                 aTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAACTAGCGGGCCCATTGGCGGTTTTTGGAGATAATT‑CCGGTCAGGGGAATCCACCGCGCCCCGATACA‑GC‑CCTT‑CCCCATTTTCCGTTGGT‑GCCCCCCGCGTT‑CTTTCCTTTTCCGCCCGCACCTCTCCTATTCAGCTCTCCCGCTGGGTGTGGGGACCGCCTTCCCTCCAGCTGCCGACAACATTTCTTACCATATTATGTACaaaaacctttgctt                                                          >  2:59514/1‑241 (MQ=255)                                                                                                                                                                                                        gAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCGTTGCT‑GCCGCCAGCGTT‑CTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTCGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATc                                                   <  2:22879/251‑1 (MQ=255)                                                                                                                                                                                                                ttCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTCGAGATCATA‑ACGGGCCGGGGAATACACCCCGACCAGATACA‑GC‑CCTT‑TCCCTTTTTCCGTTTCT‑GCCCCCCCCGGT‑CTTTCCTTTGCCGCCCGCAACTCTGCTCTCCCGCTCTCCCGCTGGGTGTGGGGGCCGACTTCCCTCATGCGGCCGACAACATTCTTCACCCCATGATGTTCaaaaacaattccttttatatcttcccccc                                           >  2:57699/1‑226 (MQ=255)                                                                                                                                                                                                                  ctcGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATA‑CCGGTCAGGGGAATCCACCGCGACCAGATACA‑AC‑CCCT‑CCGCTTTTTCCGTTTCT‑GCCGCCAGCGGT‑CTGTCCTTTTCCGCCCGCAACTCTGCTATTCAGCTCTCCGGCTGGGTGTGGGCACCGCCTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTTATATCTTCcaccac                                          >  1:136877/1‑250 (MQ=255)                                                                                                                                                                                                                    cGTTAGTCCGCCATAAATAGCGGGCCCATGGGCGGTTTTGGAAGATCATC‑CCGTTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTGGCCGTTGCG‑GCCGCCAGCGTT‑CTGTCCTTTGCCGCCAGCAACGCGGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGTCTTCCATCAGGCTGCCGACAATATTCCTGACCATAGGATGTACAAAAGCATTCGCGTTAATATCTACcaccac                                          <  2:391321/248‑1 (MQ=255)                                                                                                                                                                                                                              ccAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCGTTGCT‑GCCGCCAGCGTT‑CTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCAGATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACAGAAGGACCGTGa                             <  1:512589/251‑1 (MQ=255)                                                                                                                                                                                                                                      aGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCGTTGCT‑GCCGCCAGCGTT‑CTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGa                      <  1:596194/250‑1 (MQ=255)                                                                                                                                                                                                                                                 tGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTG‑CCGCTTTTGCCGTTGCT‑GCCGCCAGCGTT‑CTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCCCCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGc           <  2:162937/250‑1 (MQ=255)                                                                                                                                                                                                                                                    cGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCGTTGCT‑GCCGCCAGCGTT‑CTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATaa       <  1:456796/251‑1 (MQ=255)                                                                                                                                                                                                                                                          tGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTCGCCGTTGCG‑GCCGCCAGCGTT‑CTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGTCTGGTTGTGGGCACCGACTTCCATCAGGCGGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCTCCACATAAGGACCGTGACGCCTGACATTAATATGCATAACGTTg  <  2:366403/250‑1 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                               ACCAGACCGGTTTCACAGAACAAAATTAAGAACAAAATGGCATAAACCCAGACGCCGTACTCCGCGACCAGCTCTGCCAGATGCACGTCAATGTGCAGGATAAAATCGATGAGGAAATAAATCAGGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATA‑CCGGTCAGGGTAATCCACCGCGACCAGATACA‑GC‑CCTT‑CCGCTTTTGCCGTTGCT‑GCCGCCAGCGTT‑CTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATATGATGTACAAAAGCATTCGCTTTAATATCTACCACCACATAAGGACCGTGACGCGTGACATTAATATGCATAACGTTG  >  NC_000913/2434611‑2435109

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