Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,983,589	C→A	14.8%	E180* (GAG→TAG)	hemC ←	hydroxymethylbilane synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,983,589	0	C	A	14.8%	71.5 / 3.6	27	E180* (GAG→TAG)	hemC	hydroxymethylbilane synthase
Reads supporting (aligned to +/- strand):  ref base C (12/11);  new base A (2/2);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-01

AGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGCAGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGTAACGCATCCAGACTGTCATAGTTATTGGACACAAAGGCATCGCGAGGATCTTCACGCTCACAAATAGTGAC  >  CP009273/3983355‑3983841                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                             aGATTTCGCCATCAATAAGCTCGGCGTAGCTAGCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGCAGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCTGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAatgatg                                                                                                                                                                                                                      <  1:117832/275‑1 (MQ=255)                                                                                                            tggtgATTCATCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGATTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGATCAGGCGAGTGCCGGCGTTGCCGCGCAGGGAGCGGATAATCAGATcc                                                                                                                                           >  1:514714/1‑243 (MQ=255)                                                                                                            tggtgATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACAGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGAc                                                                                                                                                                                                                                                         >  2:290875/1‑133 (MQ=255)                                                                                                                gATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTaa                                                                                                                                                                                                                                                     <  1:290879/133‑1 (MQ=255)                                                                                                                                                                         ccACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGTAAGGCATCCAGACTGTCATAGTTATTg                                              >  1:177086/1‑275 (MQ=255)                                                                                                                                                                           aCCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTCAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAg                                                                                                                                                                                      >  1:53740/1‑137 (MQ=255)                                                                                                                                                                           aCCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTAtt                                                                                                                                                                                                            >  2:58490/1‑115 (MQ=255)                                                                                                                                                                           aCCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCcg                                                                                                                                                                   >  2:296607/1‑156 (MQ=255)                                                                                                                                                                           aCCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGAGCGGATAATCAGATc                                                                                                                                            >  2:381497/1‑179 (MQ=255)                                                                                                                                                                           aCCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGTAACGCATCCAGACTGTCATAGTTATTGGa                                            >  1:295979/1‑275 (MQ=255)                                                                                                                                                                               cACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTcc                                                                                                                                                                                                                                                             <  1:97013/62‑1 (MQ=255)                                                                                                                                                                               cACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGcc                                                                                                                                                                                                        <  1:58490/115‑1 (MQ=255)                                                                                                                                                                               cACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGa                                                                                                                                                                                  <  2:53740/137‑1 (MQ=255)                                                                                                                                                                               cACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCTGCGCAGGGAGCGGATAATCAGATCCGGa                                                                                                                                        <  1:381505/179‑1 (MQ=255)                                                                                                                                                                               cACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCa                                                                                                                                                               <  1:296611/156‑1 (MQ=255)                                                                                                                                                                                   ttGTCCTACCGCCGGAAGAGCAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCGAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGAGCGGATAATCAGATCCGGAAGACGTTCAGCCAGTTGGCACTGGCAGCTTAAACAGGAAGACCCGACGACACCGACTGGCTGGTACGCACCCAGACAGTCCATCTGATTGGACACaaata                                    >  2:173712/1‑273 (MQ=255)                                                                                                                                                                                   ttGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTACAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTa                                                                                                                                                                                                              >  2:380465/1‑105 (MQ=255)                                                                                                                                                                                       ccTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTACAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCg                                                                                                                                                                                                          <  1:380473/105‑1 (MQ=255)                                                                                                                                                                                           ccgccgGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCg                                                                                                                                                                                                                                        >  1:329909/1‑71 (MQ=255)                                                                                                                                                                                              ccgGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTACAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTa                                                                                                                                                                                                    >  2:454774/1‑104 (MQ=255)                                                                                                                                                                                               cgGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTa                                                                                                                                                                                                                                    <  2:329904/71‑1 (MQ=255)                                                                                                                                                                                                  aaGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTACAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCa                                                                                                                                                                                                <  1:454783/104‑1 (MQ=255)                                                                                                                                                                                                        aagctcGGGTGGGAACGCGGCGCGAAGACGGGACTCCAGAACTTAACGTTTTAGTCCGGCTACGGCAAGAATGAGGGCATCGTATGCGCCGGTATCCAGTTTGCGCAGGCGAGTGCAGACGTGGCCGCGCAGGGAGCGGATAATCAGATCCGGACGGCGGTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGTAACGCATCCAGACTGTCATAGTTATTGGACACAAAGGCATCGCGAGGATCTTCACGCt               <  1:298926/272‑1 (MQ=255)                                                                                                                                                                                                            tcGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGCAACGCACCCCGCCTGCCCTCGGTCTTGGACAACAAAGCATCCCGCGGATCTTCACGCtga             >  2:238034/1‑271 (MQ=255)                                                                                                                                                                                                                      aCGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGGGCAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGTAACGCATCCAGACTGTCATAGTTATTGGACACAAAGGCATCGCGAGGATCTTCACGCGCACAAATAGTGAc  <  1:238037/274‑1 (MQ=255)                                                                                                                                                                                                                      aCGCGGCGCGAATACGGGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGTAACGCATCCAGACTGTCATAGTTATTGGACACAAAGGCATCGCGAGGATCTTCACGCTCACAAATAGTGAc  <  1:173714/274‑1 (MQ=255)                                                                                                                                                                                                                                   aCGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGa                                                                                                                                                           >  1:301414/1‑108 (MQ=255)                                                                                                                                                                                                                                       gACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGAGCgg                                                                                                                                                       <  2:301410/108‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                             AGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGCGCAGTGCAGTTTCGTGGTGATTCAGCGCGGCAAGCAGCTCGCGAGTGCGTGAATCATCAAGGCGGCATTCAATACCCACCGCACCTTGTCCTACCGCCGGAAGAGAAATCTCGGGTGGCAACGCGGCGCGAATACGTGACTCCAGACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGTAACGCATCCAGACTGTCATAGTTATTGGACACAAAGGCATCGCGAGGATCTTCACGCTCACAAATAGTGAC  >  CP009273/3983355‑3983841

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑