Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779_cfb	2,004,055	0	A	C	28.6%	88.0 / 16.3	49	intergenic (‑129/+192)	fliA/fliC	RNA polymerase, sigma 28 (sigma F) factor/flagellar filament structural protein
Reads supporting (aligned to +/- strand):  ref base A (32/3);  new base C (0/14);  total (32/17)
Fisher's exact test for biased strand distribution p-value = 1.01e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAAC  >  NC_007779_cfb/2003809‑2004301                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        cccnccnggcgGGCAGGCGAACCGGCAGGCGAAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTACAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAAcc‑c                                                                                                                                                                                                                                                     <  2:1098679/242‑1 (MQ=255)  cACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGTGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAAcc‑cc                                                                                                                                                                                                                                                    >  1:47211/1‑250 (MQ=255)        cgcgGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGACTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGt                                                                                                                                                                                                                                              >  1:577005/1‑250 (MQ=255)         gcgGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTg                                                                                                                                                                                                                                             >  1:913676/1‑250 (MQ=255)             gCAGTCGACCCGGGAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAGCGCGGCCACAGCGAGTGTTTATCCATTCCGCCGTCAGCGGTATAGGGTGAATTCACGATAAACAGCCAGGCGTTATATGTGTTACCGGCATTATTATCCGTTTCTGCAGGGTTTTAAATCGGACGATTAGGGGGAAAAAGGGGGGGTTATTTGGGGGTTACGGGAAAATTCC‑GGGGAAAAAAAACCC‑CCGCCGGTGGCgt                                                                                                                                                                                                                                         <  2:81578/250‑2 (MQ=255)                        tgaaGGCGCAGTGCTTGGGGACGAACCAGCGGGAAATAACGTTGCCACAGGGAGTGTTTATCCATTACCCCTCCAGCGGTATAGAGTGATTTCACGAAAAACGGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCGGGGTTTTTAATCGGACGATTGGTGGGGAAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAACCC‑CCGCCGGTGGCGGGGAAGCACGtt                                                                                                                                                                                                                              <  2:1279364/247‑1 (MQ=255)                        tGCGGGCGNAATGCTTCTTGACGAGCCGGCGGGACATAACGTTGCGCCAGCGAGTGTTTATCCATGCCACCTGAAGCGGTAGAGAGTGAATTCACTAGAACCAGCCCTGCGTTATATGAGTTACCGGCATGATTAGCCGTTTCTGCAGGTTTTTTAATCGGAGGATTAGTGGGTGAAGTGGGGGTTTATTTGGGGGTTACAGGTAATTTCC‑AGCCGGAAAAAACCC‑CCGCCGGTGGCGGGGAAGCACGtt                                                                                                                                                                                                                              <  2:656093/250‑1 (MQ=255)                                      ttCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGGAGAACGTTTCTGACAAATTtcg                                                                                                                                                                                                                >  2:178720/1‑250 (MQ=255)                                                  cAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGcc                                                                                                                                                                                                    >  2:295134/1‑250 (MQ=255)                                                     tgggACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGCGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGa                                                                                                                                                                                                 >  1:1257196/2‑250 (MQ=255)                                                                                      cATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACaaaaa                                                                                                                                                                >  1:230909/1‑250 (MQ=255)                                                                                       aTTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAAt                                                                                                                                                               >  1:366230/1‑250 (MQ=255)                                                                                           cagcTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGCTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCAAGTTGCTGACAAATTGCGCTTTATGTTTCCGGATGCGGCGTAAACGCCTTATCCGGGCTACAAAAAtgtcc                                                                                                                                                           >  2:1157158/4‑248 (MQ=255)                                                                                                 cAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTc                                                                                                                                                     >  1:1270091/1‑250 (MQ=255)                                                                                                         taGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAAtt                                                                                                                                             >  1:1024770/1‑250 (MQ=255)                                                                                                          ngattnatttCGCGATAAACAGCCCTGCTTTATATGGGTTATCGGCATGATTAGCCTTTTCTGCGGGGTTTTTAATGGGACGATTGGTGGGTAAAAGGAGGGTTTTTTTGGGGGTTAGAGTTAAATCCG‑AGGAAAAAAAAACCC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                                                                                                                            <  2:785939/242‑1 (MQ=37)                                                                                                          ggagTGATTCCACGATAAACTGCCCGGCGTTATATGGGTTTGCGGCTTGATTCTCCGTTCCGGCAGTGTTTTTAATCGGACGATTATTGGGAAAAATGAGGGTTTTTTTGGGGGTAACAGATAAATCCC‑AGAAAAAAAAAACCC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTg                                                                                                                                            <  2:230909/249‑1 (MQ=37)                                                                                                                     cagATAAAAAGCCCTGCGTTATATGAGTTGTCGATATAATTCTCTGTTGTTGGGGTGTTTTTAAGGGGACGATTAGGGGTTAAAAGGAGGGTTTTTTTGGGGGTAAAAGGAAAATCCA‑AGGAAAAAAAAACCC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACtt                                                                                                                                 <  2:1270091/248‑1 (MQ=25)                                                                                                                        aTAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTa                                                                                                                              >  1:1052572/1‑250 (MQ=255)                                                                                                                         aaaaCAGCCTTGCGTTATAGGAGTGAGCGGCATGATTATCCGTTTCTGCGGGGTTTTTAATCGGACGATTAGTGGGGAAAAGGAGGGGTTTTTTGGGGGTTACAGAAAAATCCC‑AGGAAAAAAAAACCC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAg                                                                                                                             <  2:1024770/249‑1 (MQ=38)                                                                                                                               tttcgtttctataGGATTTATCGGCATTATTTCCCTTTTCTGCAGGTTTTTAAACGGGGCGTTTAGTGGGTAAAATGAGGGTTTTTTTGGGGGTTACAGG‑AAATTCCAAGGCAAAAAAAACCC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGa                                                                                                                       <  2:495771/241‑1 (MQ=25)                                                                                                                                         atatGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGaagatcggaag                                                                                                             >  1:495771/1‑240 (MQ=255)                                                                                                                                                ttATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGAGGGTGAAATGAGGGGGTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAcgcgcagcgctacag                                                                                                      >  2:1034922/1‑245 (MQ=255)                                                                                                                                                       cATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGCTATTTGGGGGGTACTGGTAAATTCC‑CGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGTGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGAAAttt                                                                                               >  2:226847/1‑250 (MQ=255)                                                                                                                                                            ttATCCGTCTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGt                                                                                          >  1:1195141/1‑250 (MQ=255)                                                                                                                                                                  gTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATGAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGt                                                                                    >  1:280734/1‑250 (MQ=255)                                                                                                                                                                      cTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGATAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGt                                                                                >  1:594314/1‑250 (MQ=255)                                                                                                                                                                             gTTTTTAATCGGACGATAGGGGGGGAAAAGGGGGGGTTATTGGGGGGTTACAGGAAAATTCC‑AGGAAGAAAAAACCC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCGGTCAGTCTCAGtt                                                                         <  1:178720/250‑1 (MQ=37)                                                                                                                                                                               ttttAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTaa                                                                       >  1:1080100/1‑250 (MQ=255)                                                                                                                                                                                   aaTCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAg                                                                   <  1:295134/250‑1 (MQ=255)                                                                                                                                                                                     tCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGt                                                                 >  1:208448/1‑250 (MQ=255)                                                                                                                                                                                      cGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAg                                                                           <  1:80419/239‑1 (MQ=255)                                                                                                                                                                                      cGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAg                                                                           >  2:80419/1‑239 (MQ=255)                                                                                                                                                                                        gACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTAc                                                              >  1:265739/1‑250 (MQ=255)                                                                                                                                                                                            aTTAGTGGGTGAAATGAGGGGTTATTTGAGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACg                                                          >  1:1293626/1‑250 (MQ=255)                                                                                                                                                                                              tgggggnggAAATGGGGGTTTTTTTTGGGGTTAAAAGTAAAATCC‑AAGAAAAAAAAAACCACCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGa                                                         <  2:1052572/242‑1 (MQ=37)                                                                                                                                                                                              aggTGGGTGAAAGGGGTGTTTTTTTGGGTGTACCAGTAAAATCCC‑AGGAAAAAAAAACCC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATGAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGAt                                                        <  2:280734/248‑1 (MQ=37)                                                                                                                                                                                                    gTGAAATGAGGGGGTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACCAATTGCGCTTTATGTTGCCGGGTGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGAGAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTCACACGATTAAccc                                                  >  2:1285990/1‑250 (MQ=255)                                                                                                                                                                                                         aTGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCGGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTCGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTgcag                                             >  1:115032/1‑250 (MQ=255)                                                                                                                                                                                                                tttATTTGGGGGTAAAAGTAAAATCCC‑AGGAAAAAAAAACCC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGAc                                      <  2:366230/249‑1 (MQ=38)                                                                                                                                                                                                                  tttttGGGGGTAACAGAAAAATACA‑GGGAAAAAAAAACCC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAg                                    <  2:208448/248‑1 (MQ=38)                                                                                                                                                                                                                  tATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAg                                    >  1:919632/1‑250 (MQ=255)                                                                                                                                                                                                                  tATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTCGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCACAGTTAATCAGTTTACAACGATTAACCCTGCCGCAGAGACAg                                    >  2:464262/1‑250 (MQ=255)                                                                                                                                                                                                                        ggggTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACctg                              >  1:1016648/1‑250 (MQ=255)                                                                                                                                                                                                                        aggctataaagTAAATACG‑AGGAAGAAAAACTCC‑CCGCCGGTGGCGGGTAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACctg                              <  2:1293626/240‑1 (MQ=38)                                                                                                                                                                                                                          ggTTAGAGGTAAATCCC‑AGGAAAAAAAAACCC‑CCGCCGGTGGCGGGGAAGCACGTCGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACctgct                            <  1:464262/250‑1 (MQ=255)                                                                                                                                                                                                                                  gTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGATGCTGACAAATTTCGCTTTATGTTGCCGGACGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATCAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTAcc                    >  2:981109/1‑250 (MQ=255)                                                                                                                                                                                                                                     aaTTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTTCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACGGAACCTGCTGCGGTACCTgg                 >  2:214291/1‑250 (MQ=255)                                                                                                                                                                                                                                      ctacc‑AGGAAAAAAAACCCC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGGCAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGt                <  2:1332969/247‑1 (MQ=255)                                                                                                                                                                                                                                        tCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGGCAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGtag              >  1:1332969/1‑248 (MQ=255)                                                                                                                                                                                                                                          ccaaganaaaaacacCC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAg             <  2:1080100/242‑1 (MQ=255)                                                                                                                                                                                                                                            aGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCGGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCt           >  1:181864/1‑250 (MQ=255)                                                                                                                                                                                                                                                caaaaaaacCCGCCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCGGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCtttt        <  2:181864/248‑1 (MQ=255)                                                                                                                                                                                                                                                    aaaaaCC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGc            >  2:292545/1‑241 (MQ=255)                                                                                                                                                                                                                                                    aaaaaCC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGc            <  1:292545/241‑1 (MQ=255)                                                                                                                                                                                                                                                     aaaaCC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCtttt        >  2:377735/1‑244 (MQ=255)                                                                                                                                                                                                                                                     aaaaCC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCtttt        <  1:377735/244‑1 (MQ=255)                                                                                                                                                                                                                                                     aaaaCC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAac  >  1:368638/1‑250 (MQ=255)                                                                                                                                                                                                                                                        aCC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGt                       >  2:393520/1‑226 (MQ=255)                                                                                                                                                                                                                                                        aCC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGt                       <  1:393520/226‑1 (MQ=255)                                                                                                                                                                                                                                                        aCC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCt           >  1:1226956/1‑238 (MQ=255)                                                                                                                                                                                                                                                        aCC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCt           <  2:1226956/238‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        CCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGCAGGGTTTTTAATCGGACGATTAGTGGGTGAAATGAGGGGTTATTTGGGGGTTACAGGTAAATTCC‑AGGCAGAAAAAAACC‑CCGCCGGTGGCGGGGAAGCACGTTGCTGACAAATTGCGCTTTATGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACAAAAATGTGCAAATTCAATAAATTGCAATTCAACTTGTAGGCCTGATAAGCGCAGCGCATCAGGCAATTTGGCGTTGCCGTCAGTCTCAGTTAATCAGGTTACAACGATTAACCCTGCAGCAGAGACAGAACCTGCTGCGGTACCTGGTTAGCTTTTGCCAAC  >  NC_007779_cfb/2003809‑2004301

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