Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007779_cfb	4,371,274	(A)5→3	intergenic (‑9/+296)	dcuA ← / ← aspA	C4‑dicarboxylate antiporter/aspartate ammonia‑lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779_cfb	4,371,271	0	A	.	95.1%	147.7 / 2.8	41	intergenic (‑6/+300)	dcuA/aspA	C4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base . (32/7);  total (32/9)
Fisher's exact test for biased strand distribution p-value = 4.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01
*	NC_007779_cfb	4,371,272	0	A	.	95.1%	147.8 / 2.8	41	intergenic (‑7/+299)	dcuA/aspA	C4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base . (32/7);  total (32/9)
Fisher's exact test for biased strand distribution p-value = 4.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTG  >  NC_007779_cfb/4371027‑4371510                                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                                               gTGATCGC‑ATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaatagc                                                                                                                                                                                                                                           <  1:1430076/250‑5 (MQ=255)       gCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGa                                                                                                                                                                                                                                    <  2:109466/250‑1 (MQ=255)         aaTCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATa                                                                                                                                                                                                                                  <  1:1355947/250‑1 (MQ=255)            cTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAAc                                                                                                                                                                                                                                >  1:1375160/1‑249 (MQ=255)            cTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAAc                                                                                                                                                                                                                                <  2:1375160/249‑1 (MQ=255)              gcgcGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGc                                                                                                                                                                                                                             <  2:1474389/250‑1 (MQ=255)                   gATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCGGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAACGATAGGCAGGAAATAGGTCACGATCGGGGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGGTGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGa                                                                                                                                                                                                                        >  1:1180563/1‑250 (MQ=255)                       ccACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaat                                                                                                                                                                                                                                              >  1:1220257/1‑225 (MQ=255)                       ccACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaag                                                                                                                                                                                                                                              <  2:1220257/226‑2 (MQ=255)                             cAGTAGAAAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCTTGATNATGGAGAt                                                                                                                                                                                                              >  2:791419/1‑250 (MQ=255)                              aGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATg                                                                                                                                                                                                             >  1:116531/1‑250 (MQ=255)                                   aCAGGGGAGGACATGGTTAACCGCCTTGTCCCTTCGAAACTCGAGCGAGAAGGGGGGGTGGCGCCAGAGAGAGGTGGCCGGTCCCAGAAAAGATAGTCAGAAAATAGGTCACGAGCGGTGCGAGGATCGTGATGTATTCCGGGCCGCGGCGCAGAAGCTTTTCTGTCGGATAACCCAGATAGCCCAGCCCCCCAGAAACCTGAAGGGCAGAAAA‑‑AGCGGCGATAACCGCCATGATAATGGAGATGACATc                                                                                                                                                                                                        <  2:277078/250‑1 (MQ=255)                                                aGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATgaga                                                                                                                                                                                           >  1:1470530/1‑247 (MQ=255)                                                         gCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGAt                                                                                                                                                                                                                                   <  1:668065/201‑1 (MQ=255)                                                         gCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGAt                                                                                                                                                                                                                                   >  2:668065/1‑201 (MQ=255)                                                          ccTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTaa                                                                                                                                                                                 >  1:1029140/1‑250 (MQ=255)                                                               ttCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaata                                                                                                                                                                                                                                             >  2:1541733/1‑185 (MQ=255)                                                               ttCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaata                                                                                                                                                                                                                                             <  1:1541733/187‑3 (MQ=255)                                                               ttCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCa                                                                                                                                                                            >  1:777091/1‑250 (MQ=255)                                                                tCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCaa                                                                                                                                                                           >  1:1154268/1‑250 (MQ=255)                                                                  cTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAAGCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGGCATCGAACGGGATGTTACCGGGTTTAACGCCAATa                                                                                                                                                                         >  1:871008/1‑250 (MQ=255)                                                                  cTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATa                                                                                                                                                                         >  1:1582781/1‑250 (MQ=255)                                                                         aCTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCaa                                                                                                                                                                  >  1:1457836/1‑250 (MQ=255)                                                                           ttCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCCGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGa                                                                                                                                                                >  1:155190/1‑250 (MQ=255)                                                                                    aaCTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAAcc                                                                                                                                                                                                                               <  1:1336479/178‑1 (MQ=255)                                                                                    aaCTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAAcc                                                                                                                                                                                                                               >  2:1336479/1‑178 (MQ=255)                                                                                        ggCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGATATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCa                                                                                                                                                   >  1:405264/1‑250 (MQ=255)                                                                                        ggCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCa                                                                                                                                                   >  1:132687/1‑250 (MQ=255)                                                                                          cAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATATTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaata                                                                                                                                                                                                                                             >  2:903996/1‑158 (MQ=255)                                                                                          cAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaata                                                                                                                                                                                                                                             <  1:903996/160‑3 (MQ=255)                                                                                                cGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCt                                                                                                                                           >  1:307475/1‑250 (MQ=255)                                                                                                 gCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAAc                                                                                                                                                                                                                                <  1:1367961/164‑1 (MQ=255)                                                                                                 gCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAAc                                                                                                                                                                                                                                >  2:1367961/1‑164 (MQ=255)                                                                                                      agagaTGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGCGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGGTAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAAACGCCATGATAATGGAGATGACATCGAACCGGCTGTTACCGGGGTTACGCCCAATAGCGTCAATAACAAGTATCCCCAACCCGCCCGCaaaa                                                                                                                                     >  2:967201/1‑250 (MQ=255)                                                                                                                                                 aCGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCa                                                                                          >  1:448892/1‑250 (MQ=255)                                                                                                                                                      cGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCaaa                                                                                     >  1:878787/1‑250 (MQ=255)                                                                                                                                                                     gATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCt                                                                      >  1:175338/1‑250 (MQ=255)                                                                                                                                                                                   tACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATAtt                                                        >  1:1090230/1‑250 (MQ=255)                                                                                                                                                                                                cTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGtt                                           >  1:361568/1‑250 (MQ=255)                                                                                                                                                                                                 ttttCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGgtt                                          >  1:1576956/1‑250 (MQ=255)                                                                                                                                                                                                      tgtctgATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTNNTNTaa                                     >  2:1056660/1‑250 (MQ=255)                                                                                                                                                                                                        tctgATGAACCAGACAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGACAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACAAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACa                                   >  1:450959/1‑250 (MQ=255)                                                                                                                                                                                                           gatgaACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGt                                >  1:1053812/1‑250 (MQ=255)                                                                                                                                                                                                                            aGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTAt               >  2:1489910/1‑250 (MQ=255)                                                                                                                                                                                                                             gACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTAtt              >  2:850908/1‑250 (MQ=255)                                                                                                                                                                                                                                 gCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTa          >  1:208774/1‑250 (MQ=255)                                                                                                                                                                                                                                      cAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTa     >  1:1149911/1‑250 (MQ=255)                                                                                                                                                                                                                                       aaCCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACa                                   >  2:420490/1‑219 (MQ=255)                                                                                                                                                                                                                                       aaCCTGCATGGC‑‑AAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACa                                   <  1:420490/219‑1 (MQ=255)                                                                                                                                                                                                                                           cctgcagggaaaaaaTGGCGGGAATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTg  <  1:1056660/240‑1 (MQ=255)                                                                                                                                                                                                                                                     ||                                                                                                                                                                                                                                               GTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGGTTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTTGCCAGTACCAGCAAAGATAGTCAGGAAATAGGTCACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTG  >  NC_007779_cfb/4371027‑4371510

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