Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	623,865	5	.	C	41.9%	6.2 / 91.5	43	intergenic (‑131/‑244)	fepB/entC	ferrienterobactin ABC transporter periplasmic binding protein/isochorismate synthase 1
Reads supporting (aligned to +/- strand):  ref base . (10/15);  new base C (5/13);  total (15/28)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

GGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGTACAGCAGACCATGGCAACACTTGCG  >  NC_000913/623569‑624159                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                                          ggggtccgatggtggAAACAGTACGCGGCGCCTGGCTTTCCAGCGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTTTTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTGGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑ccc                                                                                                                                                                                                                                                                                                                 <  1:551973/291‑1 (MQ=255) gggTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑ccc                                                                                                                                                                                                                                                                                                                 <  1:1048187/300‑1 (MQ=255) gggTGACGCTGGTGAAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑ccc                                                                                                                                                                                                                                                                                                                 <  1:980856/300‑1 (MQ=255)        gttggaggAATAAATCCGCTGCGGCTGGCTTTCCGTTGTAGTGGTGCCACGGCTGTCAGTAATCTGACGCGGCAAGTAAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGGGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑ccc                                                                                                                                                                                                                                                                                                                 <  1:1134767/291‑1 (MQ=255)                aaaCAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGAACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACACCAGCCCTCa                                                                                                                                                                                                                                                                                      >  2:71591/1‑300 (MQ=255)                  acaatgcgccggcgcgtgcgttCCAGCGTATGTGTGCCACGGCTGTGAGTAATCGGGCGGGGCCGGTCAGCGACCTGAACTGCGGCTTTGCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGGGGTAGCGCGGGGCGAGTCCCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGGCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGccc                                                                                                                                                                                                                                                                                         <  2:444121/275‑1 (MQ=255)                          cTGCGGCGGTCTTTCCAGGGTATGTGTGTCACGGCGGTCAGTAAACGGACGCGGCCAGTCAGGGGCGTGAACTGCGGCTATTCCGGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGGGCGGGGCGATTCTCACAAATCAGCTTACTGTTATTAAAAGGGTTAAGGGCGTAATGACAAATTCGAAAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAg                                                                                                                                                                                                                                                                             <  2:463090/299‑1 (MQ=38)                               gctggctTTCCAGGGTATGTGTGCCACGGCTTGCAGGAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCAccc                                                                                                                                                                                                                                                                                   <  2:700319/288‑1 (MQ=255)                                 tggctTTCCAGTGTATGTGTGCCACGGTTGTCAGTAAGCTGCCGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTGGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCAccccagccc                                                                                                                                                                                                                                                                             <  1:1201310/292‑7 (MQ=255)                                  ggctTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCATCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGCGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAC‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCAccccagccctcaccctg                                                                                                                                                                                                                                                                     >  2:971441/1‑285 (MQ=255)                                      ttCCAGTGGATGGGTGCCACGGTGGTAAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCGGAAAGAAAAAGTCCTGTTAATAGAAGGGCGTTGGGGTAGAGCGGGGCGAGTCTCCCAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATGCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCAccccagccctcaccctggaaggg                                                                                                                                                                                                                                                               <  2:1208769/301‑21 (MQ=255)                                         cAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCAccccagccctcaccctggaaggg                                                                                                                                                                                                                                                               >  2:271057/1‑278 (MQ=255)                                             gtATGGGGGCCACGCTTGCCAGTAATCTGACGCGGCCCGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCGGTTAATAGAAGGGCGTTGCGGTAGGGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAGGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGGGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCAccccagccc                                                                                                                                                                                                                                                                             <  1:228692/280‑7 (MQ=255)                                                   gtgCAACGGCTGTCAGTAAGCTGACGCGGCAAGTAAGCGGCCTGAACTGCGGCTATACCTGAAAGCAAAAGTACTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCAAGTCTCACAAATCAGCTTCCTGTTATTAATAAAGTTAAGGGCGTAATAAAAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCAccccagccctcaccctggaagggagaggg                                                                                                                                                                                                                                                         <  1:868002/294‑27 (MQ=38)                                                                  gTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑ccc                                                                                                                                                                                                                                                                                                                 >  2:802180/1‑235 (MQ=255)                                                                  gTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑ccc                                                                                                                                                                                                                                                                                                                 <  1:802201/235‑1 (MQ=255)                                                                        tGACGCGGCCAGGCAGCGGCCTGAACTGCGGCTGTTCCGACAGGCAAAAGTCCTGTTAATGGGAGGGCGTTGCGGTAGAGCGGGGCAGGTTTAACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAAGTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGGGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACAt                                                                                                                                                                                                                                            <  2:86724/298‑1 (MQ=255)                                                                          aCGGGGCCGGTCAGCGGCCTGAACTGCGGCTATTCCGGAAAGCAAAAGTCCTGTTGATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGGGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACAt                                                                                                                                                                                                                                            <  1:724324/296‑1 (MQ=255)                                                                              ggCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATcac                                                                                                                                                                                                                                              <  1:85180/278‑1 (MQ=255)                                                                                   gcaaGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAGGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCaat                                                                                                                                                                                                                     <  2:206983/295‑1 (MQ=255)                                                                                        cGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGGGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATcc                                                                                                                                                                                                                           <  1:557358/299‑1 (MQ=255)                                                                                                      gCTATTCCTGAAAGCAAAAGGCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCAt                                                                                                                                                                                                   <  1:71596/297‑1 (MQ=255)                                                                                                               gAAAGCAAAAGTGCTGTTAATGAGCGGGCGGTGGGGTAGAGCGGGGCGAGTCGAAAAAATCAGGTGCCTGTGATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAACCCGTCCCCTCGACCCTGTGGGGAGAGGGTTAGAGTGAGGGGGACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGTGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTc                                                                                                                                                                                                     <  1:953467/298‑1 (MQ=255)                                                                                                                        aGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACcac                                                                                                                                                                                          >  1:813167/1‑300 (MQ=255)                                                                                                                           ccTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTCAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAAcc                                                                                                                                                                            >  2:276549/1‑299 (MQ=255)                                                                                                                                     gAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAg                                                                                                                                                                 >  2:335296/1‑300 (MQ=255)                                                                                                                                             ttGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAACCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGt                                                                                                                                                                      >  2:1180968/1‑299 (MQ=255)                                                                                                                                                                 tctcACAAATCACCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTAGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAagag                                                                                                                                               <  1:720196/290‑1 (MQ=255)                                                                                                                                                                 tctaACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACAt                                                                                                                                        <  1:335311/297‑1 (MQ=255)                                                                                                                                                                  cggcaacgtcagatgtgttaaagagacagggTTAAGGGCGTAAGAACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGGGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCt                                                                                                                                                     <  2:308401/266‑1 (MQ=255)                                                                                                                                                                      caAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAagag                                                                                                                                               >  1:66977/1‑297 (MQ=255)                                                                                                                                                                        aaTCAGCTTCCTGTTATTAATAAGGTTAAGTGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAagag                                                                                                                                               <  2:66973/295‑1 (MQ=255)                                                                                                                                                                          tCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAACCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACAt                                                                                                                                        <  1:1181006/300‑1 (MQ=255)                                                                                                                                                                               ttCGTGTTATGAATAAGGTTAAGAGGGGAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGGGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTg                                                                                                                                    <  2:813146/299‑1 (MQ=255)                                                                                                                                                                               ttCCTGTTATGAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTGGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTg                                                                                                                                    <  1:6884/299‑1 (MQ=255)                                                                                                                                                                                    gTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTAGGCGAAGCCATCAATAATGCTTCTCATTTTCATTTTAACCACAACCAGATGCAACCCCGAGTTGAAGATTTCGTTACCCATACGGTTGAGATAGTgcctgattgc                                                                                                                           >  2:1259222/1‑283 (MQ=38)                                                                                                                                                                                              aGGTTAAGGGCGTAATGACAAATGCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTGAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAgg                                                                                                                    <  1:591367/300‑1 (MQ=255)                                                                                                                                                                                               ggTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCt                                                                                                                                                     >  1:308415/1‑266 (MQ=38)                                                                                                                                                                                               ggTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGCGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCCTCACCCCAGCCATCACCCCAGCCCTCACCCTGGACGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATCATGCTTCTCATTTTCATTGTAACCACAACCAGATTCAAACCCGAATTGCCGATTGCGTTACCTCTAGCGTTGACATAGTGCACGTTTGCTTTTAGttt                                                                                                                  >  2:217477/1‑301 (MQ=255)                                                                                                                                                                                                       gCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAAt                                                                                               <  1:372156/300‑1 (MQ=255)                                                                                                                                                                                                           aaggaaaaaTGCGAAAAAGGGCAAAATCCGGCCCCTCGCCCCGTTGGGGAGAGGGGTAGGGGGAGGGGAACGGGCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGGGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATg                                                                                              <  2:1081394/291‑1 (MQ=255)                                                                                                                                                                                                                  aatgtgACAAGGCGACCAACCCGGGCCCGCGCCCATTTGGGAGGAGGGTTAGGTTGGGGGGAACAGCCAGCACTGGTGCGAACAGTAACCCTCACCGCAGCCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATaaa                                                                                                <  2:1103141/295‑1 (MQ=255)                                                                                                                                                                                                                     tCGACAGCGCGCACCATCAGTCCCCGCGACGCTGGGGGGAGAGGTTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCGCAGCCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATAAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGa                                                                                             <  2:998644/300‑1 (MQ=255)                                                                                                                                                                                                                          aaaGCGCACAATCCGTCCCCGCGCCCCTTGGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCa                                                                                       <  1:971466/301‑1 (MQ=255)                                                                                                                                                                                                                                    acccGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAACCCTCACCCCAGCCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCCGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCCTTATTAAAg                                                                              >  1:1420288/3‑300 (MQ=255)                                                                                                                                                                                                                                                        ggggAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGa                                                >  2:887632/1‑298 (MQ=255)                                                                                                                                                                                                                                                               gggTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑CCACCCCAGCCCTCACCCTGGCAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCCGATGCAACCCCGAGTTGCAGATTGCGTTACCTCATGAGTTGACATAGTGCGCGTTTGCATTTAGGTTAGCGACCGAAAATATAAATGATCATCATTATTAAAGCCTTTATCATCTTGTGGAGGATGATATGGTTACGTc                                          >  1:733658/1‑297 (MQ=255)                                                                                                                                                                                                                                                                                     cagcACTGGTGCGAACACTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACAGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGTAcagc                      <  1:1259264/295‑1 (MQ=255)                                                                                                                                                                                                                                                                                                    cATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGTCATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGTACAGCAGACCATGGCAACACTTgcg  >  1:685253/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                                          GGGTGACGCTGGTGGAAACAATACGCTGCGGCTGGCTTTCCAGTGTATGTGTGCCACGGCTGTCAGTAATCTGACGCGGCCAGTCAGCGGCCTGAACTGCGGCTATTCCTGAAAGCAAAAGTCCTGTTAATAGAAGGGCGTTGCGGTAGAGCGGGGCGAGTCTCACAAATCAGCTTCCTGTTATTAATAAGGTTAAGGGCGTAATGACAAATTCGACAAAGCGCACAATCCGTCCCCTCGCCCCTTTGGGGAGAGGGTTAGGGTGAGGGGAACAGCCAGCACTGGTGCGAACATTAA‑‑‑‑‑‑‑‑‑‑‑‑CCC‑‑‑‑‑‑‑‑‑‑‑‑TCACCCCAGCCCTCACCCTGGAAGGGAGAGGGGGCAGAACGGCGCAGGACATCACATTGCGCTTATGCGAATCCATCAATAATGCTTCTCATTTTCATTGTAACCACAACCAGATGCAACCCCGAGTTGCAGATTGCGTTACCTCAAGAGTTGACATAGTGCGCGTTTGCTTTTAGGTTAGCGACCGAAAATATAAATGATAATCATTATTAAAGCCTTTATCATTTTGTGGAGGATGATATGGATACGTCACTGGCTGAGGAAGTACAGCAGACCATGGCAACACTTGCG  >  NC_000913/623569‑624159

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑