Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	1,820,438	G→T	17.3%	P248Q (CCG→CAG)	astE ←	succinylglutamate desuccinylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	1,820,438	0	G	T	17.3%	40.1 / 2.1	23	P248Q (CCG→CAG)	astE	succinylglutamate desuccinylase
Reads supporting (aligned to +/- strand):  ref base G (13/6);  new base T (2/2);  total (15/8)
Fisher's exact test for biased strand distribution p-value = 5.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

GCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCG  >  CP009273/1820187‑1820707                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                              gCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACTGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACAGAGGGGCGGGGTTCTAACGATACCGac                                                                                                                                                                                                                                                         >  1:227872/1‑274 (MQ=255)                                  ttttttCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGGGTTCGCACGATACCGACACTCTCACCACATCGCAGCGCCGCCATTgcgtg                                                                                                                                                                                                                        >  2:219522/1‑272 (MQ=255)                                                                                   aGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACAGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCCGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGCGCAACGCtt                                                                                                                                                                      >  1:207445/1‑275 (MQ=255)                                                                                         tAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCCCACCAGATAGCAGGGCAGACATTGCGCGGGCAGTTACTTCAAACTTGCGAAGAACGTTTTGCCCAAAAGGCAACCGTTTACCaa                                                                                                                                                                >  1:45790/1‑275 (MQ=255)                                                                                                            tGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCTGCGGTGTTCTCACGATACCGACACTCTCACCAGATagca                                                                                                                                                                                                                                        >  1:66217/1‑184 (MQ=255)                                                                                                               gTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCGTTGCCAAGTTCAAGAGAACAGGCCaacac                                                                                                                                          >  1:136441/1‑271 (MQ=255)                                                                                                                ttACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCTGCGGTGTTCTCACGATACCGACACTCTCACCAGATAgcagcgc                                                                                                                                                                                                                                    <  2:66216/184‑1 (MQ=255)                                                                                                                                       tCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAg                                                                                                                                                                                                                                              >  1:395117/1‑151 (MQ=255)                                                                                                                                           gcgcCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATag                                                                                                                                                                                                                                          <  2:395111/151‑1 (MQ=255)                                                                                                                                            cgcCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCCGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAACATCGTTTTGCCCAACGGGCAACGCTTAGCCAAGTTCCAGAGTACAGGCCATCGCGCCAAAATGTCAGGCTCTGACATgtgt                                                                                                               >  2:634101/1‑273 (MQ=255)                                                                                                                                            cgcCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGAACAGGCCAGCGCGCCACAATGTCTGGCGCTGaatgtgagaa                                                                                                             >  2:100975/1‑267 (MQ=255)                                                                                                                                                      gTTCCTTTCTCAAACGGCATAACATTCAGCGTGTCACTTGCCATCTTCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGAGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGAt                                                                                                                                                                                                                                                               >  1:361434/1‑119 (MQ=255)                                                                                                                                                          cTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCg                                                                                                                                                                                                                                                           <  2:361429/119‑1 (MQ=255)                                                                                                                                                                   aCGGCATACACTTCCGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGGGGCGAGTAATTTGCGAAACCACCCGATAACGGGGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATgg                                                                                      <  1:219525/275‑1 (MQ=255)                                                                                                                                                                     ggCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAgcggcg                                                                                                                                                                                                                                                                            <  1:434428/91‑1 (MQ=255)                                                                                                                                                                                 aGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCTGCGGTGTTCTCACGATACCGACActc                                                                                                                                                                                                                                                     >  2:459041/1‑102 (MQ=255)                                                                                                                                                                                     tgtCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCTGCGGTGTTCTCACGATACCGACACTCTCAc                                                                                                                                                                                                                                                 <  1:459047/102‑1 (MQ=255)                                                                                                                                                                                      gtCACGTGCCATAGGCATTTCGAAGGGCGGCGAGGGGCGAGTAATTTGCGAACCCCCCCGATAACGGAGCGGCGGTGTTCTTACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGc                                                                    <  2:227869/274‑1 (MQ=255)                                                                                                                                                                                          cTTGCCATATGCATTTCGAAGGACGTCGAGTTGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTAACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCccgc                                                               <  1:100976/275‑1 (MQ=255)                                                                                                                                                                                             gCCATATGCATTTCGAGGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAAGGATTTGCCAAGTTCAAGCNTACAGGc                                                                                                                                                >  1:552361/1‑191 (MQ=255)                                                                                                                                                                                                 tctgaaTTTCGAAGGGCGGCGAGTGGCGAGTAATTTGCGACACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAAGTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAgc                                                                                                                                            <  2:552354/186‑1 (MQ=255)                                                                                                                                                                                                            aaGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAgcggcg                                                                                                                                                                                                                                                                            <  1:223872/52‑1 (MQ=255)                                                                                                                                                                                                                                                        gaaaggcggGGTTATCACAATACCGACACCCTCACCAGATAGCAGCGCAGCAATTGGGTTGGCATTTACTGAAAATTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCGGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCg  <  2:636114/270‑1 (MQ=255)                                                                                                                                                                                                                                                         aGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGc                                                                                                                                                                   >  2:531781/1‑112 (MQ=255)                                                                                                                                                                                                                                                         aGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATTTCTGGCGCTGAAATg                                                                                                                  >  2:277092/1‑161 (MQ=255)                                                                                                                                                                                                                                                         aGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATg                                                                                                                  >  1:237609/1‑161 (MQ=255)                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                              GCAGAATGATTTCTGCGGGTAAGTATTAGCTTATTTTTTCGAGCATTAATCCCGCGCGTAATCCCAACGCTACCAACGGATTAGGGAATAACACATACTCTACATCATGGGTTACGGTAAAACGTTCCTCTCCGTCCTGCGCCAGCAATGTTCCTTTCTCAAACGGCATAAAATTCAGCGTGTCACTTGCCATATGCATTTCGAAGGACGGCGAGTGGCGAGTAATTTGCGAAACCACCCGATAACGGAGCGGCGGTGTTCTCACGATACCGACACTCTCACCAGATAGCAGCGCAGCAATTGCGCTGGCAGTTACTGCAAACTGGCGAAGATCGTTTTGCCCAAAGGGCAACGCTTTGCCAAGTTCCAGCGTACAGGCCAGCGCGCCAAAATGTCTGGCGCTGAAATGGGTAAACGTACCACCAGGTTCCTGATGGAACACCAGCGCCTCCAGCCCCGCCGCACCCAGCCACGTCAGAAATTTCTCGTCCCAGGGAATGTCGCGTTGCGGTAATACACCG  >  CP009273/1820187‑1820707

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