Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,065,598	0	G	A	28.6%	68.1 / 10.2	42	L57L (CTG→TTG)	yihO	putative sulfoquinovose transporter
Reads supporting (aligned to +/- strand):  ref base G (24/6);  new base A (1/11);  total (25/17)
Fisher's exact test for biased strand distribution p-value = 2.85e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATCGCACCATACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑GAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTACGGGCTTCTATCTC  >  NC_000913/4065348‑4065838                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 atcGCACCATACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑g                                                                                                                                                                                                                                                  <  1:579875/251‑1 (MQ=255)    gCACCATACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑Gaa                                                                                                                                                                                                                                                >  2:174168/1‑250 (MQ=255)       ccATACGAGCAGGTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑GAAAGAt                                                                                                                                                                                                                                            <  1:86252/251‑1 (MQ=255)                cagatcatcaGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGAGGTGGCGATAAACTGGAGCGTTTCAATTAACGCCGCAGGAACCGCTGCATCTCAAATAACGGTGCGAAATTTGCCTTTAAGCCCAATA‑TTTTT‑CCGCGAGTCGAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGGAAAAAA‑C‑TTCGCCACCA‑GACAGATAATcccacc                                                                                                                                                                                                                                   >  2:714475/5‑248 (MQ=255)                                   aaaGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑GAAAGATAATCCCGCCATAGTAAGCAGGCATCCCt                                                                                                                                                                                                                >  1:725333/1‑251 (MQ=255)                                            acataaaaaaCAGCGGGGGGGCAATCGTCGTTTAAACCGGAAAAAAAAAGGGGGTGGCAATAAACGGGAGCGTTGAAATTAACGCCGCCGGAACCGCTGAATATAAAATAAAGGGACAAAATTTGCCTTTAGGCCAAATA‑TTTTC‑CCGCGGGTCGATTAAAAAC‑CCGGGGAGCATATCGGAAAC‑CGCGGTAAAAAC‑C‑TCCGCCCCCA‑GAAAGATAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCg                                                                                                                                                                                                       <  1:408912/244‑1 (MQ=255)                                                catgaaCAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑GAAAGATAAt                                                                                                                                                                                                                                         >  1:550815/5‑213 (MQ=255)                                                catAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCCCCA‑GAAAGATAAt                                                                                                                                                                                                                                         <  2:550815/213‑1 (MQ=255)                                                 atAAACAGAGCGGTAGCAATCGTCGTTTTAACCGGCAAACAACAGGTGATGGCGATAAACTGTAGCGTTGCCATTAACGCCGCCGGAACCACAGCATATACAATAAAGGGACGAAATTTGCCTTTAGCCCCAATA‑TTTTT‑CCTCGAGTCGAGTAAAACT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑ATCGCAGCCA‑AAAAGATAATCCCGCCATAATAAGCAGGCATCACTAATTCATCggaat                                                                                                                                                                                                   >  2:729586/1‑247 (MQ=255)                                                         cgcgGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑GAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTg                                                                                                                                                                                           >  2:101709/1‑250 (MQ=255)                                                                  aatcttgttgttaaacgcaaaaaaaaaGGTGGTGGGAATAAACTGGTGCGTTGAATTTACCGCCGCCGGACCCGCTGAATATAAAAAAAGGGGTCAAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCAAGGCGAGTAAAAAT‑CCGGTGAG‑ATTAGGGTAAAACTCGTTAAAAAC‑C‑TTCCCCCCAA‑AAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACa                                                                                                                                                                                 <  1:816636/228‑1 (MQ=255)                                                                                  ggCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑GAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGAATAAAATTTGAGGAAATAAAGAGTACCGATGCag                                                                                                                                                                   >  2:574685/1‑249 (MQ=255)                                                                                   gCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCAGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGCGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑GAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGGGTACCGATGCAAAGc                                                                                                                                                                >  2:615884/1‑251 (MQ=255)                                                                                     aaaaaaagggtggtggaaaaaaaCGGGAGCGTTACAATAAACGCCGCCGGAACCGCTGAATAAAAAATAAAGGGACGAAATTTGCTTTTAGGCCAAATA‑TTTTC‑CCGGGAGTCGAGAAAAAAT‑CGGGTGAGCATATGGGTAAC‑CGCGGAAAAAAA‑C‑TCCGCCACAA‑AAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCAt                                                                                                                                                              <  1:200261/231‑1 (MQ=255)                                                                                          aaaGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑GAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAAt                                                                                                                                                         >  1:11633/1‑251 (MQ=255)                                                                                              gtggtgGCGATAAGCTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑GAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATAGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCg                                                                                                                                                      >  2:107199/1‑250 (MQ=255)                                                                                                  tggCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑GAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCg                                                                                                                                                 >  1:20587/1‑251 (MQ=255)                                                                                                  gggCGAAAACCTGGAGCGTTGCAATTAACGCCGCCGAAACCTCTGCATATAAAAAAAAGGAACAAAATTTCCTTTTAGGCCAAATA‑TTTTC‑CCGCGAGTCGAAAAAAAAC‑CCGGGGAGCATTCCGGTAAA‑CGCGGAAAAAAC‑C‑TTCCCCCCCA‑GAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATGCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCg                                                                                                                                                 <  1:113985/250‑1 (MQ=255)                                                                                                    gCGAAAAACGGGAGCGTTGCAATTAACGCCCCCGGACCCGCTACAAAAAAAATAAAGGGACAAAATTCGCCTTTAGGCCCAATA‑TTTTC‑CCGCGAGGCGATAAAAAAC‑CCGGTGAGAATATCGGTAAA‑CGCGGAAAAAAC‑C‑TCCCCCACCA‑AAAAGAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGAc                                                                                                                                               <  1:761274/251‑1 (MQ=255)                                                                                                             tGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑GAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCccat                                                                                                                                      >  1:753662/1‑251 (MQ=255)                                                                                                                aGCGTTACAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAACAAA‑C‑TTCGCCACCA‑GAAAGATAATCCCGCCATAGTCAGCAGGCATCCCTAATTCATCGGTATAAACTTTGAGGAGATACCGGGTACCGATGCAGAGCCTTAAATTCGAACCGACGTCGCccatacc                                                                                                                                   >  2:323599/1‑251 (MQ=255)                                                                                                                         attaccccctgcgGAACCGCTGCAAATAAAATAAAGGTACTAAATTTGCCTTTAGTCACAATT‑TTTTT‑CCGAGAGGCGATTAAAAAC‑CCGGTGAACTTATCTGTAAA‑CGCGGTAAAAAA‑C‑CCCGCCCCAA‑AAAAGATACTCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGAt                                                                                                                          <  1:280584/240‑1 (MQ=255)                                                                                                                           tAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑GAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGCTGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATCGGCGAttt                                                                                                                        >  1:77116/1‑251 (MQ=255)                                                                                                                               gccgccGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGAGTACAA‑T‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑GAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTcc                                                                                                                   >  1:394023/1‑251 (MQ=255)                                                                                                                                 cgccgGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑GAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTccc                                                                                                                  >  1:616214/1‑251 (MQ=255)                                                                                                                                       aaCCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑GAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGc                                                                                                                >  1:247806/1‑247 (MQ=255)                                                                                                                                         ccGCTGAATAAAAAAGAAAGGAACAAAATTGGCCTTAAGGCC‑AATAATTTTC‑CCGGGAGTAAAAAAAAAAC‑CCGGTGGGCATA‑CCGTAAATCGCGGAAAAAAC‑C‑TCCCCCACAA‑AAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTc                                                                                                          <  2:417893/251‑1 (MQ=255)                                                                                                                                          cgctgtaaaaaaaataaaGGGACTAAATTTCCCTTTAGGCCCACTA‑TTTTC‑CCGCGATGCGTGTAAAAAC‑CCGGTGAGCATTTCGGAAAA‑CGCG‑TAAAAAA‑CCTTCCCCCCCA‑AAAAGATAATCCCTCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCa                                                                                                         <  1:649351/241‑1 (MQ=255)                                                                                                                                           gCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGAGAGTCCAGTAAAAAT‑CCGGTGCGAATATCGGCAAA‑CGCGGTAACACA‑C‑TTCGCCACCA‑GAAAGCTAGTACCGCCATCGCAAGCAGTAATCCCTAATCAATCGGTATAAAATTTAAGGAGATACAGAGTACCGATGCAGCGCATTAAATTAGACCCGCCGTCGCCCATCCCATAGGCGATTTTTTCCCGAAGaatc                                                                                                          >  1:529406/1‑246 (MQ=255)                                                                                                                                               cATATAAAAAAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑GAAAGATAATCCCGCCATAGTAAGAAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATAAAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGtttt                                                                                                    >  2:533290/1‑251 (MQ=255)                                                                                                                                                 tataAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TGTTT‑CCGCGAGTCAAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTAGCCACCA‑GAACGATAATCCCGCCATAGTAAGAAGGCATACCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCAAACCGACGTCGCCCATACCATAAGCGATTTTTTccagacgtttcaggttta                                                                                                   >  1:669381/1‑234 (MQ=255)                                                                                                                                                      aataaaGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑GAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGCTTTTTTCCCGCAGATTCAGTTTTAATGTc                                                                                              >  1:290089/1‑250 (MQ=255)                                                                                                                                                           aGGGACAAAATTTGC‑TTTTGGCCCAATT‑TTTTTTCCGCGGGTCGA‑AAAAAATGCCGGTAGGCATA‑CGGGAAA‑CCGGGAAAAAAA‑C‑CTCCCCCCCAAGAAAAAAAACCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGa                                                                                        <  1:209469/251‑1 (MQ=255)                                                                                                                                                                       tGCCTTTAGCCCAAATA‑TTTTC‑CCGCGAGACGAATAAAAAC‑CCGGTGAGAATATGGGAAAA‑CGCG‑TAAAAAA‑CCTTCGCCCCAA‑AAAAAATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGa                                                                             <  1:101709/250‑1 (MQ=255)                                                                                                                                                                                         atgttgccGCGAGTAGAGAAAAAAC‑CGGGGGAGCATATCGAAAAA‑CGCGGAAAAAAA‑C‑TGCGCCACCA‑AAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTtata                                                            <  1:728024/245‑1 (MQ=255)                                                                                                                                                                                             ttcccgggattggattaaacaccgggTGAGCATATCGGTAAA‑CG‑GGTAAAAAAAC‑TTCGCCCCCA‑AAAAGAAAATACCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAAc                                                          <  1:61628/226‑1 (MQ=255)                                                                                                                                                                                                        gagtAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑GAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGcg                                             >  1:273481/1‑251 (MQ=255)                                                                                                                                                                                                             aaaaT‑CCGGTGAGCATACCGGAAAC‑CGCGATAAAAAA‑C‑TCCGCCCCCA‑AAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGtt                                        <  1:252225/251‑1 (MQ=255)                                                                                                                                                                                                                 t‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑GAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAAt                                                                                                 >  1:811087/1‑190 (MQ=255)                                                                                                                                                                                                                     ggTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑GAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTc                                                                                                          >  2:474826/1‑178 (MQ=255)                                                                                                                                                                                                                     ggTGAGCATATCGGTAAA‑CGCGAAAAAAAC‑C‑TCCCCCACAA‑AAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTc                                                                                                          <  1:474826/178‑1 (MQ=255)                                                                                                                                                                                                                                       acggggTAAAAAA‑C‑TTCCCCACCA‑GAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGGCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTa                <  1:311231/247‑1 (MQ=255)                                                                                                                                                                                                                                                       ttCGCCACCA‑GAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTACGGGCTTCTAtctc  >  1:301523/1‑251 (MQ=255)                                                                                                                                                                                                                                                                  gAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTCATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAgg                                                >  2:230191/1‑195 (MQ=255)                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 ATCGCACCATACGAGCAGTTCATCAGACTGTATGAAAGTCCAAACATCATAAACAGCGCGGTGGCAATCGTCGTTTTAACCGGCAAACAAAAGGTGGTGGCGATAAACTGGAGCGTTGCAATTAACGCCGCCGGAACCGCTGCATATAAAATAAAGGGACGAAATTTGCCTTTAGGCCCAATA‑TTTTT‑CCGCGAGTCGAGTAAAAAT‑CCGGTGAGCATATCGGTAAA‑CGCGGTAAAAAA‑C‑TTCGCCACCA‑GAAAGATAATCCCGCCATAGTAAGCAGGCATCCCTAATTCATCGGTATAAAATTTGAGGAGATACAGAGTACCGATGCAGAGCATTAAATTCGAACCGACGTCGCCCATACCATAGGCGATTTTTTCCCGCAGATTCAGTTTTAATGTCAGTGGATTATGGTCAGACATAATCACTCCTTATAACGCCCCGCAGGGCGTCGTTTATTATTATTACGCCGTGCGTTTACGGGCTTCTATCTC  >  NC_000913/4065348‑4065838

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