Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	3,544,630	G→A	11.3%	H401H (CAC→CAT)	malP ←	maltodextrin phosphorylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	3,544,630	0	G	A	11.3%	91.0 / 4.4	37	H401H (CAC→CAT)	malP	maltodextrin phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (16/16);  new base A (2/2);  total (18/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 = 7.12e-01

AGCCGCCAGTGCCGGGTTGCACTGTTTGATCCAGCGACGTGGGGTAATACCGTTGGTGACGTTATGGAATTTGTTCGGCCATAGCTGGTGATATTCCGGGAACAGATCTTTCACCACCAGATCCGAGTGCAGCGCCGCAACACCGTTCACCGCGAAACCGCCAACCACACACAGGTTCGCCATATGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGCCCACACTTTTTCATCGCCCGGCCAGGTTTTCTCTACCAGCGTTTTAAAGCGAGTATTAATTTCGTTAATAATCTGCATGTGGCGCGGCAGTAAGCCTTTCACCAGTTTCACATCCCAGCGTTCCAGCGCTTCTGGCATCAGGGTATGGTTGGTGTAAGCGAAAGTTTTGCTGGTAATGGCCCAGGCGTCATCCCAGCTCATCTGGTGCTCATCGATCAGCACGCGCAGCAGTTCTGGAATCGCGATAGT  >  CP009273/3544432‑3544899                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                              agCCGCCAGTGCCGGGTTGCACGGTTTGATCCAGCGACGTGGGGTAATACCGTTGGTGACGTTATGGAATTTGTTCGGCCATAGCTGGTGATATTCCGGGAACAGATCTTTCACCACCAGATCCGAGTGCAGCGCCGCAACACCGTTCACCGCGAAACCGCCAACCACACACAGGTTCGCCATATGCCCTTGTTTGTCGTGCACTACCGCCAGTTTGTCCCCCACTTCTTAATCGCCCGGCCAGGTTGTCTCTACCCGCGTTTTAAAGCtattat                                                                                                                                                                                                   >  2:634713/1‑275 (MQ=255)                                    gACGTGGGGTAATACCGTTGTTGACTTTATGGAATTTGTTCGGCCATAGCTGGTGATATTCCGGGAACAGATCTTTCACCACCAGATCCGAGTGCAGCGCCGCAACACCGTTCACCGCGAAAGCGCCAACCACACACAGGTTCGCCATAGGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGCCCACACTTTTTCATCGCCCGGCCAGGTTTTCTCTACCAGCGTTTTAAATCGAGTATTAATTTCGTTAATAATCTGCATGTGGCGCGGCAGt                                                                                                                                                                <  1:634724/275‑1 (MQ=255)                                                     ttGGTGACGTTATGGAATTTGTTCGGCCATAGCTGGTGATATTCCGGGAACAGATCTTTCACCACCAGATCCGAGTGCAGCGCCGCAACACCGTTCACCGCGAAACCGCCAACCACACACAGGTTCGCCATATGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGCCca                                                                                                                                                                                                                                                        >  1:346723/1‑170 (MQ=255)                                                         tGACGTTATGGAATTTGTTCGGCCATAGCTGGTGATATTCCGGGAACAGATCTTTCACCACCAGATCCGAGTGCAGCGCCGCAACACCGTTCACCGCGAAACCGCCAACCACACACAGGTTCGCCATATGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGCCCACACt                                                                                                                                                                                                                                                    <  2:346718/170‑1 (MQ=255)                                                          gACGTTATGGAATTTGTTCGGCCATAGCTGGTGATATTCCGGGAACAGATCTTTCACCACCAGATCCGAGTGCAGCGCCGCAACACCGTTCACCGCGAAACCGCCAACCACACACAGGTTCGCCATATGCACTTGTTTGTCATGCACCACCGCCAGTTTGGCACACACTTTTTCATCGCCCGGCCAGGTTTTCCCTACCAGCGTTTTAAAGCGAGTATTAATCTCGTTAATAATCTTCATGTGGCGCTNTATGAATCCTTCCACCATTCTcaca                                                                                                                                           >  2:443690/1‑274 (MQ=255)                                                                              gcccttgtcgttcatatTCCGTGACCAGATCTTTCACCACCAGATCCGCGTGCATCGCCGCAACTCCGTTCACCGCGAAACCGCCAACCACTCACAGGTTCGCCATATGCACTTGTTTGTCCTGCACCACCTCCAGTTTGGCCCACACTTTTTCATCGCCCGGCCAGGTTTTCTCTACCATCGTTTTAAAGCGAGTATTAATTTCGTTAATAATCTGCATGTGGCGCGGCAGTAAGCCTTTCACCAGTTTCACATCCCAGCGTTCCAGCGCTTCt                                                                                                                      <  1:168461/262‑1 (MQ=255)                                                                                                                    accaGATCCGAGTGCAGCGCCGCAACACCGTTCACCGCGAAACCGCCAACCACACACAGGTTCGCCATATGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGCCCACACTTTTTCATCGCCCGGCCAGGTTTtct                                                                                                                                                                                                                           >  2:361905/1‑136 (MQ=255)                                                                                                                        gATCCGAGTGCAGCGCCGCAACACCGTTCACCGCGAAACCGCCAACCACACACAGGTTCGCCATATGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGCCCACACTTTTTCATCGCCCGGCCAGGTTTTCTCTAc                                                                                                                                                                                                                       <  1:361910/136‑1 (MQ=255)                                                                                                                         ctCCGAGTGCAGCGCCGCAACACCGTTCACCGCGAAACCGCCAACCACCCACAGGTTCGCCATATGCACTTGTTTGTCATGCNCCACCGCCNTTTTGGCCCACACTTTTTCATCGCCCGGCCAGGTTTTCTCTACCAGCGTTTTAAAGCGAGTATTAATTTCGTTAATAATCTGCATGTGGCGCGGCAGTAAGCCTTTCACCAGTTTCACATCCCAGCGTTCCAGCGCTTCTGGCATCAGTGTATGGTTGGTGTAAGCGAAAGTTTTGCTGGTaa                                                                           <  1:443696/274‑1 (MQ=255)                                                                                                                              aGTGCAGCGCCGCAACACCGTTCACCGCGAAACCGCCAACCACACACAGGTTCGCCATATGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGCCCACACTTTTTCATCg                                                                                                                                                                                                                                           >  2:377528/1‑110 (MQ=255)                                                                                                                                  cAGCGCCGCAACACCGTTCACCGCGAAACCGCCAACCACACACAGGTTCGCCATATGCACTTGTTTGTCGTGCACCACCGCCAGttt                                                                                                                                                                                                                                                              >  1:309600/1‑87 (MQ=255)                                                                                                                                  cAGCGCCGCAACACCGTTCACCGCGAAACCGCCAACCACACACAGGTTCGCCATATGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGCCCACACTTTTTCATCGCCCg                                                                                                                                                                                                                                       <  1:377533/110‑1 (MQ=255)                                                                                                                                  cAGCGCCGCAACACCGTTCACCGCGAAACCGCCAACCACACACAGGTTCGCCATATGCACTTGTTTGTCATGCACCACCGCCAGTTTGGCCCACACTTTTTCATCGcc                                                                                                                                                                                                                                         >  2:467227/1‑108 (MQ=255)                                                                                                                                     cgccgcAACACCGTTCACCGCGAAACCGCCAACCACACACAGGTTCGCCATATGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGCCCACACTTTTTCATCGCCCGGCCAGGTTTTCTCTACCAGCGTTTTAAAGCGAGTATTAATTTCGTTAATAATCTGCATGTGGCGCGGCAGTAAGCCttt                                                                                                                                                        >  2:259110/1‑186 (MQ=255)                                                                                                                                      gccgcAACACCGTTCACCGCGAAACCGCCAACCACACACAGGTTCGCCATATGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGcc                                                                                                                                                                                                                                                          <  2:309595/87‑1 (MQ=255)                                                                                                                                      gccgcAACACCGTTCACCGCGAAACCGCCAACCACACACAGGTTCGCCATATGCACTTGTTTGTCATGCACCACCGCCGGTTTTTCCCACACTTTTTCATCGCCCGGc                                                                                                                                                                                                                                     <  1:467234/108‑1 (MQ=255)                                                                                                                                         gcAACACCGTTCACCGCGAAACCGCCAACCACACACAGGTTCGCCATATGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGCCCACACTTTTTCATCGCCCGGCCAGGTTTTCTCTACCAGCGTTTTAAAGCGAGTATTAATTTCGTTAATAATCTGCATGTGGCGCGGCAGTAAGCCTTTCAcc                                                                                                                                                    <  1:259114/186‑1 (MQ=255)                                                                                                                                              acCGTTCACCGCGAAACCGCCAACCACACACAGGTTCGCCATATGCACTTGTGTGTCGTGCACCAGCGCCAGTGTGGCCCACACTTTTTCATCGCCCGGCCAGGtttt                                                                                                                                                                                                                             <  1:353467/108‑1 (MQ=255)                                                                                                                                                  ttCACCGCGAAACCGCCAACCACACACAGGTTCGCCATATTCACTTGTTTGTCGTGCACCACCGCCAGTTTTGCCCACCCTTTTTCATCGCCCGGCCAGGTTTTCTCTACCAGCGTTTTAAAGCGAGTATTAATTTCGt                                                                                                                                                                                          >  1:578069/1‑139 (MQ=255)                                                                                                                                                  ttCACCGCGAAACCGCCAACCACACACAGGTTCGCCATATGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGCCCACACTTTTTCATCGCCCGGCCAGGTTTTCTCTACCAGCGTTTTAAAGCGAGTATTAAtt                                                                                                                                                                                              >  2:581047/1‑135 (MQ=255)                                                                                                                                                  ttCACCGCGAAACCGCCAACCACACACAGGTTCGCCATATGCACTTGTTTGTCGTGCACCAACGCCAGTTTGGCCCACACTTTTTCATCGCCCGGCCAGGTTTTCTCTACCAGCGTTTTAAAGCGAGTATTAATTTCGTTAATAATCTGCATGTGGCGCGGCAGTAAGCCATTCACCAGTTTCACATCCCAGCGTTCCAGCGCTTCTGGCATCAGGGTATGGTTGGTGTAAGCGAAAGTTTTGCTGGTAAt                                                                          >  1:418746/1‑251 (MQ=255)                                                                                                                                                    cACCGCGAAACAGCCAACCACACACAGGTTCGCCATATGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGCCCACACTTTTTCATCGCCCGGCCAGGTTTTctct                                                                                                                                                                                                                         >  2:528373/1‑106 (MQ=255)                                                                                                                                                      ccGCGAAACCGCCAACCACACACAGGTTCTCCATATGCACTTTTTTGTCGTGCACCACCGCCAGTTTGGCCCACACTTTTTCATCGCCCGGCCAGGTTTTCTCTACCAGCGTTTTAAAGCGAGTATTAATTTCGTtaat                                                                                                                                                                                      <  2:578058/139‑1 (MQ=255)                                                                                                                                                      ccGCGAAACCGCCAACCACACACAGGTTCGCCATATGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGCCCACACTTTTTCATCGCCCGGCCAGGTTTTCTCTACCAGCGTTTTAAAGCGAGTATTAATTTCGt                                                                                                                                                                                          <  1:581058/135‑1 (MQ=255)                                                                                                                                                        gcgAAACCGCCAACCACACACAGGTTCGCCATATGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGCCCACACTTTTTCATCGCCCGGCCAGGTTTTCTCTACCa                                                                                                                                                                                                                     <  1:528384/106‑1 (MQ=255)                                                                                                                                                                               gTTCGCCATATGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGCCCACACTTTTTCa                                                                                                                                                                                                                                              >  2:423082/1‑58 (MQ=255)                                                                                                                                                                               gTTCGCCATATGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGCCCACACTTTTTCATCGCCCGGCCAGGTTTTCTCTACCAGCGTTTTAAAGCGAGTATTAATTTCGTTAATAATCTGCAtgt                                                                                                                                                                           >  1:238104/1‑125 (MQ=255)                                                                                                                                                                                   gCCATATGCACTTGTTTGTCGTGCACCACCGCCAGTTTgg                                                                                                                                                                                                                                                            >  2:219604/1‑40 (MQ=255)                                                                                                                                                                                   gCCATATGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGCCCACACTTTTTCATCGc                                                                                                                                                                                                                                          <  1:423088/58‑1 (MQ=255)                                                                                                                                                                                   gCCATATGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGCCCACACTTTTTCATCGCCCGGCCAGGTTTTCTCTACCAGCGTTTTAAAGCGAGTATTAATTTCGTTAATAATCTGCATGTGgcg                                                                                                                                                                       <  2:238100/125‑1 (MQ=255)                                                                                                                                                                                       tatGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGCCca                                                                                                                                                                                                                                                        <  1:219608/40‑1 (MQ=255)                                                                                                                                                                                       tatGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGCCCACACTTTTTCATCGCCCGGCCAGGTTTTCTCTACCAGCGTTTTAAAGCGAGTATTAATTTCGTTAATAATCTGCATGTGGCGCGGGAGTAAGCCTTTCACCAGTTTCACATCCCAGCGTTCCAGCGCTCCTCGCATCCGGCAATGGTTGGTGGAAGCGTAAGTTTCGCTGGTAATGGCCCAGGAGTAATCCAAGCTCATCGGACGCACCTCGATCAGCACGCgcagccgctccgcc             >  2:624344/1‑266 (MQ=255)                                                                                                                                                                                             cTTGTTTGTCGTGCACCACCGCCAGTTTGGCCCACACtt                                                                                                                                                                                                                                                   >  1:126934/1‑39 (MQ=255)                                                                                                                                                                                                 tttgtCGTGCACCACCGCCAGTTTGGCCCACACTTTTTc                                                                                                                                                                                                                                               <  2:126933/39‑1 (MQ=255)                                                                                                                                                                                                 ttggtCGTGCCCCTCCGCCAGTTTGGCCCACACGTTTTCATCGCCCGGCCAGGTTTGCTCTACCAGCGTTTTAAAGCGAGTATTAATTTCGTTAATAATCTGCATGTGGCGCGGCAGTAAGCGTTTCACCAGTTTCACATCCCAGCGTTCCAGCGCTTCTGGCATCAGGGTATGGTTGGTGTAAGCGAAAGTTTTGCTGGTAATGTCCCAGGCGTCATCCCAGCTCATCTGGTGCTCATCGATCAGCACGCGCAGCAGTTCTGGAATCGCGATAg   <  1:624355/272‑1 (MQ=255)                                                                                                                                                                                                  ttgtCGTGCACCACCGCCAGTTTGGCCTACACTTTTTCATCGCCCGGCCAGGTTTTCTCTACCAGCGTTTTAAAGCGAGTATTAATTTCGTTAATAATCTGCATGTGGCGCGGCAGTAAGCCTTTCACCAGTTTCACATCCCAGCGTTCCAGCGCTTCTGGCATCAGGGTATGGTTGGTGTAAGCGAAAGTTTTGCTGGTAATGGCCCAGGCGTCATCCCAGCTCATCTGGTGCTCATCGATCAGCACGCGCAGCAGTTCTGGAATCGCGATAGt  <  1:6210/275‑1 (MQ=255)                                                                                                                                                                                                      cGTGCACCACCGCCAGTTTGGCCCACACTTTTTCATCGCCCGGCCAGGTTTTctct                                                                                                                                                                                                                         >  2:99393/1‑56 (MQ=255)                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                              AGCCGCCAGTGCCGGGTTGCACTGTTTGATCCAGCGACGTGGGGTAATACCGTTGGTGACGTTATGGAATTTGTTCGGCCATAGCTGGTGATATTCCGGGAACAGATCTTTCACCACCAGATCCGAGTGCAGCGCCGCAACACCGTTCACCGCGAAACCGCCAACCACACACAGGTTCGCCATATGCACTTGTTTGTCGTGCACCACCGCCAGTTTGGCCCACACTTTTTCATCGCCCGGCCAGGTTTTCTCTACCAGCGTTTTAAAGCGAGTATTAATTTCGTTAATAATCTGCATGTGGCGCGGCAGTAAGCCTTTCACCAGTTTCACATCCCAGCGTTCCAGCGCTTCTGGCATCAGGGTATGGTTGGTGTAAGCGAAAGTTTTGCTGGTAATGGCCCAGGCGTCATCCCAGCTCATCTGGTGCTCATCGATCAGCACGCGCAGCAGTTCTGGAATCGCGATAGT  >  CP009273/3544432‑3544899

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