Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_007779_cfb	3,693,457	+A	noncoding (1417/1542 nt)	rrsC ←	16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_007779_cfb	3,693,457	1	.	A	97.4%	66.4 / NA	38	noncoding (1417/1542 nt)	rrsC	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base A (20/17);  minor base G (0/1);  total (20/18)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGTGTCCTGGGCCTCTAGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTC‑CCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGTGGCATTCTGATCCACGA  >  NC_007779_cfb/3693216‑3693523                                                                                                                                                                                                                                                   |                                                                   gCGGTGTCCTGGGCCTCTAGACGAAGGGGACGTATCAGTCGGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCGCTCACCATGgt                                                             <  1:787304/250‑1 (MQ=14)  cGGTGTCCTGGGCCTCTAGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGgtg                                                            >  2:205251/1‑250 (MQ=14)  cGGTGTCCTGGGCCTCTAGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGgtg                                                            <  1:357722/250‑1 (MQ=14)   ggTGTCCTGGGCCTCTAGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGgtgt                                                           >  1:1262016/1‑250 (MQ=14)      gtCCTGGGCCTCTAGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGNNGNGAc                                                        >  2:971107/1‑250 (MQ=14)       tCCTGGGCCTCTAGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACg                                                       <  1:1079677/250‑1 (MQ=14)       tCCTGGGCCTCTAGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACg                                                       <  1:1269968/250‑1 (MQ=14)           gggCCTCTAGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCg                                                   >  1:1067074/1‑250 (MQ=14)             gCCTCTAGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCGGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTc                                                                     >  1:975326/1‑230 (MQ=14)             gCCTCTAGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCGGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTGTGCAACCCACTc                                                                     <  2:975326/230‑1 (MQ=14)              ccTCTAGACGAAGGGGACGAATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGgtg                                                <  2:680006/250‑1 (MQ=14)                tctAGACGCAGGGGCCGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGgtgtg                                              >  2:691746/1‑250 (MQ=14)                 ctAGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGgtgtgt                                             >  1:571253/1‑250 (MQ=14)                  tAGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTa                                            <  1:848953/250‑1 (MQ=14)                   aGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGgtgggtca                                           >  2:1082930/1‑248 (MQ=14)                   aGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTAc                                           >  2:427970/1‑250 (MQ=14)                    gACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACa                                          <  2:144899/250‑1 (MQ=14)                     aCGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACaa                                         >  1:1351020/1‑250 (MQ=14)                      cGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAg                                        >  2:1265376/1‑250 (MQ=14)                       gAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAgg                                       <  1:1265376/250‑1 (MQ=14)                         aGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGcc                                     >  1:244284/1‑250 (MQ=14)                          ggggACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGccc                                    <  1:413282/250‑1 (MQ=14)                          ggggACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGccc                                    <  1:691746/250‑1 (MQ=14)                           gggACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTcacca                                                                 >  1:1006049/1‑216 (MQ=14)                           gggACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTcacca                                                                 <  2:1006049/220‑5 (MQ=14)                                gTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAAc                              <  1:1082930/250‑1 (MQ=14)                                  aTCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGt                            >  1:108616/1‑250 (MQ=14)                                    cAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTAt                          >  1:440496/1‑250 (MQ=14)                                      ggcNGNNTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTc                        <  2:950931/248‑1 (MQ=14)                                       tCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCa                       <  2:405838/250‑1 (MQ=14)                                        cTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAGCGTATTCAc                      <  2:1067074/250‑1 (MQ=14)                                            ttCGCAAGACGCCTTGCCTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCCTGTTACGACTTCACCCCAGTCATGAATCGCAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCGCCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGTg                  <  2:1274493/250‑1 (MQ=9)                                              cGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGTGGc                >  1:785809/1‑250 (MQ=14)                                                cAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGTGGCAt              <  2:244284/250‑1 (MQ=14)                                                 aaGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGTGGCAtt             >  1:383528/1‑250 (MQ=14)                                                    aCGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGTGGCATTCTg          >  2:859673/1‑250 (MQ=14)                                                      gCCTCGCTTTTCACTTCCCATCAGACAATCTTTGTGAGCACTGCATAGAACGCTTCTTTAAGGCAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAATGTGGTAAGTGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGACGTGGAACGCTTTCACCGTGGTATNCTGan        >  2:792665/1‑249 (MQ=11)                                                       ccTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGCTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGTGGCCTTCTGATc       >  1:1047979/1‑250 (MQ=14)                                                        tttGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGTGGCATTCTGAGcc      <  1:1076585/249‑1 (MQ=14)                                                         ttGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGTCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGAGGCATTCTGATcaa     >  2:1076585/1‑248 (MQ=9)                                                          tGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGTGGCATTCTGATCCAc    >  1:256458/1‑250 (MQ=14)                                                           gCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGTTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGTGGCATTCTGATCCACg   >  2:171015/1‑250 (MQ=14)                                                            cTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGTGGCATTCTGATCCACGa  <  1:971107/250‑1 (MQ=14)                                                                                                                                                                                                                                                   |                                                                   GCGGTGTCCTGGGCCTCTAGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTC‑CCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGTGGCATTCTGATCCACGA  >  NC_007779_cfb/3693216‑3693523

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