Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,423,338	0	G	C	47.7%	52.7 / 78.1	45	intergenic (+26/‑62)	ydaW/rzoR	putative uncharacterized protein YdaW/putative prophage outer membrane lipoprotein RzoR
Reads supporting (aligned to +/- strand):  ref base G (11/12);  new base C (9/12);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 7.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.39e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGAC  >  NC_000913/1423099‑1423579                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  ggcaatCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGagatgcagatgc                                                                                                                                                                                                                                         <  1:63265/250‑1 (MQ=255)        aGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGt                                                                                                                                                                                                                                 <  2:18160/251‑1 (MQ=255)            aaGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATg                                                                                                                                                                                                                             >  1:374028/1‑251 (MQ=255)                      gATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTtgctgc                                                                                                                                                                                                                      <  1:27034/248‑1 (MQ=255)                            cgtaactaaGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTACTCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGCAGCGTCAGCGTGATGTTGCTGCGCTCGATGc                                                                                                                                                                                                             >  2:40786/9‑251 (MQ=255)                                       gCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTAGCTCGGCTCTGGTTATCGTCAGCGTCGTCTGCCTGGCATGGGGTGTTAATCATTACCGTGATAACGCCATGGCCTACAAAGAACGGCGTGTCAAAAAAGCCAGAGAACTGAAGCTGTCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTTCGCTCTATGCAAAATACACGa                                                                                                                                                                                                  <  1:465108/251‑1 (MQ=255)                                            aaaTTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAgg                                                                                                                                                                                               >  2:140958/1‑249 (MQ=255)                                               tttCGTTAATGCATACCGGGCGATAATTTACCAGCAGAGCCTACAGCAGCATCTTAATCTACGCTCTGGTTATCTGCCTCGTCGTCTGCATGTAATGGGGTGTTACTCATTACCGTGATACCGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCACATACTGAAATGCAGAAGCGTCAGCGTGCTGTTGCTGCGCTctcgtccaaatacacgatgtagtt                                                                                                                                                                                           >  2:271026/1‑227 (MQ=18)                                                ttCGTTACTGCATACCGGGCGATAAATTACCAGCAGGGCCGAAAGCAGCATCGGTAGCTCCGCTCTGGTTATCTGCATCGTCGTCTGCTTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCTCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGAGGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTa                                                                                                                                                                                          <  2:252325/250‑1 (MQ=255)                                                   gTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACATCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGAt                                                                                                                                                                                                               <  1:14222/226‑1 (MQ=255)                                                   attACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGCAGCTGGCGAACGCAGCAATTACTGAGATGCAGATGCGTCAGCGTGGTGTTGCTGCGCTCGAt                                                                                                                                                                                                               >  2:14222/2‑226 (MQ=255)                                                         gCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCg                                                                                                                                                                                 <  2:374028/250‑1 (MQ=255)                                                                             ccAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTg                                                                                                                                                            >  1:338609/1‑251 (MQ=255)                                                                                             caccgcgattATCTCCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGcc                                                                                                                                             <  1:315943/241‑1 (MQ=11)                                                                                                    ggTATCTCCGCTCTGGTTATCTGCATCGTCGGCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGgtcg                                                                                                                                     <  2:7281/251‑1 (MQ=255)                                                                                                       acctcCGCTCTTGTTATCTGCATCGTCGTCTGCCTGTCATGCGCTTTTAATCATTACCGTGATAACGCCATTACATACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATTCTCTGCGTGATGATGTTGCCGCTGgtcgtc                                                                                                                                   <  2:82844/248‑1 (MQ=11)                                                                                                        tctcCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGc                                                                                                                                                           >  2:250338/1‑225 (MQ=11)                                                                                                        tctcCGCTCTGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGc                                                                                                                                                           <  1:250338/225‑1 (MQ=11)                                                                                                               tctGGTTATCTGCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGgt                                                                                                                                       >  1:445828/1‑238 (MQ=11)                                                                                                                          gCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGTACATCAAAg                                                                                                                 <  2:351283/249‑1 (MQ=255)                                                                                                                          gCATCATCGTCTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAAGGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCAGTTGCACATCaa                                                                                                                   <  2:42037/247‑1 (MQ=11)                                                                                                                              cgtcgtCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCt                                                                                                           >  1:233888/1‑251 (MQ=255)                                                                                                                              cggcgtcTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCATAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTTATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTc                                                                                                            <  2:59507/247‑1 (MQ=255)                                                                                                                                gcggcTGCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGGACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCGCGATGCAAAATACACGAAGGAGTTAGCTGAGGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTGGCCGCGGGTCGTCGTCGGTTGCACATCAAAGCAGTCTgt                                                                                                         <  2:416473/247‑1 (MQ=11)                                                                                                                                      gCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTg                                                                                                                           <  1:319026/227‑1 (MQ=11)                                                                                                                                      gCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTg                                                                                                                           >  2:319026/1‑227 (MQ=11)                                                                                                                                      gCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAACGCAGTCTgtccgtc                                                                                                    >  2:167163/1‑250 (MQ=11)                                                                                                                                      gCCTGTCATGGGCTGTTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTgtcagt                                                                                                     <  1:34395/249‑1 (MQ=11)                                                                                                                                               gggctgtTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGaa                                                                                          <  1:185381/247‑1 (MQ=11)                                                                                                                                                  ctgtTAATCATTACCGTGATAACGCCATTACCTACAAAGCCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGcc                                                                                       >  1:89558/2‑251 (MQ=255)                                                                                                                                                                              cGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGcaca                                                                                                                       <  2:310630/191‑1 (MQ=39)                                                                                                                                                                              cGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGcaca                                                                                                                       >  1:310630/1‑191 (MQ=39)                                                                                                                                                                                ccTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAACATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAACTGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTgt                                                                                                         >  1:308218/1‑203 (MQ=255)                                                                                                                                                                                ccTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTgt                                                                                                         <  2:308218/203‑1 (MQ=38)                                                                                                                                                                                     aaagcccAGCGAGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTcccc                                                    >  1:59356/7‑251 (MQ=2)                                                                                                                                                                                       agcccAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGACCGCGGCAATTACTGACATGCAGATGCGTCAGCGTAATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTcccccc                                                  <  1:132739/247‑1 (MQ=9)                                                                                                                                                                                       agcccAGCGCGACAAAAATGCCAGAGAACGGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTcccccc                                                  <  1:269797/247‑1 (MQ=9)                                                                                                                                                                                        gAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCg                                                 >  1:201214/1‑251 (MQ=35)                                                                                                                                                                                          ccaacacgACAAAAAAGCCAGAGACCTAAAGATGGCGCACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGACATTGCTGCGCTCGCTGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAACCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGCCGTGGATAATGAAGCCTCCCCCCGAc                                               >  2:281720/7‑251 (MQ=11)                                                                                                                                                                                           cAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACt                                              <  1:375598/251‑1 (MQ=34)                                                                                                                                                                                                    aaaaaTGCCAGAGAACGGAAGCTGGCGAACGCGGCAATTACGGACATGCAGAGGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCGGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCCCCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGacac                                     <  1:10535/251‑1 (MQ=11)                                                                                                                                                                                                      aaaTGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTc                                                                            >  2:117206/1‑210 (MQ=14)                                                                                                                                                                                                      aaaTGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACACCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTc                                                                            <  1:117206/210‑1 (MQ=9)                                                                                                                                                                                                                       aGCTGGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGACGCTGGTCGTCGTCGGTTGAACATAAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCACCCGACTGGAAGACACCGCTGAACGGGATTATTTc                  >  2:172523/1‑251 (MQ=11)                                                                                                                                                                                                                              gAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGATGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAACGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATAATTTCACCCtac           >  1:164048/1‑249 (MQ=16)                                                                                                                                                                                                                                       aaTTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGAc  >  1:401895/1‑251 (MQ=38)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                  GGCAATCAGCAAAGGTCGTCTGATTCGCGTAAATCAGGGCGGTAAATTTCGTTACTGCATACCGGGCGATAATTTACCAGCAGAGCCGAAAGCAGCATCGGTATCTCCGCTCTGGTTATCTGCATCGTCGTCTGCCTGTCATGGGGTGTTAATCATTACCGTGATAACGCCATCGCCTACAAAGAACAGCGCGACAAAAATGCCAGAGAACTGAAGCTGGCGAACGCGGCAATTACTGAGATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTGATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGGGAACGAC  >  NC_000913/1423099‑1423579

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