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	93.1%	45.0 / ‑4.6	29	noncoding (1417/1542 nt)	rrsC	16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base A (17/10);  minor base T (0/2);  total (17/12)
Fisher's exact test for biased strand distribution p-value = 1.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTC‑CCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGTGGCATTCTGATCCACGATTACTAGCG  >  NC_007779_cfb/3693209‑3693532                                                                                                                                                                                                                                                          |                                                                            tGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTca                                                                             <  2:36289/250‑2 (MQ=14)        gCGGTGTCCTGGGCCTCTAGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGgt                                                                      <  2:549230/250‑1 (MQ=14)               ccTGGGCCTCTAGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCTCCATGGTGTGACgg                                                               <  2:480714/250‑1 (MQ=14)                      ctctAGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGgtgt                                                        <  2:1231326/250‑1 (MQ=14)                        ctAGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGgtgtgt                                                      >  1:971880/1‑250 (MQ=14)                         tAGACGAAGGGGCCGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTCTGCAACCCACTCACCATGgt                                                                      >  2:1079065/1‑233 (MQ=9)                         tAGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTCTGCAACCCACTCACCATGgt                                                                      <  1:1079065/233‑1 (MQ=14)                                 ggggACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGCAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGccc                                             >  1:749691/1‑250 (MQ=14)                                    gACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCGACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCggg                                          >  2:453488/1‑250 (MQ=14)                                        tATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACg                                      <  1:509590/250‑1 (MQ=14)                                           cAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGtgt                                   <  1:431403/250‑3 (MQ=14)                                           cAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTAt                                   >  1:1436095/1‑250 (MQ=14)                                              tCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCa                                >  2:686310/1‑250 (MQ=14)                                                 gCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGgt                                                          <  1:1269734/221‑1 (MQ=14)                                                 gCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGgt                                                          >  2:1269734/1‑221 (MQ=14)                                                 gCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCg                             >  1:65333/1‑250 (MQ=14)                                                  cTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGAGTGAGCACTGCAAAGAACGCCTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGTCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAATGCCCGGGAACGTATTccctgt                            >  2:662201/1‑245 (MQ=255)                                                   ttCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTATCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGAANTCAACGTg                           >  2:796560/1‑250 (MQ=9)                                                    tCGCAAGACGCCTTGCCTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAAGCCACTCACCATGGTGTGACGGGc                                                             <  2:999050/215‑1 (MQ=9)                                                    tCGCAAGACGCCTTGCCTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGc                                                             >  1:999050/1‑215 (MQ=14)                                                      gCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGcc                                              >  2:1113409/1‑228 (MQ=14)                                                      gCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGcc                                              <  1:1113409/228‑1 (MQ=14)                                                        aaGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGTGGCAtt                      <  1:686310/250‑1 (MQ=14)                                                        aaGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGTGGCAtt                      >  1:222086/1‑250 (MQ=14)                                                          gACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGTGGCATTCt                    >  1:465427/1‑250 (MQ=14)                                                                  gCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCTCCATGGTGTGACGGGCGGTGTGTACAAGGcc                                              <  1:583307/216‑1 (MQ=14)                                                                  gCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACACGGcc                                              >  2:583307/1‑216 (MQ=14)                                                                   cTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGTGGCATTCTGATCCACGa           >  1:465366/1‑250 (MQ=14)                                                                           ttCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGAGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGTGGCATTCTGATCCACGATTACTAGc   >  1:1136663/1‑250 (MQ=14)                                                                            tCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTCACCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCGCCGTGGCATTCTGATCCACGATTACTTGCg  <  2:1436095/250‑1 (MQ=11)                                                                             ccTATCAGACAATCTGTGTGAGCGCTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTcacc                                                                           <  1:687786/176‑4 (MQ=11)                                                                             ccTATCAGACAATCTGCGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTcacc                                                                           >  2:687786/1‑173 (MQ=11)                                                                                                                                                                                                                                                          |                                                                            TGAAAGGGCGGTGTCCTGGGCCTCTAGACGAAGGGGACGTATCAGTCTGCTTCGCAAGACGCCTTGCTTTTCACTTCCTATCAGACAATCTGTGTGAGCACTGCAAAGAACGCTTCTTTAAGGTAAGGAGGTGATCCAACCGCAGGTTCCCCTACGGTTACCTTGTTACGACTTCACCCCAGTCATGAATCACAAAGTGGTAAGCGCCCTCCCGAAGGTTAAGCTACCTACTTCTTTTGCAACCCACTC‑CCATGGTGTGACGGGCGGTGTGTACAAGGCCCGGGAACGTATTCACCGTGGCATTCTGATCCACGATTACTAGCG  >  NC_007779_cfb/3693209‑3693532

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