Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,808,232	0	G	T	34.5%	58.5 / 54.0	55	intergenic (‑296/‑259)	PP_4234/dsbD‑II	bifunctional aldehyde oxidase/xanthine dehydrogenase/thiol/disulfide interchange protein 2
Reads supporting (aligned to +/- strand):  ref base G (16/20);  new base T (7/12);  total (23/32)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAA  >  NC_002947/4808087‑4808380                                                                                                                                                  |                                                                                                                                                      gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCttcgcg                                                                                                                                                   >  1:875378/1‑150 (MQ=255) gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCttcgcg                                                                                                                                                   <  1:2646854/150‑1 (MQ=255)  cGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgg                                                                                                                                                  <  1:1320825/150‑1 (MQ=255)  cGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgg                                                                                                                                                  <  2:5183081/150‑1 (MQ=255)  cGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgg                                                                                                                                                  <  2:1502812/150‑1 (MQ=255)  cGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgg                                                                                                                                                  <  1:1414711/150‑1 (MQ=255)  cGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgg                                                                                                                                                  <  1:3439994/150‑1 (MQ=255)  cGGAATACGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgg                                                                                                                                                  <  2:1299940/150‑1 (MQ=38)   ggAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcggg                                                                                                                                                 >  1:4933078/1‑150 (MQ=255)    gAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggt                                                                                                                                                <  1:4047152/150‑2 (MQ=255)    gAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggt                                                                                                                                                >  1:2006748/1‑149 (MQ=255)        tCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGTTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGAACGAGCGACGCACGCGCGGCGCTTCgcgggtaaac                                                                                                                                            >  1:4845310/1‑145 (MQ=37)        tCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggtaaac                                                                                                                                            <  1:466819/150‑6 (MQ=38)        tCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggtaaac                                                                                                                                            <  1:3500404/150‑6 (MQ=38)          gcgcTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggtaaaccc                                                                                                                                          <  1:2529105/150‑8 (MQ=38)          gcgcTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggtaaaccc                                                                                                                                          >  1:5278728/1‑143 (MQ=38)          gcgcTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcgggtaaaccc                                                                                                                                          >  1:3160562/1‑143 (MQ=38)                                 gCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcggg                                                                                                                                                 <  2:4605441/120‑1 (MQ=37)                                 gCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcggg                                                                                                                                                 >  1:4605441/1‑120 (MQ=37)                                                                gCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcggg                                                                                                                                                 >  2:4479276/1‑89 (MQ=18)                                                                gCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGCGGCGCTTCgcggg                                                                                                                                                 <  1:4479276/89‑1 (MQ=18)                                                                                                                              gggtttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCca                     <  2:4610903/143‑1 (MQ=255)                                                                                                                              gggtttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCca                     >  2:633489/8‑150 (MQ=255)                                                                                                                               tttttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTAAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa                    <  1:44963/144‑1 (MQ=255)                                                                                                                               ggtttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACCGGCGCTGAAAATCATCCGGCAAACACCACAAAATTCCAACGCCGTCCACTCAAAATTCCAACTCCACCCCCTCAAAATTCCAAAGCCGCCTCCCCTCGAAAATCcaa                    <  1:414302/144‑1 (MQ=38)                                                                                                                               ggtttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa                    <  1:1385034/144‑1 (MQ=255)                                                                                                                               ggtttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa                    <  1:1385261/144‑1 (MQ=255)                                                                                                                               ggtttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa                    <  2:4130403/144‑1 (MQ=255)                                                                                                                               ggtttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa                    <  2:3562735/144‑1 (MQ=255)                                                                                                                               ggtttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCACCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaa                    <  1:3060801/144‑1 (MQ=255)                                                                                                                                gtttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaac                   >  1:2014486/6‑150 (MQ=255)                                                                                                                                gtttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaac                   >  2:3215699/6‑150 (MQ=255)                                                                                                                                  ttacccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCcaacaa                 <  1:207836/147‑1 (MQ=255)                                                                                                                                      ccgcGAAGCGCCGCGTGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa             >  1:3027388/1‑150 (MQ=255)                                                                                                                                      ccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCa                                             <  2:1745122/118‑1 (MQ=255)                                                                                                                                      ccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCa                                             >  1:1745122/1‑118 (MQ=255)                                                                                                                                      ccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGACAATCCAACAACaaa             >  2:792358/1‑150 (MQ=255)                                                                                                                                      ccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa             >  2:579069/1‑150 (MQ=255)                                                                                                                                      ccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa             >  1:4138613/1‑150 (MQ=255)                                                                                                                                      ccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa             <  1:1892584/150‑1 (MQ=255)                                                                                                                                      ccgcGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAAGTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa             >  2:792308/1‑150 (MQ=255)                                                                                                                                       cgcgAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAg            >  2:1701217/1‑150 (MQ=255)                                                                                                                                        gcgAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGa           <  2:4275630/150‑1 (MQ=255)                                                                                                                                        gcgAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGa           >  1:4646725/1‑150 (MQ=255)                                                                                                                                        atcaaGCGCCGCGCGGGCGGCGC‑TTGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGa           >  2:3702872/4‑150 (MQ=255)                                                                                                                                           gaGCGCCGCCCGCGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa             >  1:2493303/2‑145 (MQ=255)                                                                                                                                           gaGCGCCGCCCGCGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa             <  2:2493303/144‑1 (MQ=255)                                                                                                                                           gaGCGCCGCCCGCGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGAccc        >  2:5063454/2‑150 (MQ=255)                                                                                                                                           aaGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGAccc        <  1:2309079/150‑1 (MQ=255)                                                                                                                                            aGCGCCGCGCGGGCGGCGC‑TCTATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCt       <  1:4011462/150‑1 (MQ=255)                                                                                                                                            aGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCt       <  2:132253/150‑1 (MQ=255)                                                                                                                                            aGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCt       <  1:1343935/150‑1 (MQ=255)                                                                                                                                            aGCGCCGCCCGCGCGGCGCTTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGAccc        >  1:1434265/1‑150 (MQ=255)                                                                                                                                             gcgcCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTg      <  2:2014486/150‑1 (MQ=255)                                                                                                                                             gcgcCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTg      <  1:324934/150‑1 (MQ=255)                                                                                                                                             gcgcCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTg      <  1:5150661/150‑1 (MQ=255)                                                                                                                                             gcgcCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTg      >  1:1053771/1‑150 (MQ=255)                                                                                                                                              cgccgcCCGCGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTc                                               <  1:658813/108‑1 (MQ=255)                                                                                                                                              cgccgcCCGCGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTc                                               >  2:658813/1‑108 (MQ=255)                                                                                                                                               gccgcGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCa    >  1:774238/1‑150 (MQ=255)                                                                                                                                               gccgcGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCa    >  2:193091/1‑150 (MQ=255)                                                                                                                                                 cgcgcgGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCaaa  <  2:1424/150‑1 (MQ=255)                                                                                                                                                 cgcgcgGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCaaa  <  2:2647367/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                      GCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGC‑TCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAA  >  NC_002947/4808087‑4808380

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