Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,161,142	0	C	A	24.0%	113.0 / 11.8	50	Q26K (CAA→AAA)	fabR	DNA‑binding transcriptional repressor FabR
Reads supporting (aligned to +/- strand):  ref base C (35/3);  new base A (0/12);  total (35/15)
Fisher's exact test for biased strand distribution p-value = 3.75e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGC  >  NC_000913/4160892‑4161390                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                          cGGCCGTTGCTTTGTCTGTTTCTTTCTTTTGCCCTGACGTTTTATTGCATTTTTATCTACGATTCACTCCGTTTAAGCGGCAGGGTTCCGCTTTTCGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATTTTTTATTGCGTTACCGTTCATTCACAATACTGGATCAATCAAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACAACGTGATTGGCGTAA‑GA‑GCGcaac                                                                                                                                                                                                                                                           <  2:115982/251‑1 (MQ=255)  ggCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGcaac                                                                                                                                                                                                                                                           >  1:332699/1‑250 (MQ=255)    ccTTTGCTGTGTCTGTTTCTTCCCTTTGCCCTGACGTTTTATTTGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTTCGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTGCTTTTGGCAA‑GG‑A‑TTCAGTCATCGTGATGGGCGTAA‑GA‑GCGCAACaaa                                                                                                                                                                                                                                                        <  2:459863/251‑1 (MQ=255)     ctttgctttgTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACaaa                                                                                                                                                                                                                                                        >  1:124071/1‑250 (MQ=255)                      cTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAAACC‑GCCGtt                                                                                                                                                                                                                                     >  2:477357/1‑251 (MQ=255)                        tACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCg                                                                                                                                                                                                                                   >  1:34472/1‑251 (MQ=255)                               ccccTGACGTTTTATTGGATTTTTATCGACGTTACCCCCCGTTAAAGCGGCGGTTTTCCGTTGTAGGAAAAAGACCCGGCCAAAGAATTGAAGTAAATATTTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAACCCAGTATTTTCATTCTCTGTTATAGTGCC‑AGCATTACTTTTGGAAA‑GA‑A‑TTCAGCCATCGTGATGGG‑GTAAAGA‑GCGCAAAAAAAAAAAAAAACCC‑GCCGTTCGCtggtgg                                                                                                                                                                                                                            <  2:124071/250‑1 (MQ=255)                                                gATTTTTCTCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCa                                                                                                                                                                                                            >  1:494056/1‑250 (MQ=255)                                                                 ctctCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTaa                                                                                                                                                                                                       >  2:509029/1‑238 (MQ=255)                                                                 ctccccGTTTAGGCGGCAGGTTTCCGCTGTACGTAAAAGACCCGGCCAAAGAATTGCAGAAAATATGTTTTATTGCGTTACCGTTCATTCACAATCCTGGAGCAATCCAGTATGTTCATTCTCTGGTATGGTGCC‑AGCAGTACTTTGGGCAA‑GG‑A‑TTCAGACATCGTGAGGGGCGTAA‑GA‑GCGCAACAAAAAAAAAAACCCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTaa                                                                                                                                                                                                       <  1:509029/234‑1 (MQ=255)                                                                        tttAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGcc                                                                                                                                                                                   >  1:20031/1‑251 (MQ=255)                                                                           aaGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGc                                                                                                                                                                                >  2:321523/1‑251 (MQ=255)                                                                               ggCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTg                                                                                                                                                                             >  1:589328/1‑250 (MQ=255)                                                                                  agggttCCGATGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGAGTTACCGTTCATTCACAATACTGGAGCAATACAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGtt                                                                                                                                                                         >  2:94442/5‑251 (MQ=255)                                                                                   ngTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGttt                                                                                                                                                                        >  2:264453/2‑251 (MQ=255)                                                                                       tCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGt                                                                                                                                                                    >  2:210668/1‑251 (MQ=255)                                                                                          ggaaaacggaaaaaacccGGCCAAAAAATTGGAGAAAATATGTTTTATGGCTTTACCGTTAATTCACAAGACTGGAACAATCCAGTATTTTCATTCTCTGTTAAAG‑GCCGATCAGTACTTTTGGAAA‑GA‑A‑TTCAGACAGCTTGATGGGCAAAA‑GG‑GCACAAAAAAAAAAAAAACCCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGaa                                                                                                                                                                 <  2:318432/235‑1 (MQ=255)                                                                                             gTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTg                                                                                                                                                              >  1:348994/1‑251 (MQ=255)                                                                                                                agatcttgacagaaaaTATTTTTTATTGCGTTACCGTCAATTAAAAATGTGGGAGCAACCCAATATGTTCATTCTCGGGTATAGTGCC‑AGCAGTATTTGTGACAA‑GG‑A‑TTCAGACATCGTTAGGGGGGAAA‑GA‑GCGAAAAAAAAAAAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATGAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTg                                                                                                                                           <  2:513111/238‑1 (MQ=255)                                                                                                                    aaTTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTcc                                                                                                                                       >  2:534952/1‑251 (MQ=255)                                                                                                                     aTTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCAtt                                                                                                                                            >  1:342320/1‑245 (MQ=255)                                                                                                                      ttGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGtt                                                                                                                                                                                                                                     >  1:10740/1‑155 (MQ=255)                                                                                                                      ttGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTg                                                                                                                                           >  1:513111/1‑245 (MQ=255)                                                                                                                          agtaaattgtttttATTGCTTTCCCGTTCATCAACAATCCTGGAGAAATCCATAATGTTCTTCCTCTGGTATGGGGCA‑GGCATTA‑TTTTGGGAA‑GGGA‑TCAAGCCATCGTGAGGGGCGAAA‑GG‑GCGAAAAAAAAAAAAAAACCCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACctc                                                                                                                                 <  2:32654/241‑1 (MQ=255)                                                                                                                                           gggttcctgtccattcaAAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGAAGTATTTTTGGCAA‑GG‑A‑TTCAGACATCGTGTGGGGCTTAA‑GA‑GCCCAAAAAAAAAAAAAACCCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTc                                                                                                                 <  2:75013/239‑1 (MQ=255)                                                                                                                                             gTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCcgc                                                                                                              >  1:34053/1‑251 (MQ=255)                                                                                                                                                  cGTTCATTAACAATACTGTAGAAATCAAGTATGTTCTTTCTCTGGTATAGGGCC‑AGCAGTACTTTTGGAAA‑GG‑A‑TTAAGACATCGGGAGGGGCGTAA‑GA‑GCGCAAAAAAAAAAAAAACCCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGt                                                                                                         <  2:586000/251‑1 (MQ=255)                                                                                                                                                          cacaATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGGGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAAc                                                                                                 >  2:356950/1‑251 (MQ=255)                                                                                                                                                           acaATCCTGGAGCAACCCGGTATGTTCTTTCTCGGGAGTAGAGCC‑ACCAGTACTTTTGAAAA‑GG‑A‑TACAGACATCGTAGGGGGCG‑AA‑GAAGC‑CAAAAAAAAAAAAAAACCCCGCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACt                                                                                                <  2:589328/251‑1 (MQ=255)                                                                                                                                                                 cTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCt                                                                                          >  2:85140/1‑251 (MQ=255)                                                                                                                                                                   ggAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTg                                                                                         >  1:126658/1‑250 (MQ=255)                                                                                                                                                                       aaaTCCAGTTTTTTCATTCTCTGGTATAGTCCC‑AGCAGTACTTTGGGCAA‑GG‑A‑TTCAGAAATCGTGATGGGCGTAA‑GA‑GCGAAAAAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCAt                                                                                    <  1:210668/250‑1 (MQ=255)                                                                                                                                                                              gTATGTTCATTCTCTGGTATATGGCC‑AG‑AATTTTTTTGGGAA‑GG‑AATCCAGACATGGGGGTGGGAGTAA‑GA‑CCGAAAAAAAAAAAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTtgat                                                                             <  2:348994/251‑1 (MQ=255)                                                                                                                                                                                  tttCATTCTCGGGTATGGTGCC‑AGCAGTA‑TTTTGGAAAAGG‑A‑TAAAGACATCGTGGTGGGCGAAA‑GA‑GCACAAAAAAAAAAAAAACCCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATgag                                                                          <  2:126658/249‑1 (MQ=255)                                                                                                                                                                                          tctGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGttt                                                                  >  1:223636/1‑250 (MQ=255)                                                                                                                                                                                             ggTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATg                                                              >  1:51559/1‑251 (MQ=255)                                                                                                                                                                                                ataGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGCGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGCCGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCt                                                            >  2:435579/1‑250 (MQ=255)                                                                                                                                                                                                       c‑agcagTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAAc                                                    >  1:451610/1‑251 (MQ=255)                                                                                                                                                                                                       c‑agcagTACTTTTGGCAA‑GG‑A‑TTAAGACATCGTGAGGGGCGAAA‑GA‑GCGAAACAAAAAAAAAAACCCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAAc                                                    <  1:534952/251‑1 (MQ=255)                                                                                                                                                                                                             gagttacttttggagagatatCAGTCATGGTGAGGGGCATAA‑GA‑GCGCAAAAAAAAAAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTc                                                  <  2:332699/231‑1 (MQ=255)                                                                                                                                                                                                                     ggCAA‑GG‑A‑TTCAAACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGgc                                       >  1:426550/1‑251 (MQ=255)                                                                                                                                                                                                                       cAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGgcgc                                     >  1:515494/1‑251 (MQ=255)                                                                                                                                                                                                                        aa‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGcgcg                                    >  1:279537/1‑251 (MQ=255)                                                                                                                                                                                                                           gg‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGAc                                                                                                     >  2:584534/1‑184 (MQ=255)                                                                                                                                                                                                                            g‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTc                                  >  1:560504/1‑250 (MQ=255)                                                                                                                                                                                                                                 tCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCttttt                                                                                                                            >  2:507071/1‑157 (MQ=255)                                                                                                                                                                                                                                 tCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCCTGCGCCAGGCGCGTCAGCg                              >  2:253934/1‑251 (MQ=255)                                                                                                                                                                                                                                 ccAGACATCGTGATGGGCGTAA‑GA‑GCGCAAAAAAAAGAAAAACCCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCttttt                                                                                                                            <  1:507071/156‑1 (MQ=255)                                                                                                                                                                                                                                   aGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGt                             >  2:340471/1‑250 (MQ=255)                                                                                                                                                                                                                                               gggCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCg                >  2:316562/1‑251 (MQ=255)                                                                                                                                                                                                                                                ggCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCgg               >  1:586817/1‑251 (MQ=255)                                                                                                                                                                                                                                                      a‑ga‑gCGCAACAAACAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAgtgg         >  2:567790/1‑250 (MQ=255)                                                                                                                                                                                                                                                           gcacaacaaAAAAAAAAACCCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGAt      <  1:193416/248‑1 (MQ=255)                                                                                                                                                                                                                                                             acaacaaAAAAAAAAACCCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATcc    <  1:382319/250‑1 (MQ=255)                                                                                                                                                                                                                                                               aacaaAAAGAAAAAACCC‑GCCGTTCGCTGGTGGCAGCCGCATTTAGCCAATAAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCCTGGCTGATGAGAGCGGTTTAATGCTACGCCAACTCACGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGCTCCGc  >  2:493951/1‑251 (MQ=255)                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                          CGGCCTTTGCTTTGTCTGTTTCTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGC  >  NC_000913/4160892‑4161390

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