Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,808,226	0	A	G	45.6%	16.1 / 71.5	57	intergenic (‑290/‑265)	PP_4234/dsbD‑II	bifunctional aldehyde oxidase/xanthine dehydrogenase/thiol/disulfide interchange protein 2
Reads supporting (aligned to +/- strand):  ref base A (17/14);  new base G (9/17);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 1.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCT  >  NC_002947/4808083‑4808375                                                                                                                                                |                                                                                                                                                       aTGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGgcgctt                                                                                                                                                  >  2:5222054/1‑148 (MQ=255)  tGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGgcgcttc                                                                                                                                                 <  1:938043/150‑4 (MQ=255)    ggCGGTATTCGCGCTTTGGCGGGGAGGGGTTTGGCGGGAGGTTGTCAGCGGCGGGGCGGTGTTCGCCGGGGGATTTTCAGCGCCAGGGAGGGGGTCGCCGGGTGAGTTTCAGCGCCGGGGAGATCGAGCGCCGCCCGCGCGgcgcttcgc                                                                                                                                               <  2:3095126/150‑6 (MQ=25)    ggCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGgcgcttcgc                                                                                                                                               <  1:1559537/150‑6 (MQ=255)    ggCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGgcgcttcgc                                                                                                                                               <  2:3094860/150‑6 (MQ=255)    ggCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGgcgcttcgc                                                                                                                                               >  1:2518593/1‑145 (MQ=255)     gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCttcgcg                                                                                                                                              >  2:4368533/1‑150 (MQ=255)     gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCttcgcg                                                                                                                                              >  1:4479964/1‑150 (MQ=255)     gCGGAAGTCGCGCTGTGGCGTGTCGGCGTTTGGCGGGAGTCGCTCAGCGGCGGGGAGGTGTTCGGCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCttcgcg                                                                                                                                              <  2:932380/150‑1 (MQ=21)      cGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCCGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgg                                                                                                                                             <  2:5334280/150‑1 (MQ=38)      cGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgg                                                                                                                                             >  1:84920/1‑150 (MQ=255)      cGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgg                                                                                                                                             <  2:3799384/150‑1 (MQ=255)      cGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgg                                                                                                                                             <  2:2509015/150‑1 (MQ=255)      cGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgg                                                                                                                                             >  1:5210387/1‑150 (MQ=255)      cGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgg                                                                                                                                             <  1:4695419/150‑1 (MQ=255)      cGGAAGTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgg                                                                                                                                             <  1:3281712/150‑1 (MQ=38)        gAATTCGCGCTTTGGCGGGTAGGTGTTGGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggt                                                                                                                                           <  1:656003/150‑2 (MQ=38)         aaTTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggta                                                                                                                                          <  1:578509/150‑3 (MQ=255)            tCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggtaaac                                                                                                                                       <  2:4196784/150‑6 (MQ=38)            tCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggtaaac                                                                                                                                       <  1:1277806/150‑6 (MQ=38)             cgcgCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggtaaacc                                                                                                                                      <  2:4479964/150‑7 (MQ=38)              gcgcTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggtaaaccc                                                                                                                                     <  2:84920/150‑8 (MQ=38)              gcgcTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggtaaaccc                                                                                                                                     >  2:5260215/1‑143 (MQ=38)              gcgcTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggtaaaccc                                                                                                                                     >  1:2675950/1‑143 (MQ=38)                                                                         gggATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCttcgcg                                                                                                                                              >  1:1056198/1‑82 (MQ=18)                                                                         gggATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCttcgcg                                                                                                                                              <  2:1056198/82‑1 (MQ=18)                                                                                                                                  gggtttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCca                >  2:580092/8‑150 (MQ=255)                                                                                                                                  gggtttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCca                >  2:580050/8‑150 (MQ=255)                                                                                                                                   ggtttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCACCGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcac               <  2:4291124/144‑2 (MQ=255)                                                                                                                                   ggtttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa               <  1:5223853/144‑1 (MQ=255)                                                                                                                                   ggtttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa               <  2:3702953/144‑1 (MQ=255)                                                                                                                                   ggtttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa               <  2:4885388/144‑1 (MQ=255)                                                                                                                                   ggtttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa               <  2:1626354/144‑1 (MQ=255)                                                                                                                                   ggtttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa               <  1:2632577/144‑1 (MQ=255)                                                                                                                                   ggtttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa               >  1:676622/7‑150 (MQ=255)                                                                                                                                   ggtttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa               <  1:5224209/144‑1 (MQ=255)                                                                                                                                    gtttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaac              >  1:1377781/6‑150 (MQ=255)                                                                                                                                     tttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaaca             <  2:3550687/146‑1 (MQ=255)                                                                                                                                       tacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaacaac           <  2:3549111/148‑1 (MQ=255)                                                                                                                                        acccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTACACTCAAAATTCCAACGCCGCACCCTCAAAATACCAACGCCCCCTCCCCTCGAAAATCccacaaaa          >  1:5206965/2‑148 (MQ=39)                                                                                                                                         cccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaacaacaa         <  1:1637445/150‑1 (MQ=255)                                                                                                                                          ccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAACTTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa        >  2:4836195/1‑150 (MQ=255)                                                                                                                                          ccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCt                          >  1:3485732/1‑132 (MQ=255)                                                                                                                                          ccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCt                          <  2:3485732/132‑1 (MQ=255)                                                                                                                                          ccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa        <  2:530766/150‑1 (MQ=255)                                                                                                                                          ccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa        >  1:395667/1‑150 (MQ=255)                                                                                                                                          ccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa        >  2:2085771/1‑150 (MQ=255)                                                                                                                                          ccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa        >  1:708421/1‑150 (MQ=255)                                                                                                                                          ccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa        <  1:32657/150‑1 (MQ=255)                                                                                                                                          ccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa        >  1:3965354/1‑150 (MQ=255)                                                                                                                                          ccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa        >  1:3965116/1‑150 (MQ=255)                                                                                                                                          ccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAAGTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa        >  1:2575703/1‑150 (MQ=255)                                                                                                                                           gatcgaGCGCCGCCCGCGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCGCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAg       >  1:1338422/6‑150 (MQ=255)                                                                                                                                           cgcgAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCACACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAg       >  1:1420974/1‑150 (MQ=255)                                                                                                                                            gcgAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGa      <  2:922926/150‑1 (MQ=255)                                                                                                                                            gcgAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGa      >  2:1969038/1‑150 (MQ=255)                                                                                                                                            gcgAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGa      >  1:4579773/1‑150 (MQ=255)                                                                                                                                            atcgaGCGCCGCCCGCGCGGCGCTTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAg       >  2:530041/5‑150 (MQ=39)                                                                                                                                             cgAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGAc     >  1:1097975/1‑150 (MQ=255)                                                                                                                                                aGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCt  <  1:2413025/150‑1 (MQ=255)                                                                                                                                                aGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCt  <  2:2480708/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                       ATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCT  >  NC_002947/4808083‑4808375

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑