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	.	94.9%	136.4 / 4.7	40	intergenic (‑6/+300)	dcuA/aspA	C4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base . (33/4);  total (33/7)
Fisher's exact test for biased strand distribution p-value = 2.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01
*	NC_007779_cfb	4,371,272	0	A	.	92.5%	135.5 / 5.7	40	intergenic (‑7/+299)	dcuA/aspA	C4‑dicarboxylate antiporter/aspartate ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base A (0/3);  new base . (33/4);  total (33/7)
Fisher's exact test for biased strand distribution p-value = 3.54e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

CGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGG‑TTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTT‑GCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG  >  NC_007779_cfb/4371025‑4371518                                                                                                                                                                                                                                                          ||                                                                                                                                                                                                                                                       cGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGG‑TTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGACAGTGTCGCCAGAGAGATGTT‑GCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCGGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaat                                                                                                                                                                                                                                                      <  1:1333941/250‑2 (MQ=255)    tGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGG‑TTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGGGATGTT‑GCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaatagc                                                                                                                                                                                                                                                   >  1:981138/1‑246 (MQ=255)      aTCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCCAGG‑TTTAACGCCTTGTTCCTTCGCAACTTCAGCGATACCTGGCATTGCCGCAAGGGGGATTTT‑GCCAGTTCCAGCAACGATATCCGGGAAATAGGT‑CACGAGGGGTGGGAGGAGCGTGAGGTATTTGGGGTTACGGCGCACAAGCTTTTCTGTCTGATGACCCAGTAAGTCCAAACCGCCAGCAACCTGCAGGGC‑‑AAAAATAgcgc                                                                                                                                                                                                                                               >  2:913217/1‑249 (MQ=255)       tCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGG‑TTTAACGCCTTGTTCCTTCGCAACGTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTT‑GCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGaaatagcggc                                                                                                                                                                                                                                                <  2:300535/250‑8 (MQ=255)        cGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGG‑TTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTT‑GCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAgcggcg                                                                                                                                                                                                                                             <  1:913217/250‑1 (MQ=255)          cAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGG‑TTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTT‑GCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTTCATGGCAG‑‑AAATAGCGGCGAt                                                                                                                                                                                                                                           <  2:800635/250‑1 (MQ=255)                                    gACAGCGGACGGCAAGG‑TTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTT‑GCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACAt                                                                                                                                                                                                                 >  1:1094347/1‑250 (MQ=255)                                         tccatatcAAGGTTTTACCGCTTTGTTCCTTCGAAACTTAGGGAAAAACTGGCAGGGTCGCAAGAGAGTTGTG‑GCAAGAACAAGAAAAGATAGCAAGAAAATGGTC‑CACGTTGGGTGCGGGGATCGTGATGTATTGCGGGTTACGGCGCAGCATTTTTTCTGTCTGATAACCCAAATAGCCAAGCCCGCCAGAACCCTGCA‑GGGA‑AAAAATAGCGGCAATAACCGCCATGATAATGGAGATGACATCGaa                                                                                                                                                                                                             <  2:96169/243‑1 (MQ=255)                                          ggACGGCAAGG‑TTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTT‑GCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGTTCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACg                                                                                                                                                                                                           >  1:254042/1‑250 (MQ=255)                                               gCAAGT‑NTTAACGTCTTGTCCCTTCCCAACATCAGAAAAACCTGGGGGGGACGCGGGAGAGATGTG‑GCCGATACAAAAAAAGATAGTCGGAAAATAGGG‑CACGAGGGGGGGGAGGTTCGTTAGGTATTGCGGGTTGGGGCGCAGCAGTTTTTCTGTCTGATAAAACAGATAGCCCAGCCCCCCAGCAACCTGCGTGGC‑‑AAAAAGGGGGGCGAAAGCCGCCATGATAATGGAGATGACATCGAACGGGATg                                                                                                                                                                                                      <  2:1098662/250‑1 (MQ=255)                                                 aaGG‑TTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTT‑GCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGaaagatcggaa                                                                                                                                                                                                    >  1:96169/1‑241 (MQ=255)                                                         aaCGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTT‑GCCAGTACCCGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAgcggcg                                                                                                                                                                                                                                             >  2:967223/1‑202 (MQ=255)                                                         aaCGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTT‑GCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACAGCCAGCAACCTGCATGGCAG‑‑AAATAgcggcg                                                                                                                                                                                                                                             <  1:967223/202‑1 (MQ=255)                                                           cGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTT‑GCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGttt                                                                                                                                                                                           >  1:487419/1‑250 (MQ=255)                                                              cTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTT‑GCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCCTGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAAc                                                                                                                                                                                        >  1:863379/1‑250 (MQ=255)                                                                    ccTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTT‑GCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGAt                                                                                                                                                                                                                                           >  2:639777/1‑193 (MQ=255)                                                                    ccTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTT‑GCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGAt                                                                                                                                                                                                                                           <  1:639777/193‑1 (MQ=255)                                                                            aCTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTT‑GCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCaa                                                                                                                                                                          >  1:1268268/1‑250 (MQ=255)                                                                                aagaaatacgtGGCAGTGCCGCAAGAGAG‑TGTTCGCCTGTCCCGGAAAATATATTCA‑GAAATTGGTCCCCGTGGGGGGCGAGTATGTTGGTTTTTTGCGGGTTACGCGGCAGCAGTTTTTGTTTCTGATAACCCAAATGGCCCACCCCCCCGGAAACCTGCA‑GG‑ACGAAAAGAGGGGCAATAACCGTCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAAc                                                                                                                                                                      <  2:927637/240‑1 (MQ=255)                                                                                                   cGCCAGAGAGATGTT‑GCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCt                                                                                                                                                   >  1:374796/1‑250 (MQ=255)                                                                                                        gagagaTGTT‑GCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCaaa                                                                                                                                              >  1:432916/1‑250 (MQ=255)                                                                                                             aTGTT‑GCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCaa                                                                                                                                         >  1:199773/1‑250 (MQ=255)                                                                                                              tGTT‑GCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAAt                                                                                                                                        >  1:263106/1‑250 (MQ=255)                                                                                                                       gTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTc                                                                                                                                >  1:719574/1‑250 (MQ=255)                                                                                                                         aCCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGAATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTccc                                                                                                                              >  1:39457/1‑250 (MQ=255)                                                                                                                                aaGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATc                                                                                                                       >  1:16291/1‑250 (MQ=255)                                                                                                                                aaGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATCTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATc                                                                                                                       >  2:541478/1‑250 (MQ=255)                                                                                                                                       tCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAACCTNgcgcn                                                                                                                >  2:784955/1‑249 (MQ=255)                                                                                                                                             aaTAGGT‑CACGAACGGTGCGAGGCTCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGATTAACGCCAATTGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAATACAAATACCTAATCCCCACAAGCGCGCGCCAaagca                                                                                                          >  2:932268/1‑248 (MQ=255)                                                                                                                                                     cACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCaaagatcggaag                                                                                                   >  1:134001/1‑240 (MQ=255)                                                                                                                                                            ggTGCGAGGATCGTGATGTATTTCGGGTTACTGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCaaaa                                                                                            >  1:974244/1‑250 (MQ=255)                                                                                                                                                                 gAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTatg                                                                                       >  1:604400/1‑250 (MQ=255)                                                                                                                                                                 gAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGACGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTatg                                                                                       >  1:657840/1‑250 (MQ=255)                                                                                                                                                                           aTGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTaagatcggaagag                                                                             >  1:652947/1‑242 (MQ=255)                                                                                                                                                                            tGTATTTCGGGTAACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCTCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCCGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATTATTAGTTCTAc                                                                            >  2:1028264/1‑250 (MQ=255)                                                                                                                                                                             gTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCACCAGCAAAACTATGATGAGTTCTACa                                                                           >  1:63191/1‑250 (MQ=255)                                                                                                                                                                                     ggTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCat                                                                   >  1:226725/1‑250 (MQ=255)                                                                                                                                                                                        tACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCCAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTANNANAtt                                                                >  2:957435/1‑250 (MQ=255)                                                                                                                                                                                             cgcAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGcct                                                           >  1:776817/1‑250 (MQ=255)                                                                                                                                                                                                 gcagcTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAGACTATGATGAGTTCTACAACTAGCATATTAGccttcct                                                       >  2:647707/1‑250 (MQ=255)                                                                                                                                                                                                                gatgaACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCTGGTTTAACGCCAATATCGTCAAGAACCAGCACCCCCAATCCGCCTGCAACACCAATACCTATTCCCCCCAATCTGGCGCCCAAGATGATCGCCAGCAAAATTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTATCAAGt                                        >  1:1287174/1‑250 (MQ=255)                                                                                                                                                                                                                         gATAGTCCAGACCGCCAGCAACCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAg                               >  1:581296/1‑250 (MQ=255)                                                                                                                                                                                                                                            aaCCTGCATGGCAG‑‑AAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTAc            >  1:711675/1‑250 (MQ=255)                                                                                                                                                                                                                                            aaCCTGCATGGCAG‑‑AAATAGAGGCGATAATCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACTAGCACCCCAAATCAGCCTGCAACACCAATACCTATTCCCCTCAATCTGGCGCCCAAGAAGACCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATAATCCTTCTTTGTTTTTTAACTAGTTGTTATTAGATTGTTATTTTTAAGTTAc            >  2:1296809/1‑250 (MQ=255)                                                                                                                                                                                                                                                        aaaaaaaaAGCGTAGAACACCGCCAAGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCAAAAg  <  2:776817/248‑1 (MQ=255)                                                                                                                                                                                                                                                          ||                                                                                                                                                                                                                                                       CGGTGATCGCAATCTGCGCGGATACCACTGCAGTAGACAGCGGACGGCAAGG‑TTTAACGCCTTGTTCCTTCGCAACTTCAGCGATAACTGGCAGTGTCGCCAGAGAGATGTT‑GCCAGTACCAGCAAAGATAGTCAGGAAATAGGT‑CACGATCGGTGCGAGGATCGTGATGTATTTCGGGTTACGGCGCAGCAGCTTTTCTGTCTGATGAACCAGATAGTCCAGACCGCCAGCAACCTGCATGGCAGAAAAATAGCGGCGATAACCGCCATGATAATGGAGATGACATCGAACGGGATGTTACCGGGTTTAACGCCAATAGCGGCAAGAACCAGCACCCCCAATCCGCCTGCAAAACCAATACCTATTCCCCCCAATCTGGCGCCCAAGAAGATCGCCAGCAAAACTATGATGAGTTCTACAACTAGCATATTAGCCTTCCTTGTTTTTTAACAAGTTGATATTAGATTGTTATTTTTAAGTTACTGCTCACAAG  >  NC_007779_cfb/4371025‑4371518

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