Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,636,413	0	C	A	73.0%	61.7 / 46.5	74	intergenic (‑46/‑343)	nohA/ynfO	putative prophage DNA‑packaging protein NohA/DUF3950 domain‑containing protein YnfO
Reads supporting (aligned to +/- strand):  ref base C (7/13);  new base A (24/30);  total (31/43)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGT‑CATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACG  >  NC_000913/1636165‑1636646                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                          ccttcAACGTCCCGGCTCAGCTTTTCGTTCCCAATTGCAGCCCCCCTTCCGTCTTACCGCCTAATAACGCCGGCAGAGG‑CATAAAGTACCTCATTACCCTTGCCACCGCCTCGAAGAACTGGCACTCCCTGTTCCTGCCAGTTCTTAATGGTACGGATGTTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGGTGACTCCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGACACGACCGAGGCCCa                                                                                                                                                                                                                                          <  2:268115/249‑1 (MQ=9)          ccccGGCGCCGGTTGTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATACCGGCGGCAGAGT‑CATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCgct                                                                                                                                                                                                                                <  2:200913/250‑3 (MQ=255)           cccGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGT‑CATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCgctt                                                                                                                                                                                                                               >  1:350283/1‑248 (MQ=25)            ccgccgctgCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATTTTATAACGGCGGCAGAGT‑CATAAAGCCCCTCATTACCCTGGCCACCCCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGATGCCAAAAAGCTCgcctt                                                                                                                                                                                                                              <  2:238969/243‑5 (MQ=11)             cGGCGCAGCTTTTCGTTTCCAATTCCAGCATCCCTTTCGGCATCCCATCTTATAACGGCGGCAGAGT‑CATAAAGCACCCCATTCCCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTc                                                                                                                                                                                                                             <  2:190034/251‑1 (MQ=11)                       tttCGTTCTCAATTTCACCAGCCCTTTCGGCATACCAACTTAAAACGCCGGCAGAGG‑CTAAAAGCTCCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGACAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTg                                                                                                                                                                                                                     <  2:182638/249‑1 (MQ=11)                          cGTTATCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGT‑AATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACAGAAAATGTCAGCCAGATGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGCGAAAGGAAAAGACCGAGGCCAAAAAGCTCGCTTTCAGCACATGTCGtt                                                                                                                                                                                                                >  1:313200/1‑251 (MQ=11)                           gTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGT‑CATAAAGCACCTCATTACCATTGCCACCGCCTCGCAGAAAGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCCACGAAAATGTCAGCCAACTGCTTTTTGTTGACTTCCATTGTTCATTTCACGGCCAAAAACAAGGAAAGGAAACAACAGAGGCCAAAAAGCTCGCTTTCAGCCCCTGTCGttt                                                                                                                                                                                                               >  2:370421/1‑251 (MQ=11)                               ccAATGTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGT‑CATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCctt                                                                                                                                                                                                           <  1:303454/250‑1 (MQ=19)                                      cagactCCCTTTCGGCATCCCATTTTAAAACGTCGGGAGAGT‑CATACATCCCCCCATTACCCTTCCCACCGCCCCGCAGAACGGGCATTCCCTTTGCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAAGGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTGTCTTTTc                                                                                                                                                                                                    <  2:327513/246‑1 (MQ=11)                                         cATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGG‑CATAAAGCCCCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCGCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAgg                                                                                                                                                                                                 <  2:196937/251‑1 (MQ=255)                                               tttcggattaccagcgtataACGGCGGCAGAGT‑CATAAAGCCCCCCATTCCCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGACAAAAAGCTCGCTGTCAGCACCTGTCGTTGCCTTTCTTTTCAGAGGGTa                                                                                                                                                                                            <  2:170211/234‑1 (MQ=11)                                                 tCGGCATACCATCTTATAACGGCGGCAGAGT‑CATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTAtt                                                                                                                                                                                          <  1:246018/250‑1 (MQ=25)                                                  cGGCATTCCATCTTATAACGGCGGCAGAGT‑CATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTAttt                                                                                                                                                                                         <  2:375811/250‑1 (MQ=255)                                                  cGGCATCCCAACTTATAACGGCGGCAGAGT‑CATAAAGCACCTCATTCCCCTTGCCACCGCCTCGCAGAACGGGCTTTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTAtttt                                                                                                                                                                                        <  2:60265/251‑1 (MQ=255)                                                  cGGCAAACCATTGTATTACGGCGGCAGAGT‑CATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTAtttt                                                                                                                                                                                        <  1:257672/251‑1 (MQ=255)                                                   ggCATACCATCTTATAACGGCGGCAGAGT‑CATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGGAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTAttttt                                                                                                                                                                                       >  1:352694/1‑250 (MQ=255)                                                   ggCATACCATCTTATAACGGCGGCAGAGT‑CATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTAtttt                                                                                                                                                                                        >  1:270363/1‑250 (MQ=25)                                                        aCCATCTTATAACGGCGGCAGAGT‑CATAAAGCCCCCCATTACCTTGGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTtaaata                                                                                                                                                                                   <  2:104186/250‑1 (MQ=255)                                                          cATCTTATAACGGCGGCAGAGT‑CATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATAATCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCCTTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACCGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTATTTTCAGAGGGTATTTtaaataaa                                                                                                                                                                                 >  1:352602/1‑250 (MQ=11)                                                               tataACGGCGGCAGAGT‑CATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGGACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCTGGGGGTATTTTACATAAAAACAtt                                                                                                                                                                           >  2:99732/1‑251 (MQ=255)                                                                 aaaCGGCGGCAGAGTACA‑AAAGCCCCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCTTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACATGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTaa                                                                                                                                                                         <  2:103853/250‑1 (MQ=11)                                                                  aaCGGCGGCAGAGT‑CATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGCATGGGACGGATACTCGCACCGAAAATGTCAGCCCGCTGCTTTTTGTTTACTTCCATTGTTCATTCCACGGCCCAAAACTGAGAAAGGAAACGACAGAGGCCAAAAAACTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTTAATAAAAACATTaa                                                                                                                                                                         >  2:318523/1‑250 (MQ=11)                                                                        ggcAGAGT‑CATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTatg                                                                                                                                                                   >  1:150177/1‑248 (MQ=16)                                                                           agagC‑CATAAAGACCCTCATTTCCCTTGCCCCCGCCTCGCAGAACGCGCATCTCCTGTCTCTGCCATTTCTTAATGGTACGTATGCTCGCACCTAAAATGTCAGCCAGCCGCTTTTTGTTGACTTCCATTGTTCATCCCACGGCCAAAAACTGAGAATGGAAACGACAGAGTCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGATGGCATTTTAAATAAAAACATTAAGTTatgtcga                                                                                                                                                               <  1:261362/251‑7 (MQ=2)                                                                             agT‑CATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAACACATTAAGTTAcgacga                                                                                                                                                               >  2:100089/1‑249 (MQ=39)                                                                                        cACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGa                                                                                                                                                    >  1:40267/1‑250 (MQ=16)                                                                                          ccTCATTCCCCCTTCCACCGCCCCGCAGAACGGGCATTCCCTTTTCCTGCCAGTTCTTAAGGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTGTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGaaac                                                                                                                                                 <  2:328604/251‑2 (MQ=11)                                                                                               gaccccTTGCCACCGCCACGCAGAACGGGCATTCCCTGGTCCTGCCAGTGCTGAATGGTACGGATACTCGCACCGAAAATGGCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATACCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTATTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCtt                                                                                                                                            <  2:192907/248‑1 (MQ=14)                                                                                                 aCCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTaa                                                                                                                                          <  1:99732/251‑1 (MQ=255)                                                                                                                   cAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCa                                                                                                                                                                                                   <  1:264369/176‑1 (MQ=17)                                                                                                                   cAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCa                                                                                                                                                                                                   >  2:264369/1‑176 (MQ=17)                                                                                                                     gAACTGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGGTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGACTTTTAAAAAAAAACATTAAGTTATGACGAAGAAGAACGGAAAGGCCTTAAACCGGAAATTTTTCATaaa                                                                                                                       >  2:337777/1‑250 (MQ=11)                                                                                                                          ggCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTACTTTCTTTTCAGGGGGTATTTTAAATAAAAACCTTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAACTTTTCATAAATAGCg                                                                                                                  >  2:297386/1‑250 (MQ=255)                                                                                                                             aTTCCCTTTTCCTGCCAGTTCTGAATGGTACGTATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCCTTCCACGGCCAAAAACAGAGAAAGGAACCGACAGAGGCCCAAAAGCTCGTTTTCAGCCCCTGTCGTTTCCTTTCTTTACAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGCAGAAGAACGAAAATGCATTAAACCGGAAAATTTTAATAAATAGCGaaa                                                                                                               >  1:53571/1‑250 (MQ=14)                                                                                                                                 ccttttCCGGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCcg                                                                                                          <  2:224807/247‑1 (MQ=11)                                                                                                                                   tGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCTATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACATAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAACACCcgcg                                                                                                        >  2:83967/1‑251 (MQ=11)                                                                                                                                    tttCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCcgcg                                                                                                        <  2:159702/249‑1 (MQ=255)                                                                                                                                       ccTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGACAGGCAACGCCAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGACGAAGACCGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGt                                                                                                    >  2:280706/1‑251 (MQ=11)                                                                                                                                        cTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTc                                                                                                   >  1:11023/1‑251 (MQ=255)                                                                                                                                         tGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTcc                                                                                                  >  1:355672/1‑250 (MQ=255)                                                                                                                                         tGCCAGTTCTGAATGGAACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAACTAAAAACATTAAGTTATGCCGAAGAAGACAGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGACAACCAGCGAGGTcg                                                                                                  >  2:229652/1‑251 (MQ=11)                                                                                                                                          gCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTcg                                                                                                  <  2:285984/250‑1 (MQ=255)                                                                                                                                              gTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTGCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGGCGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCgccgcc                                                                                             <  1:261790/251‑1 (MQ=11)                                                                                                                                                          aCGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTc                                                                                 <  1:129624/251‑1 (MQ=255)                                                                                                                                                             gATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCgg                                                                               >  1:119996/1‑250 (MQ=255)                                                                                                                                                                            aaTGTCAGCCAGATGCTCTTTGTTGACTTCCATAGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAATAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCTCCGGAAAGGAcc                                                               <  2:111600/251‑1 (MQ=255)                                                                                                                                                                                    ccagcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGt                                                       >  1:347761/1‑251 (MQ=255)                                                                                                                                                                                      agcTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTg                                                      >  1:323822/1‑250 (MQ=255)                                                                                                                                                                                          gctTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAAt                                                 <  1:221027/251‑1 (MQ=255)                                                                                                                                                                                             ttttGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAACCGCCTTAAACCGGAAAATTTTCAt                                                                                                                          >  2:270902/1‑175 (MQ=11)                                                                                                                                                                                             ttttGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCAt                                                                                                                          <  1:270902/175‑1 (MQ=11)                                                                                                                                                                                              tttGTTGACTTCCATTGTTCATTCCTCGGCCAAAAACATAGAAATGAAACGACAGAGGCCAAAAAGCTCGCGTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTAGTTTAAATAAAAACAGTAAGTTAAGACGAAGAATAACGGAAACGCCTTAAACCGGAAAATTTTCATAAACAGCGAAAACCCGCGAGGACGCAGCCCCGTATCCTGTCGGATCAACTGAAATGACCCGTAAAGTGATAATGAt                                              <  2:239284/250‑1 (MQ=14)                                                                                                                                                                                                tgttgCCTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTat                                           >  1:41438/1‑251 (MQ=255)                                                                                                                                                                                                 gttgACGTTCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTatc                                          <  1:301682/251‑1 (MQ=14)                                                                                                                                                                                                 gttaaCTTCCATTGTTCATTCCACGGCCAAAAACTGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTatc                                          <  1:318523/247‑1 (MQ=14)                                                                                                                                                                                                       ttCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGTCAGAGGCCAAAACGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCTCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTAc                                    >  2:37321/1‑251 (MQ=14)                                                                                                                                                                                                          cATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACat                                  >  1:91352/1‑250 (MQ=255)                                                                                                                                                                                                              gTTCATTCCCCGGCCAAAAACAGAGATAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGATGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATca                              <  2:257391/250‑1 (MQ=14)                                                                                                                                                                                                              gTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGGTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCTGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATca                              <  2:302844/250‑1 (MQ=14)                                                                                                                                                                                                                 cATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGTGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGAtt                                             <  1:116981/232‑1 (MQ=14)                                                                                                                                                                                                                 cATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGTGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGGGATAATGAtt                                             >  2:116981/1‑232 (MQ=11)                                                                                                                                                                                                                    tCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAAcgt                        <  2:356000/250‑1 (MQ=255)                                                                                                                                                                                                                      caatgCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCAGAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAAcgtgcg                     <  2:373280/247‑1 (MQ=255)                                                                                                                                                                                                                         ggCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTTGTTTTCAGCACCTGACGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTGACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCTGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTgg                  <  2:41438/251‑1 (MQ=255)                                                                                                                                                                                                                          gCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGa                 <  1:333683/251‑1 (MQ=255)                                                                                                                                                                                                                           ccAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAg                <  1:253113/251‑1 (MQ=255)                                                                                                                                                                                                                           ccAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAg                <  2:11023/251‑1 (MQ=255)                                                                                                                                                                                                                              aaaaCAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGcc             >  2:102440/1‑251 (MQ=255)                                                                                                                                                                                                                                   agagAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCCCCGCCCCGTAACCTTTCGGATCCCCGGAAAGGACCCGTAAAGGGATAATGATTATCATCAACATATCACAACGTGCGTGGAGGCCATCa         >  1:332109/1‑250 (MQ=255)                                                                                                                                                                                                                                    gagaAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCaaa       <  2:336673/251‑1 (MQ=255)                                                                                                                                                                                                                                       aaaGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAAcc     >  1:340214/1‑250 (MQ=255)                                                                                                                                                                                                                                        aaGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCAc   <  2:347761/251‑1 (MQ=255)                                                                                                                                                                                                                                         agtaaaCGACAGAGGCCAAAAAGCCCGCTTTCAGCACCTGGAGTTTCCTGTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCTCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACg  <  2:271687/248‑1 (MQ=14)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                          TCTTCAACTTCCCGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGT‑CATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACG  >  NC_000913/1636165‑1636646

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