Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,684,736	0	A	C	32.4%	129.7 / 10.3	69	T160P (ACC→CCC)	tus	DNA replication terminus site‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (40/6);  new base C (0/22);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 9.40e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TTATTTATTCAACAACAGTCCGAAAATCGCAGCAGCAAGGCCGCTGTCCGTCTGCCTGGCGTGTTGTGTTACCAGGTCGATAACCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAACGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAG  >  NC_000913/1684487‑1684969                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          ttatttatTCAACAACAGTCCGAAAATCGCAGCAGCAAGGCCGCTGTCCGTCTGCCTGGCGTGTTGTGTTACCAGGTCGATAACCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCAc                                                                                                                                                                                                                                          >  2:135488/1‑251 (MQ=255)   atttattCAACAACAGTCCGAAAATCGCAGCAGCAAGGCCGCTGTCCGTCTGCCTGGCGTGTTGTGTTACCAGGTCGATAACCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCg                                                                                                                                                                                                                                        >  2:741675/1‑251 (MQ=255)       nttCAACAACAGTCCGAAAATCGCAGCAGCAAGGCCGCTGTCCGTCTGCCTGGCGTGTTGTGTTACCAGGTCGATAACCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCt                                                                                                                                                                                                                                    >  2:344137/2‑251 (MQ=255)       attCAACAACAGTCCGAAAATCGCAGCAGCAAGGCCGCTGCCCGTCTGCCTGGCGTGTTGTGTTACCAGGTCGATAACCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTc                                                                                                                                                                                                                                     <  2:545646/250‑1 (MQ=255)         tCAACAACAGTCCGAAAATCGCAGCAGCAAGGCCGCTGTCCGTCTGCCTGGCGTGTTGTGTTACCAGGTCGATAACCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGc                                                                                                                                                                                                                                  >  2:690249/1‑251 (MQ=255)             caacaGTCCGAAAATCGCAGCAGCAAGGCCGCTGTCCGTCTGCCTGGCGTGTTGTGTTACCAGGTCGATAACCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACg                                                                                                                                                                                                                               >  2:421902/1‑250 (MQ=255)              aacaGGCCGAAAATCGCAGCAGCAAGGCCGCTGTCCGTCTGCCTGGCGTGTTGTGTTACCAGGTCGATAACCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGAc                                                                                                                                                                                                                             <  2:425451/251‑1 (MQ=255)               acaGTCCGAAAATCGCAGCAGCACGGCCGCTGTCAGACTGCCTGGCGTGTTGAGTTACCAGGCCGATAACCCTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGAcc                                                                                                                                                                                                                            >  2:692069/1‑251 (MQ=255)                 aGTCCGAAAATCGCAGCAGCAAGGCCGCTGTCCGTCTGCCTGGCGTGTTGTGTTACCAGGTCGATAACCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGAcccc                                                                                                                                                                                                                          <  1:415236/251‑1 (MQ=255)                             gcagcagcaAGGCCGCTGTCCGTCTGCCTGGCGTGTTGTGTTACCAGGTCGATAACCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGc                                                                                                                                                                                                              >  2:642749/1‑251 (MQ=255)                              cagcagcaAGGCCGCTGTCCGTCTGCCTGGCGTGTTGTGTTACCAGGTCGATGACCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGACTGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACttt                                                                                                                                                                                                                  >  1:59800/1‑246 (MQ=255)                                  agcaAGGCCGCTGTCCGTCTGCCTGGCGTGTTGTGTTACCAGGTCGATAACCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTttgg                                                                                                                                                                                                         >  1:389379/1‑251 (MQ=255)                                     aaGGCCGCTGTCCGTCTGCCTGGCGTGTTGTGTTACCAGGTCGATAACCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTttggttg                                                                                                                                                                                                      >  2:290368/1‑251 (MQ=255)                                     aaGGCCGCTGTCCGTCTGCCTGGCGTGTTGTGTTACCAGGTCGATAACCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTttggttg                                                                                                                                                                                                      >  2:699375/1‑251 (MQ=255)                                         ccGCTGTCCGTCTGCCTGGCGTGTTGTGTTACCAGGTCGATAACCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCt                                                                                                                                                                                                  >  1:60353/1‑251 (MQ=255)                                           gCTGTCCGTCTGCCTGGCGTGTTGTGTTACCAGGTCGATACCCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCtc                                                                                                                                                                                                 >  2:313172/1‑249 (MQ=255)                                           gCTGTCCGTCTGCCTGGCGTGTTGTGTTACCAGGTCGATAACCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCta                                                                                                                                                                                                 >  2:639245/1‑250 (MQ=255)                                            cgttccgcctgcgtgggggttttggttACCAGGGCAATACCCTTCCGCAAGCAGCGTTGGTCGGT‑CATATACAC‑CACATAAAAAAACCCAAGCCCCCGTACGAGCATATCGTCACGTTT‑AAATCA‑GACCCCCCCCCCGCGGCATTTTTTGAAGGGGTGCATCGTCATTGGCCGGGGCTGATCCCCCTTAATG‑CTT‑CCGGCCCGCTCCCTGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCtaat                                                                                                                                                                                               <  2:203713/227‑1 (MQ=255)                                              gTCCGTCTGCCTGGCGTGTTGTGTTACCAGGTCGATAACCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCCTTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCtaata                                                                                                                                                                                              >  2:568162/1‑250 (MQ=255)                                                      gCCTGGGGTGTTGTGTTCCCAGGTCAATAACCTTTCGCAAGCAGCGTTGGTAAGT‑ATTTTTCAG‑CACATAAATAACCCCAAGCCCACTTTCGAGCATATCGTCAGGGTT‑GAATCA‑GACCCCCCCCCCGCGGCACTTTTTGAAGGGGTGCATCGTCATTGGCGGGGGCTGATCCCCCTTAATG‑CTT‑CCCGCACGCTCCCCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATc                                                                                                                                                                                      <  1:576204/250‑1 (MQ=255)                                                           gggtgttttgttCCCAGGTCTATAACCTTTCGCAAGCAGCGTTGG‑CAGTTAATATT‑AGACACATCAATAACCTCAAGACCGCGTTCGAGCATATCGTCCCGGTT‑AAACAA‑GAACCCCCCCCCGCGGCACTTTTTAAGTGGGTGCATCGTCTTTTGCCGGGGCTGATCCCCCTTAATG‑CTT‑CCCGCACGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAg                                                                                                                                                                                <  1:500724/249‑1 (MQ=255)                                                           gcgggttgtgttGCCATGTCGATTCCCTTTCGCAAGCAGCGTGGGCCCGT‑CATATTCAG‑TACATCAATAAACTCAAGACGACGTTCGTGCATATCGTCAAGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCCCCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCCCGCTCACCGTTCCGCACGACCCCGCCACATTACGCTTAGCATGCGCTAATAAACAGATCATTa                                                                                                                                                                                  >  2:681366/5‑249 (MQ=255)                                                           gCGTGTTGTGTTACCAGGTCGATACCCTTCCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAAAAAACACAAGACCACGTTCGAGCATATCGTCACGGTT‑AAATCA‑GAACCCCCCCCCGCGGCACTTTTTGAATGGGTGCATCGTCTTTTGCCGGGGCTGATCCCCCTTAATG‑CTT‑CCCGCCCGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAg                                                                                                                                                                                <  1:335583/251‑1 (MQ=255)                                                            cGTGTTGTGTTACCAGGTCGATAACCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTa                                                                                                                                                                                  >  2:490264/1‑248 (MQ=255)                                                               gttgttttACCAGCTCGATACCTTTTCGCAAGCAGCGTGGGTCATT‑AATTTCCAG‑CCAATAAATAAACTAAAGCCCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GACCCCCCCACCGCGGCCCTTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCCCCTTTAATG‑CTT‑CCCGCGCGCCCCCCTTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAtt                                                                                                                                                                            <  2:39658/251‑1 (MQ=255)                                                                ttgtgttACCAGGTCGATAACCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAAttt                                                                                                                                                                           >  1:602809/1‑251 (MQ=255)                                                                     ttcccAGGTCAATAACCTTTCGCAAGCAGCGTTGGTCAGT‑TATATTCAG‑CACATAAATAAACTCAAGACCCCGTTCGAGCATATCTTCACGTTA‑GAATCA‑GAACCCCCCCCCGCGGCCCGTTTTGAATGGGTGCATCGTCTTTTGCCGGGGCTGATCCCCTTTAATG‑CTT‑ACCGCACGCTCACCTTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACAt                                                                                                                                                                       <  1:421902/247‑1 (MQ=255)                                                                      tACCAGGTCAATACCCTTCCGCAAGCAGCGTTGGTCAGT‑CATATCCAC‑CACATAAAAAAACCAAAGCCCCCGTTCGAGTATATCGTCACGGTT‑GAATCA‑AACCCCCCCCCCGCGGCCCTTTTTGAATGGGTGCATCGTCATTGGCCGGGGCTGACCCCCCTTAATG‑CTT‑CCCGCCCCCTCCCCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATc                                                                                                                                                                      <  1:690249/250‑1 (MQ=255)                                                                               gATAACCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGt                                                                                                                                                            >  1:461178/1‑251 (MQ=255)                                                                                    ccTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTgg                                                                                                                                                       >  2:549287/1‑251 (MQ=255)                                                                                          tcAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGc                                                                                                                                                 >  2:55411/2‑251 (MQ=255)                                                                                            aaGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTg                                                                                                                                               >  2:298179/1‑251 (MQ=255)                                                                                              gcagcGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTgg                                                                                                                                              >  1:585372/1‑250 (MQ=255)                                                                                                agcGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGaa                                                                                                                                            >  1:367670/1‑250 (MQ=255)                                                                                                          agt‑caTATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGa                                                                                                                                  >  1:636187/1‑250 (MQ=255)                                                                                                                 atTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGCTGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATcacc                                                                                                                           >  2:280680/1‑251 (MQ=255)                                                                                                                  ttCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATcacc                                                                                                                           >  1:735524/1‑250 (MQ=255)                                                                                                                   tCAG‑CACATCAAAAAACGCAAGCCCCCGTTCGATCATATCGTCACGGTT‑AAACCA‑GACCCCCCCCCCGCGGCACTTTTTGAATGGGTGATTCTTCTTTTGCGGGGGCTGACCCCCCTTAATG‑CTT‑CCCGCACGCCCCCCTTTCTGCCCGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGTCTGAAATcacca                                                                                                                          <  1:549287/250‑1 (MQ=255)                                                                                                                    tcaccacaacaaaaaccccAAGCCCACGTTCGAGCATATCGTCGCGGTT‑GAATCA‑GACCCCCCCCCCGCGGCACGTTTTGAATGGGTGCATCGTCATTGGCCGGGGCTGATCCCCCTTAATG‑CTT‑CCCGCACGCTCCCCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATcacca                                                                                                                          <  1:344137/232‑1 (MQ=255)                                                                                                                         acaTCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAgt                                                                                                                    >  2:54879/1‑251 (MQ=255)                                                                                                                            aaaaTAAACTCAACACCACGTTCGAGCATATCGTCACGGTT‑AAATCA‑CACCCCCCCCCCGCGGCCCGTTTGGAATGGGTCCATCGTCATTTGCCGGGGCTAATCCCCTTTAATG‑CTT‑CCCGCCCGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTc                                                                                                                 <  1:699375/249‑1 (MQ=255)                                                                                                                                  aaCTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGAGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTACTAAACATATCATTAAGAATTTAAATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCg                                                                                                           >  2:605937/1‑251 (MQ=255)                                                                                                                                    cacAAGACCTGGTTCGCGCATATCGTCGCGTTT‑GATTCA‑TAACTCCCCCCGGGGGCACGTTGTGAGGGGGTGCTTCGTCATGGGGGGGGGCTGACCCCCTTTATTG‑TGT‑CCCGCCCGCCCCCCGTCCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTg                                                                                                         <  2:502588/249‑1 (MQ=255)                                                                                                                                       aagcccACGTTCGAGCATATCGTCAGGTTT‑GAATCA‑GACCCCCCCCCCGCGGCAGTTTTTGAATGGGTGCTTCTTCTTTGGCCGGGGCTGATCCCCCTTAATG‑CTT‑CCCGCCCGCCCCCCGTTCTCCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAc                                                                                                      <  1:290368/247‑1 (MQ=255)                                                                                                                                        aacccACGTTCGAGCATTTCTCCCCGGTT‑GATTCA‑AACCTCCCCCCCGCGCCCCGTTTTGATTGGGTGCATCGCCTTTTGCCGGGGCTCATCCCCCTTA‑TGTCTT‑CCCGCCCGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcg                                                                                                     <  1:727005/248‑1 (MQ=255)                                                                                                                                        aGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcg                                                                                                     >  2:76103/1‑251 (MQ=255)                                                                                                                                         gACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcgc                                                                                                    >  2:71561/1‑251 (MQ=255)                                                                                                                                           ccccGTTCGGGCATATCGCCGCG‑TT‑TAATCAAGACCTCCCCCCCGCGGCCCGTTTTGAATGGGTTCATCGGCTTTTCCCGGGGCTGACCACCTTTAATC‑CTT‑CCCGCACGCTCCCCTTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTTGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcgcgc                                                                                                  <  1:691699/248‑1 (MQ=255)                                                                                                                                           ccACGTTCGAGCATACCGTCCCGGTT‑GAATCA‑AACCCCCCCCCCGCGCCACGTTTGGAAGGGGTGCATCGTCTTTTGCCGGGGCTGATCACCCTTAATT‑CTT‑ACCCCCCGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcgcgc                                                                                                  <  1:741675/251‑1 (MQ=255)                                                                                                                                            catgTTCGAGCATATCGTCCCGGTT‑AAACCA‑GACCCCCCCCCCGCGGCACGTTTTGAAGGGGTGCTTCTTCATTTGCCGGGGCTAACCCCCCTTATTG‑CTA‑CCCGCCCGCTCCCCGTTCTGCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcgcgcg                                                                                                 <  1:596021/248‑1 (MQ=255)                                                                                                                                            cACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcgcgcg                                                                                                 >  2:314085/1‑251 (MQ=255)                                                                                                                                            cACGTCCGAGCATATCG‑CCCGGTTTGAAACA‑GACCTCCCCCCCGCGGCACTTTTGGAATGGGTGCATCGTCATTTGCGGGGGCTGTTC‑CCCTTTATT‑CTTTCCCGCCCGCTCACCTTTCTGCACTCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcgcgcg                                                                                                 <  1:59063/251‑1 (MQ=255)                                                                                                                                             aCGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTACTAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGAAcgcgcgg                                                                                                >  2:603887/1‑250 (MQ=255)                                                                                                                                             aCGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGAcgcgcg                                                                                                 >  2:583363/1‑250 (MQ=255)                                                                                                                                                     gCATATCGTCAGGGTT‑GAAACA‑GACCCCCCCCCCGCGCCACGTTTTGAATGGGTGCCTCGTCATTTGCCGGGGCTGATCCCCCTTAATG‑CTT‑CCCGCACGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGc                                                                                        <  1:408284/251‑1 (MQ=255)                                                                                                                                                             cgacgcgatagatccataccTCCCCGCCGCGGCGTATTTTGAAGGGGTGCATCTTCATTTGCGGGGGCTGACCCCCCTTTGTG‑CTT‑CCCGCACCCTCACCGTTCTCCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAg                                                                                    <  2:314761/230‑1 (MQ=255)                                                                                                                                                              cagggTT‑GAATCA‑GACCCCCCCCCCGCGGCACGTTTTGAAGGGGTGCATCGTCATTTGCCGGGGCTGACCCCCCTAAATG‑CTT‑ACCGCCCGCTCCCCTTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAAc                                                                               <  1:583363/248‑1 (MQ=255)                                                                                                                                                                ggcggtgtgctctcccgctctccccGCGGCGCGTTTTGAGGGGGTCCATCGTCTTTTGCCGGGGCTGACCCCCTTTAATG‑CTT‑CCCGCCCGCTCCCCGTTCTGCACGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAAc                                                                               <  1:373851/228‑1 (MQ=255)                                                                                                                                                                  gTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTgg                                                                            >  2:409689/1‑250 (MQ=255)                                                                                                                                                                                  ccccACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAgg                                                             >  2:521093/1‑251 (MQ=255)                                                                                                                                                                                     ccccGCGGCCTGTTTTGAAGGGGTGCCTCGTTTTTTGCGGGGCCTTCTCCCCCTTAATG‑CTC‑ACCCCCCGCCCCCCTTTCTACCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGata                                                          <  2:446469/249‑1 (MQ=255)                                                                                                                                                                                     cACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAAc                                                                               >  2:373851/1‑230 (MQ=255)                                                                                                                                                                                               gattttGAATGGGGGCACCGTCGTTTGCGGGGGCTGACCCCCCTTAAGT‑CTT‑GCCGCCCGCTCCCCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGctgccct                                                <  2:585372/249‑1 (MQ=255)                                                                                                                                                                                                         gggTGCATCGTCATTTGCCGGGGCGGATCCCCTTTAACG‑CTA‑CCCGCCCGCTCTCCGTTCGGCCCAACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAAc                                      <  1:313172/251‑1 (MQ=255)                                                                                                                                                                                                               aTCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGATTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAACGCGa                                  >  2:206065/1‑249 (MQ=255)                                                                                                                                                                                                                            gggcttgaccccccTTAATG‑CTT‑CCCGCCCCCTCCCCGTTCCCCCCGCCCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAACGCGAAGTTAAAAATCAAAc                   <  1:314085/240‑1 (MQ=255)                                                                                                                                                                                                                                      aCCCTTAATG‑CGT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCGGCCACAGAACGCGAAGTTAAAAATCAAATGTCCGTTGa          >  2:713993/1‑250 (MQ=255)                                                                                                                                                                                                                                                tgtttcccGCACCCTCCCCTTTCTGCCCGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAACGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAg  <  1:568162/244‑1 (MQ=255)                                                                                                                                                                                                                                                          aCGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAAc                                      >  2:529268/1‑204 (MQ=255)                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                          TTATTTATTCAACAACAGTCCGAAAATCGCAGCAGCAAGGCCGCTGTCCGTCTGCCTGGCGTGTTGTGTTACCAGGTCGATAACCTTTCGCAAGCAGCGTTGGTCAGT‑CATATTCAG‑CACATCAATAAACTCAAGACCACGTTCGAGCATATCGTCACGGTT‑GAATCA‑GAACTCCCCACCGCGGCACGTTTTGAATGGGTGCATCGTCATTTGCCGGGGCTGATCACCCTTAATG‑CTT‑ACCGCACGCTCACCGTTCTGCACGACCCCGCCACTTTACGCTTTGGTTGGGCTAATAAACATATCATTAAGAATTTACATCGTGATGAAGTCCTGGCACAGCTGGAAAAAAGCCTGAAATCACCACGCAGTGTCGCACCGTGGACGCGCGAGGAGTGGCAAAGAAAACTGGAGCGAGAGTATCAGGATATCGCTGCCCTGCCACAGAACGCGAAGTTAAAAATCAAACGTCCGGTGAAGGTGCAG  >  NC_000913/1684487‑1684969

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