Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,636,423	0	T	C	72.2%	53.6 / 43.1	72	intergenic (‑56/‑333)	nohA/ynfO	putative prophage DNA‑packaging protein NohA/DUF3950 domain‑containing protein YnfO
Reads supporting (aligned to +/- strand):  ref base T (7/13);  new base C (23/29);  total (30/42)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGT‑CATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGA  >  NC_000913/1636177‑1636664                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  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                                                                                                                                                                                                                     >  2:264369/1‑176 (MQ=17)                                                                                                       cAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCAAAAAGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCa                                                                                                                                                                                                                     <  1:264369/176‑1 (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)                                                                                                                                                                                                                                              aaaaaGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATg   <  2:7389/250‑1 (MQ=255)                                                                                                                                                                                                                                               aaaaGCTCGCTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGAGGGTATTTTAAATAAAAACATTAAGTTATGACGAAGAAGAACGGAAACGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATgc  <  2:72849/251‑2 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  CGGCGCAGCTTTTCGTTCTCAATTTCAGCATCCCTTTCGGCATACCATCTTATAACGGCGGCAGAGT‑CATAAAGCACCTCATTACCCTTGCCACCGCCTCGCAGAACGGGCATTCCCTGTTCCTGCCAGTTCTGAATGGTACGGATACTCGCACCGAAAATGTCAGCCAGCTGCTTTTTGTTGACTTCCATTGTTCATTCCACGGCCAAAAACAGAGAAAGGAAACGACAGAGGCCCAAAAGCTCGTTTTCAGCACCTGTCGTTTCCTTTCTTTTCAGGGGGTATTTTAAATAAAAACATTAAGTTACGACGAAGAAGAACGGAAATGCCTTAAACCGGAAAATTTTCATAAATAGCGAAAACCCGCGAGGTCGCCGCCCCGTAACCTGTCGGATCACCGGAAAGGACCCGTAAAGTGATAATGATTATCATCTACATATCACAACGTGCGTGGAGGCCATCAAACCACGTCAAATAATCAATTATGA  >  NC_000913/1636177‑1636664

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