Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,572,964	C→T	12.1%	A747T (GCA→ACA)	pqqL ←	putative zinc peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,572,964	0	C	T	12.1%	61.6 / 2.1	33	A747T (GCA→ACA)	pqqL	putative zinc peptidase
Reads supporting (aligned to +/- strand):  ref base C (11/17);  new base T (2/2);  total (14/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

ACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCG  >  NC_000913/1572824‑1573101                                                                                                                                             |                                                                                                                                          tctctgATCGTTGGCCAGTGCGACGTTAAAAGGATCGAGCGCCATACGCGTCGGCATATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGATACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGt                                                                                                                                         <  2:293978/138‑1 (MQ=255)    tAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTATCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGATTATTTTGTACTTTTAGAGCACCACATGAGTTGt                                                                                                                                      >  2:457779/1‑143 (MQ=255)     aGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGCAACCTGTGCCACCGGTTCATTTTGTTCTCTTACAGTAACCGATGCGTTGTc                                                                                                                                     >  1:374930/1‑143 (MQ=255)       aTCGTTTGCCAGTGGGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCg                                                                                                                                   <  2:255524/143‑1 (MQ=255)                   tGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTa                                                                                                                       >  2:101382/1‑143 (MQ=255)                      gACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTACCCGATGCGTAGTCCGTCGCGCGAGTTAATg                                                                                                                    >  2:84808/1‑143 (MQ=255)                      gACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATg                                                                                                                    <  1:341448/143‑1 (MQ=255)                            aaaaGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACGGGTCCGGTATTGGCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTAc                                                                                                              >  1:247923/1‑143 (MQ=255)                                 cATCGAGCGGCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCg                                                                                                         <  1:405272/143‑1 (MQ=255)                                      aGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTaa                                                                                                    <  1:500779/143‑1 (MQ=255)                                           cATACGCGTCGGCAGATTAATAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCg                                                                                               <  1:60292/143‑1 (MQ=255)                                             tACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGaa                                                                                             <  2:289694/143‑1 (MQ=255)                                                 cgTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCag                                                                                         <  1:349929/143‑1 (MQ=255)                                                        aGATTAACAGGTGTCCGGGAATCATAACGCTTCCATTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGttt                                                                                  <  2:479349/143‑1 (MQ=255)                                                                    gtCCGGGAATTATAACGCTTCCTCTGTGAAACCTGGGCCACTGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTaa                                                                      <  2:298787/143‑1 (MQ=255)                                                                                gagcGCTGCCACGGTGAAATCGCTGCCACAGGTTCATGTTGTTCTTTTACAGTAACCGATGCGTTGGCCGTCGCGCGGGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAAt                                                          <  2:163640/140‑1 (MQ=255)                                                                                 aaCGCTTCCACTGTGAAACCTGTGCTACAGTTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGGCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAAtt                                                         <  1:29441/143‑1 (MQ=255)                                                                                    gCTCCCTCTGTTAGACGTGTGCGAAACTTTCATTTCGTTCTTTTATATTAATCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAAc                                                      <  2:362432/143‑1 (MQ=255)                                                                                      ttCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGc                                                    >  1:208621/1‑143 (MQ=255)                                                                                       tCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTGTGACTGATCCTAAGTAACGCGTAATTAACGcc                                                   >  2:278889/1‑143 (MQ=255)                                                                                                 caCCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTTTCCGTCGCGCGAGTTAATCGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGt                                         <  1:30935/142‑1 (MQ=255)                                                                                                  aCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGTGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGtc                                        >  1:226739/1‑143 (MQ=255)                                                                                                  aCCTGTGGCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGTGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGtc                                        >  1:226639/1‑143 (MQ=255)                                                                                                             aGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACAt                             <  1:2124/143‑1 (MQ=255)                                                                                                                    ttttGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGAc                      >  1:339750/1‑143 (MQ=255)                                                                                                                      ttGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCCAcca                    >  2:136336/1‑141 (MQ=255)                                                                                                                       tGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAAt                   >  2:367449/1‑143 (MQ=255)                                                                                                                            ttttACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACaa              >  1:524469/1‑143 (MQ=255)                                                                                                                                       ccGATGCGTTGTCCGTCGCGCGAGTTCATGGTTTACCTGCGGCTAATGGCGACTAAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATcc   >  1:294020/1‑143 (MQ=255)                                                                                                                                        cGATGTGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCg  <  2:226639/143‑1 (MQ=255)                                                                                                                                        cGATGTGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCg  <  2:226739/143‑1 (MQ=255)                                                                                                                                        cGATGCGTTTTCCGTCGCGCGAGTTAATGGTTTACCTGGGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCg  <  2:11384/143‑1 (MQ=255)                                                                                                                                        cGATGCGTTGTCCGTCGCGTGAGTTAATGGTTTACTTGCGGCTGATGGCGAATCAGAGTGTTTGATTGATCCTATGTAACGCGTAATTAACGCCACGAGTTTTTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCg  <  2:165233/143‑1 (MQ=255)                                                                                                                                             |                                                                                                                                          ACGTAGATCTTTTGCCAGTGCGACGTTAAAAGCATCGAGCGCCATACGCGTCGGCAGATTAACAGGTGTCCGGGAATCATAACGCTTCCACTGTGAAACCTGTGCCACAGGTTCATTTTGTTCTTTTACAGTAACCGATGCGTTGTCCGTCGCGCGAGTTAATGGTTTACCTGCGGCTAATGGCGAATCAGAGTGTTTGATTGATCCTAAGTAACGCGTAATTAACGCCACGAGTTTGTCTTCTGCGACATTACCGACAATGACAAACGTGATATCCG  >  NC_000913/1572824‑1573101

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑