Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,697,798	0	C	T	25.4%	103.2 / 11.2	55	F447F (TTC→TTT)	bcsE	c‑di‑GMP‑binding protein BcsE
Reads supporting (aligned to +/- strand):  ref base C (22/19);  new base T (13/1);  total (35/20)
Fisher's exact test for biased strand distribution p-value = 9.81e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGTTTCCAGAAGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATC  >  NC_000913/3697550‑3698037                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     ggTTTCCAGAAGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTcc                                                                                                                                                                                                                                                     <  1:95976/250‑1 (MQ=255) ggTTTCCAGAAGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTcc                                                                                                                                                                                                                                                     <  1:256431/250‑1 (MQ=255)                   gtgtTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGGGCCGTTTCTATCATTCT‑GTCGGAATAACGAT‑CTTGAAA‑CC‑AGG‑TT‑GAATCCTATTTTTCCATTTGCTAC‑TGGCGAcct                                                                                                                                                                                                                                 >  1:448223/1‑249 (MQ=255)                             aTGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATAGGCT‑GGGGCTGATTTTCTCCTTTT‑GTCGGGTTAACGAT‑CCGGATA‑CC‑GCG‑TT‑GAATCAAATTTTCCCATTGCCTAC‑TGGCGACATTTTTTC‑A‑AAc                                                                                                                                                                                                                     >  1:204632/1‑251 (MQ=255)                               gccttCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AAcc                                                                                                                                                                                                                    <  1:349627/246‑1 (MQ=255)                                ccGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGGGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTTGATA‑CCCGGG‑TT‑GA‑TCATATTTTTCNCTTGCCTAC‑TCGGGACATTTTTTC‑A‑AACCgg                                                                                                                                                                                                                  >  2:226641/1‑250 (MQ=255)                                ccGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCATACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGGGCTGTTTCCTTCCTTCT‑GTTGGGTTTACGAT‑CTGGGCA‑CC‑GCG‑TG‑GAATTCTAATTTCCCACTGCCTAC‑AGGCGACCtttttctaaacagt                                                                                                                                                                                                                    >  2:419705/1‑241 (MQ=255)                                      acaacaTGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCACTACGTCCGGAACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCC‑GGTGGTGTTTCTCTCATTCT‑TTCGGATTAACGTT‑CTGGGTAACC‑G‑G‑AT‑GGATCATATTTTTCCCTTGCCTTA‑AGGCGACATTTTTTC‑C‑AAACG‑TAT‑ggat                                                                                                                                                                                                          >  1:198320/1‑249 (MQ=255)                                        aacaTGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCCGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑t                                                                                                                                                                                                                                         <  1:267615/222‑1 (MQ=255)                                           atgatgAATAACCCTCTATTACCTGCCCCCGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGGAATCGGCT‑GGTGCCGTTTCTCTCCATCT‑GTTGGATTAACGAT‑CTGGGAA‑CC‑GCG‑TT‑GAATCCTATTTTTCCCTTTACTAC‑TNGCGCCCTTTTTTC‑C‑AACCC‑TAT‑TGTCttggt                                                                                                                                                                                                     >  2:214358/1‑247 (MQ=255)                                             gatgaATAACCCTCGATTCCCTGCCCACGGTAAATGCGGTCGGGTTGCCTGACGTCCGGTACCGGGTACCCGCGTTGACAAAGCCCTGACGCTGTGTCGTCCTAACCGTTCCGGCGATATCATGACCAGTGGCGGTAATCGGCT‑TGTGCTGTTGCTCTCATTCT‑GTCTGATTAACGAT‑CTGGATA‑CC‑GCG‑TG‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTTGTTTg                                                                                                                                                                                                   <  2:470383/251‑1 (MQ=255)                                              atgaATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGGGCTGTTTCTCTCATTCT‑GTCGGGTTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAAACCTATTTTCCCATTGCCTTC‑TGGCGACATTTTCTC‑A‑AAACG‑TAT‑GGTCTGGGTTg                                                                                                                                                                                                   >  1:104159/1‑250 (MQ=255)                                               tgaATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGGAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGGTTAAAGAT‑CTGGGTA‑AC‑GCG‑GT‑GAATCAAATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGa                                                                                                                                                                                                  >  1:43863/1‑250 (MQ=255)                                                taaTAACCCCCTATTTCCGGCCCCCGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTCCCGGCGATATCCTGACCATTGGCGGTAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGaa                                                                                                                                                                                                 <  2:192122/249‑1 (MQ=255)                                                    aaCCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGCACTCGGCT‑GGTGCTGTTTCTCTCATTCT‑TTCGGGTTTACGAA‑CTGGATA‑CC‑GCT‑TT‑GAATAATATTTTCCCCTTTCCTAC‑TGGTGACATTTTTTC‑C‑CACCG‑‑ATAGGTCTGGTTTGAAgatgt                                                                                                                                                                                            >  2:133864/1‑250 (MQ=255)                                                     ccccTCTATTCCCTACCCCCGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTCCCGGGTATCCGCGTTGAACAAGCCCTGACGCTTTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAgatgat                                                                                                                                                                                           <  2:8079/250‑1 (MQ=255)                                                      cccTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTCCCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCTTGGTG‑TTTTTCTCTCATTCT‑CTCGGTTTTACGAT‑CTGGTTC‑CC‑CCG‑TTTGAAT‑ATATTTTTCCATTTCATAC‑TGGCGGCATTTTTTC‑C‑AACCT‑TAT‑GGTCTGGtttttaactgctt                                                                                                                                                                                          >  2:43036/1‑241 (MQ=255)                                                      cccTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGGGATGTTTTTCTCATTCC‑GTCGGGTTTAAGAT‑CTGGGTA‑CC‑GCG‑GT‑GGATCCTATTTTTCCCTTGCCCAC‑TGGCGGCATTTTCTC‑A‑AACCG‑TTT‑GGTtttgtttggaggtgttt                                                                                                                                                                                          >  2:55985/1‑242 (MQ=255)                                                       ccTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTTCCATTGCCTAC‑TGGCGACATTTTTTC‑A‑AACCC‑TAT‑GGTCTGGTTTTAAGATGATc                                                                                                                                                                                          >  1:388253/1‑250 (MQ=255)                                                         tctATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTACTAGGAT‑GGGGCTGTTTTTCTCATTCT‑TTCGGATTAAAGATACTGG‑TA‑CC‑GCG‑TT‑GAATCACATTTTCCCACTTCCTAC‑TGGCGCCATTTTTTC‑A‑AACCG‑TAT‑GGGCTGGGTTTAAgaggattca                                                                                                                                                                                        >  2:376970/1‑247 (MQ=255)                                                           tATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGGGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTTCCTAC‑TGGCGGCATTTTCTC‑A‑AACCG‑TAT‑GGTTTGGTTTTCAGATGATCAAATc                                                                                                                                                                                     >  1:473369/1‑251 (MQ=255)                                                            attgccTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATc                                                                                                                                                                                     <  2:202927/246‑1 (MQ=255)                                                            aTTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAgatga                                                                                                                                                                                            <  1:169509/243‑1 (MQ=255)                                                            aTTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGTGCTGGTTCTCTCATTTT‑GTCGGGTTAACGAT‑CTGGATA‑AC‑GCG‑GT‑GGATCATATTTTTCCATTGCCTAC‑TGGCGCCATTTTTTC‑C‑AACCGGTTT‑GG‑CTGGTTTGAAgatgaagatcggg                                                                                                                                                                                    >  2:169509/1‑243 (MQ=255)                                                               aCCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGCTATCATGACCATTGGCGGTAATCGGCT‑GGGGGTGTTTTTCTCATTCTTGTGGGATTA‑CGAC‑CTTGATA‑CC‑GCGGTT‑TAA‑AATATTTTTCCATTGCCTAC‑GGGCGGCGATTTC‑C‑ACAACCG‑TTT‑TGTCTCGTATGAAGATGTTCAAATccctg                                                                                                                                                                                 >  1:153585/1‑247 (MQ=255)                                                                    cccACGGTAAAGGCGTTCTGGTTGCCCTTCGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGAGCAAATCAGTGCCGAg                                                                                                                                                                            <  2:270575/251‑1 (MQ=255)                                                                    cccACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGTGCGGTTTTTCTCTTTCT‑GTCGGATTGACGAG‑CTGGATA‑CC‑GCG‑TT‑GAATAAAATTTTCCCAATACCTAC‑TGGCCAAATTTTCTC‑A‑AACCG‑TAT‑GGTATGGTTTTAACATGATCAAATTAGTGCCggg                                                                                                                                                                            >  1:18479/1‑249 (MQ=255)                                                                            aaaGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCCGCT‑GGTGCTGTTTCTCTCCTTTT‑TTCCGATTGACGAT‑CTGGATA‑CC‑GCG‑GT‑GGATCATATTTTTCCATTTCCTAC‑TTGCGGCATTTTTTC‑A‑AACCC‑GCT‑GGTCTGGTTTTAAGATGATCAACACAGTGCCGAGCTggggc                                                                                                                                                                     >  2:430110/1‑247 (MQ=255)                                                                                 cGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑C‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGGAGATgc                                                                                                                                                               >  1:77307/1‑251 (MQ=255)                                                                                        gTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGTGCTGTTTTTCTCATTCT‑GTCGGATTAAAGAT‑CTGGAAA‑CC‑GCG‑TT‑GAATCATATTTTTCCCTTTCCTAC‑TGGCGACATTTTCTC‑A‑AAACG‑TAT‑GGGCTGGGTTTAAGATGATCAAATCAGTGCCCGGCTGGTGCAGATTCgcctgct                                                                                                                                                        >  1:140980/1‑251 (MQ=255)                                                                                                         cccGGGTATCCGCGTTGAACAAGTCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATggg                                                                                                                                       <  2:110957/250‑1 (MQ=255)                                                                                                         aCCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGGCTGGTTTTAAAAAGATCAAATCAGTGCCGAGCTGGTGCAGGAGCGCTTGCCTGCCCCCGAACAATggg                                                                                                                                       >  1:110957/1‑251 (MQ=255)                                                                                                             ggTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGGGGTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑GT‑GAATCATTTTTTCCCATTGCCTAC‑TGGCGACATTTTTTC‑A‑AACCG‑TAT‑GGTCTGGGTTGAAGATGGTCAAAACAGAGGCGGGCTGGGGCCGATGCGCTTGGTTGGCCCAGAACAAAGGGGCAt                                                                                                                                   >  1:89674/1‑251 (MQ=255)                                                                                                              gTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATg                                                                                                                                  <  2:159553/251‑1 (MQ=255)                                                                                                               tATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGCAATCGGAT‑GGTGCTGTTTTTCTCCTTTT‑GTCGGATTAACGGT‑CTGGATA‑CC‑GCG‑TT‑GAATTCTATTTTTCCATTGTCTTC‑TGGGGACATTTTCTCCA‑AACCG‑‑AT‑TGTCTGGTTTGAAGATTAACAAATAAGTGCCGATCTGGTGCCGAAGCGCCTTCTTGccccccaaaactggggcactc                                                                                                                                 >  2:465519/1‑234 (MQ=255)                                                                                                                   cgcgTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGCCCATTGGCGGTAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCt                                                                                                                             <  1:425121/251‑1 (MQ=255)                                                                                                                                    gACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGt                                                                                                             <  1:358941/250‑1 (MQ=255)                                                                                                                                            gtCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGTGGTGTTTCTCTCATTCT‑GTCGGGTTAAAGCA‑CTTGATA‑AC‑GCG‑TT‑GGATCCTATTTTCCCCTTTCCTAC‑TGGCGACATTTTCTC‑C‑AACCG‑GAT‑GGTCTGGTTTGAAGATGATCAAATCCGTGCCGAGGCGGAGCAGATGCGCGTGCTTGCCCCAGAACAAAGGGGCATGCCGCTGCCTTTAACGCCAAGGTCTaaaacc                                                                                                    >  1:2067/1‑248 (MQ=255)                                                                                                                                                            cGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGTGCTGTTTCTCTCATTAC‑GTCGGGATAATGAT‑CTGGTTA‑CC‑GCG‑TT‑GAAACATATTTTTCCCTTTCCTACCTGG‑GGCATTTTCTC‑A‑CACCG‑GAT‑GGGCTGGTTTGAAGTTGATCCAATCAGAGCCGAGATGGGGCAGATGAGCTTGCTTGCCACAGAAACATAGGGCATGACGCTGCCTTTAACGCAAAGGTCTAAACCAGtccttcaatcccagg                                                                                    >  2:379636/1‑238 (MQ=255)                                                                                                                                                               cGATATCATGACCATTGGCGGTAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAAc                                                                                                      <  1:169195/230‑1 (MQ=255)                                                                                                                                                               cGATATCATGACCATTGGCGGTAATCGGCT‑GGGGCTGTTTCTCTCATTTT‑GTTGGATTAACGAT‑CTGGATA‑AC‑GCG‑GT‑GGATCATTTTTTCCCATTTCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGGAGATGATCAAATAAGTGCCGAGATGGGGCAGATGCGCTTGCTTGCCCCCGAACAATGGGGCATGCCGCTGCCTTTAACGCACAGGTCTAAAcagctctgaagggagtcgggtt                                                                                 >  2:169195/1‑230 (MQ=255)                                                                                                                                                                            aTTGGCGGTAATCGGCT‑GGTGCTGTTTCTCTCAGTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTg                                                                    <  1:242654/251‑1 (MQ=255)                                                                                                                                                                                 cGGTAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGcgac                                                               >  1:62830/1‑251 (MQ=255)                                                                                                                                                                                   gTAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGcgacga                                                             <  1:278725/251‑1 (MQ=255)                                                                                                                                                                                                            cATTCT‑GTCGGATTAACGGT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGGCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGt                                       <  2:78211/249‑1 (MQ=255)                                                                                                                                                                                                                     cGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGc                             >  1:430419/1‑251 (MQ=255)                                                                                                                                                                                                                      ggATTAACGAT‑CTGGATA‑CC‑GAG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGACGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACTAGAACCCGTGCGACTGTGAGATGATGCt                            >  1:296870/1‑251 (MQ=255)                                                                                                                                                                                                                          tAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTgg                        >  1:352034/1‑251 (MQ=255)                                                                                                                                                                                                                           aaCGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTgg                        >  1:484153/1‑250 (MQ=255)                                                                                                                                                                                                                               aT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAgcgc                   >  1:17335/1‑251 (MQ=255)                                                                                                                                                                                                                                    ggATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCtc                 <  2:62830/249‑1 (MQ=255)                                                                                                                                                                                                                                          cc‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCatg           >  1:64830/1‑250 (MQ=255)                                                                                                                                                                                                                                           c‑gcg‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TCT‑GGTCTGGTTTGCAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCatga          >  2:385015/1‑250 (MQ=255)                                                                                                                                                                                                                                             ccccgtggAATCATATTTCTCCATGGCATAC‑TGGCGACATTTACTC‑A‑AGCCG‑TAT‑GGGCTGGGTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACTATTGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATTCTGTGTAGCGCtaatcacga          <  2:33859/244‑4 (MQ=255)                                                                                                                                                                                                                                                    gAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCAt   >  1:4775/1‑251 (MQ=255)                                                                                                                                                                                                                                                    gAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCAt   <  1:483408/251‑1 (MQ=255)                                                                                                                                                                                                                                                     aaTCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGATGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATc  <  1:295079/251‑1 (MQ=255)                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     GGTTTCCAGAAGTGGGATGTGTTCTGTAATGCCGTCAACAACATGATGAATAACCCTCTATTACCTGCCCACGGTAAAGGCGTTCTGGTTGCCCTACGTCCGGTACCGGGTATCCGCGTTGAACAAGCCCTGACGCTGTGTCGCCCTAACCGTACCGGCGATATCATGACCATTGGCGGTAATCGGCT‑GGTGCTGTTTCTCTCATTCT‑GTCGGATTAACGAT‑CTGGATA‑CC‑GCG‑TT‑GAATCATATTTTCCCATTGCCTAC‑TGGCGACATTTTCTC‑A‑AACCG‑TAT‑GGTCTGGTTTGAAGATGATCAAATCAGTGCCGAGCTGGTGCAGATGCGCTTGCTTGCCCCAGAACAATGGGGCATGCCGCTGCCTTTAACGCAAAGTTCTAAACCGGTCATCAATGCCGAGCACGATGGTCGCCACTGGCGACGAATACCAGAACCCATGCGACTGTTAGATGATGCTGTGGAGCGCTCATCATGATGACCATC  >  NC_000913/3697550‑3698037

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