Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,176,209	C→A	D209Y (GAT→TAT)	cpdA ←	cAMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,176,209	0	C	A	91.7%	65.3 / NA	24	D209Y (GAT→TAT)	cpdA	cAMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (7/15);  total (9/15)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGA  >  NC_000913/3176070‑3176339                                                                                                                                            |                                                                                                                                   atCAGCGTGCTATCAGCATGTAGCTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTGGGAACAGTGCGGCTTAAACGGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAg                                                                                                                                 <  1:484716/142‑1 (MQ=255)   cAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGGCCATGCGCGGCGAGATCGAGcc                                                                                                                               >  2:479371/1‑142 (MQ=255)         gccaggaGCATGTCACCCGCGTGTACCCCCGCCGGGCCCGCTGGTAGCCAGCGTAAGGTTGGAACAGTCCCGCTTAAACCGCACCCAGCTCGCCGGCGCTGCCCGCAGGCGGCGACCATTCCAGTCGAGCTACAGCTCCtgac                                                                                                                         <  2:404780/137‑2 (MQ=255)             tccgcgtGTAACTCGAGCGCACACCAGTCGGGCGCGATCGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTACACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAat                                                                                                                     <  2:492526/137‑1 (MQ=255)                     tAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAACAAGACCGc                                                                                                             >  2:161572/1‑143 (MQ=255)                       aCTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGcac                                                                                                           <  1:242519/143‑1 (MQ=255)                         tCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAAGATGACCGAACAg                                                                                                         >  1:115144/1‑143 (MQ=255)                         tCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAg                                                                                                         >  1:534420/1‑143 (MQ=255)                            cgcggaCGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCaa                                                                                                      <  1:371956/138‑1 (MQ=255)                                 aCGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACTGTTCGACGGCGTTGTTAGCAGGCGGCGACCATTCCAGTCGAGAGACAGCTCCTGATGAATATGACCGCACAGCAAGTAtt                                                                                                 <  1:294476/143‑1 (MQ=255)                                  cGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTAttt                                                                                                <  1:378966/143‑1 (MQ=255)                                              gcgATCGTATCCAGCGGCAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGaaa                                                                                    <  1:179506/143‑1 (MQ=255)                                               cgATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACAAAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAAc                                                                                   <  1:523645/143‑1 (MQ=255)                                                    gTATCCAGCGGAACGTGGCAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGc                                                                              <  2:23582/143‑1 (MQ=255)                                                     tATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCGCAGCAAGTATTTGACGTGCGGAAACTTCGcc                                                                             >  1:107468/1‑143 (MQ=255)                                                              tAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCACCAGGCGGCGACCATTACAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTACTTGACGTGGGGAAACTTCGACAACACgaaa                                                                    >  2:405317/1‑140 (MQ=255)                                                                      gAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTc                                                            >  1:323450/1‑143 (MQ=255)                                                                      gAACAGTGCGGCTTAAACTGCACACAGGCCGCCGGCGTTGCCAGCAGGCGGCGACCACTCCATTCGAGCTACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTc                                                            <  1:294166/143‑1 (MQ=255)                                                                       agcggTGCGGCTTAAACGGCAGGCAGGTTGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGAGACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCg                                                           <  2:47019/139‑1 (MQ=255)                                                                                gCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTa                                                  <  2:244301/143‑1 (MQ=255)                                                                                    aaaCTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCa                                              >  1:496677/1‑143 (MQ=255)                                                                                      accgCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAACATGACCGCACAGCAAGTACTACACGGTCTGAAACGTCGCCTGCCCGGTATCCAGTTCGCCCGCGTTACGCAGa                                            <  1:528176/140‑1 (MQ=255)                                                                                                    gACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTAGTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTGCGCCCCCGGTACGCAGACTGAGTCGGTCGAg                              >  2:32883/1‑143 (MQ=255)                                                                                                                                tCCAGTCGAGATACAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGgatga  <  1:572966/143‑1 (MQ=255)                                                                                                                                            |                                                                                                                                   GTCAGCGTGCCATCAGCATGTAACTCGAGAGTACGCCAGCCGGGCGCGATGGTATCCAGCGTAAAGTTGGAACAGTGCGGCTTAAACTGCACACAGGTCGACGGCGTTGCCAGCAGGCGGCGACCATTCCAGTCGAGATCCAGCTCCTGATGAATATGACCGCACAGCAAGTATTTGACGTGCGGAAACTTCGCCAGCACGGTATCCAGTTCGCCCGCGTTACGCAGACTGTGTTGATCGAGCCAACTACAACCCGCAGGTAGCGGATGA  >  NC_000913/3176070‑3176339

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