Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,161,148	0	G	A	26.9%	107.5 / 14.0	52	E28K (GAA→AAA)	fabR	DNA‑binding transcriptional repressor FabR
Reads supporting (aligned to +/- strand):  ref base G (36/2);  new base A (0/14);  total (36/16)
Fisher's exact test for biased strand distribution p-value = 6.78e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTC  >  NC_000913/4160913‑4161395                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                         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)                                                                                                                                                                                                                                               aaaGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGTCGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCAcccc  >  2:141433/1‑249 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                         CTTACTTTTGCCCTGACGTTTTATTGGATTTTTATCGACGATACTCTCCGTTTAAGCGGCAGGTTTCCGCTGTACGTAAAAGAACCGGCCAAAGAATTGCAGTAAATATGTTTTATTGCGTTACCGTTCATTCACAATACTGGAGCAATCCAGTATGTTCATTCTCTGGTATAGTGCC‑AGCAGTACTTTTGGCAA‑GG‑A‑TTCAGACATCGTGATGGGCGTAA‑GA‑GCGCAACAAAAAGAAAAAACCC‑GCCGTTCGCTGGTGGAAGCCGCATTTAGCCAATTAAGTGCTGAACGCAGCTTCGCCAGCCTGAGTTTGCGTGAAGTGGCGCGTGAAGCGGGCATTGCTCCCACCTCTTTTTATCGGCATTTCCGCGACGTAGACGAACTGGGTCTGACCATGGTTGATGAGAGCGGTTTAATGCTACGCCAACTCATGCGCCAGGCGCGTCAGCGTATCGCCAAAGGCGGGAGTGTGATCCGCACCTC  >  NC_000913/4160913‑4161395

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