Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP006881	1,044,161	T→G	23.7%	S72R (AGC→CGC)	U712_05130 ←	Putative amidohydrolase yhaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP006881	1,044,161	0	T	G	23.7%	67.1 / 4.1	44	S72R (AGC→CGC)	U712_05130	Putative amidohydrolase yhaA
Reads supporting (aligned to +/- strand):  ref base T (13/20);  new base G (9/2);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 3.39e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGAGCTGCGGTGTGGCCGTCGTGGCCGCATGCATGCATGACACCAGGCACTTTGGAGGCGTAAGGGACATCTTTTTCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTGGGACTCCTAACGATTCATAATAGGAAGCAATAAATGCGGCGGTTTTTTCTTCTTGAAATGAGAGCTCAGGATACATATGA  >  CP006881/1044023‑1044291                                                                                                                                           |                                                                                                                                   gaagaGGGGCGGTGTGGCGGGCGTGGCCGCATGCATGCATGACACCCGGCAGTTTGGCGGCGTAAGGGACATATTTTTCATATTGAAGGGGGAGCCCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGct                                                                                                                                    <  4:174686/139‑1 (MQ=255)      gCTGCGGTGTGGCCGTCGTGGCCGCATGCATGCATGACACCAGGCACTTTGGAGGCGTAAGGGACATCTTTTTCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGcttcct                                                                                                                                >  7:26319/1‑138 (MQ=255)       cTGCGGTGTGGCCGTCGTGGCCGCATGCATGCATGACACCAGGCACTTTGGAGGCGTAAGGGACATCTTTTTCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGcttccgtc                                                                                                                              >  6:371799/1‑139 (MQ=255)            ggggggCGGTCGTGGCCGAATGAATGAATGACACCAGGCACTTTGGAGGCTTACGGGACTTCTTTTTCATTTTTATTGGGGAGCGCGTCAACGGCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGCTTCCTTCTatat                                                                                                                         <  4:256441/135‑1 (MQ=255)                 gCCGTCGTGGCCGCATGCATGCATGACACCAGGCACTTTGGAGGCGTAAGGGACATCTTTTTCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGCGTCGGTGTATATTTgg                                                                                                                     >  3:168985/1‑137 (MQ=255)                     tcgtGGCCGCATGCATGCATGACACCAGGCACTTTGGAGGCGTAAGGGACATCTTTTTCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTGCTaaaa                                                                                                                >  8:93299/1‑139 (MQ=255)                     ggggggCCGCATGCATGATTGCACCCCGGCATTTTGGTGGCGTAGGGGCCATCTTTTTCATCTTGAATGGGGCCCGCGGCAACGGCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTGCTaaaa                                                                                                                <  4:400938/135‑1 (MQ=255)                               atgcatgcatgcCACCAGGCACTTTGGAGGCGTAAGGGACATCTTTTTCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGCTTGCTTGGATTTTTTATAAAACCCCTCTAcc                                                                                                      >  7:293329/1‑139 (MQ=255)                                   atgcatgACACCAGGCACTTTGGAGGCGTAAGGGACATCTTTTTCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAACGCGACTCTAGGGCCGGGTTCGCTTCCTTCTAAATTTGCGAAAACCCCACAACCGcc                                                                                                   >  8:99347/1‑138 (MQ=255)                                         gACACCAGGCACTTTGGAGGCGTAAGGGACATCTTTTTCATCTTGAATGGGGCGAGCGTCAAAGTCGGCCCTCAAAGCGCCTGGAGGGCCGGGTTCGCTTCCTTCTATATTTGCTTAAACCCCGCTACCGCCCACGgtt                                                                                            >  8:36643/1‑139 (MQ=255)                                         gACACCAGGCACTTTGGAGGCGTAAGGGACATCTTTTTCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGGTAGCGTACTTCTATATTTTCTTAAAACCCACTAACGCCCACGttt                                                                                            >  2:258691/1‑139 (MQ=255)                                              cAGGCACTTTGGAGGCGTAAGGGACATCTTTTTCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGTTGGCTTGCGTGGATATTTTCTAAAACCCCTCAACCCCCAACGTTtgtggg                                                                                       >  1:86705/1‑137 (MQ=255)                                               aGGCACTTTGGAGGCGTAAGGGACATCTTTTTCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGCGTCCTTCTATATTTGCTAAAACCCCTCTACCGCCCACCGTTGTGCgg                                                                                      >  6:97056/1‑139 (MQ=255)                                                    cTTTGGCGCCGTATGGGAATTCTTTTTAATATTGATTGGGGAGAGGGGCAATGGCGGCCCTAAAAGCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTgg                                                                                 <  3:142528/139‑1 (MQ=255)                                                     tttGGAGGCGTAAGGGACATCTTTTTCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGCCTGTAGGGCCGGGTTGGCGTCGCCGTAGATATGCGAAAAACCCACTACCCCCAACCATGGTGCggctagtg                                                                                >  2:71802/1‑133 (MQ=255)                                                     tttGGAGGCGTAAGGGACATCTTTTTCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTggg                                                                                >  2:22667/1‑139 (MQ=255)                                                     tttGGAGGCGTAAGGGACATCTTTTTAATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTACGGCCGGGTTCGGGTACGTGTATAATTTCATAATCACCAAGACCCCCAACGGTTTTGCGGATTggg                                                                                >  4:368258/1‑139 (MQ=255)                                                       tGGAGGCGTAAGGGACATCTTTTTCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTGGGa                                                                               <  7:63034/138‑1 (MQ=255)                                                             cGTAAGGGACATCTTTTTCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTGGGACTCCTaa                                                                        <  1:22667/139‑1 (MQ=255)                                                             cGTAAGGGACATCTTTTTCATCTTGAATGGGGAGAGCGGCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTGGGACTCCTaa                                                                        <  8:263079/139‑1 (MQ=255)                                                               tAAGGGACATCTTTTTCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTTCTAAAACCCCACTACCGCCAACGTTTGTGCGGATTGGAACAACTAAc                                                                       >  4:332827/1‑138 (MQ=255)                                                               tAAGGGACATCTTTTTCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGGTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACCTTTGTGCGGTTTGGGACTCCTAAc                                                                       >  7:95293/1‑138 (MQ=255)                                                               tAAGGCACATCTTTTTCAGCTTGATTGGGCACACCGACAAAGGCGCCCCTCAAAGCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTGGGACTCCTAACg                                                                      <  7:111680/139‑1 (MQ=255)                                                                 cgggACATATTTTTAATCTTCACTGGGGCCCGCGAAAACGGCGGCCCCAACAGCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTGGGACTCCTAACGAt                                                                    <  6:293094/138‑1 (MQ=255)                                                                 aGGGACATCTTTTTCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTGGGACTCCTAACGAt                                                                    >  7:31717/1‑139 (MQ=255)                                                                   ggACATCTTTTTCATGTTGATTGGGGAGAGGGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGGTTGTGCGGATTGGGACTCCTAACGAtt                                                                   <  2:12346/138‑1 (MQ=255)                                                                      cATCTTTTTCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGGTCGCTTACATCTATAATTTCTAAAACCCCTCAAACCCCAACGTTTTTGCGGATTTGGACAACAAACCATTCat                                                                >  8:369808/1‑138 (MQ=255)                                                                       aTCTTTTTCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAAACCCTCTACCGCCAACGTTTGTGCGGATTGGGACTCCTAACGATTCataa                                                              >  1:267988/1‑139 (MQ=255)                                                                       aTCTTTTTCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGGTCGCGTCCGTGTATATTTTGTAAAACCCCGCAACCGCCCACGTTTGTGCGGGTTGGGACTCCTAAACATTCaaa                                                               >  7:163233/1‑136 (MQ=255)                                                                             ttCATCTTGATTGGGGAGCGGTTCAAAGTCGGCCCCAAAAGCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTGGGACTCCTAACGATTCATAATAGGa                                                         <  8:163233/138‑1 (MQ=255)                                                                             ttCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGCGTCCTTGTATATTTTCTAAAAACCCCCAACCCCCAACCTTTTTGGGGCGTTGGACGGCAAACGAAtaataaaaatca                                                        >  5:42467/1‑133 (MQ=255)                                                                                         ggggCGCGCGCCAACGGCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAACACCCCTCTACCGCCAACGTTTGTGCGGATTGGGACTCCTAACGATTCATAATAGGAAGCAATAAATgcg                                            <  6:69483/139‑1 (MQ=255)                                                                                         ggggAGAGCGTCAAAGTCGGCCCTCAACGCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTGGGACTCCTAACGATTCATAATAGGAAGCAATAAATgcg                                            <  7:93299/139‑1 (MQ=255)                                                                                          gggAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGGTCGCGTCCTTCTATATTTTCTAAAACCCCACTACCGCCAACGTTTGTGCGGATTGGGACTCCTAACGATTCATAATAGGAAGCAATAAATgcgg                                           >  1:123575/1‑139 (MQ=255)                                                                                                gCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTGGGACTCCTAACGATTCATAATAGGAAGCAATAAATGCGGCGGttt                                     >  4:199568/1‑139 (MQ=255)                                                                                                                  acaGCGACTGTAGGGCCGGGTTCGCGTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTGGGACTCCTAACGATTCATAATAGGAAGCAATAAATGCGGCGGTTTTTTCTTCTTGAAATga                     <  1:250921/135‑1 (MQ=255)                                                                                                                  aaaGCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTGGGACTCCTAACGATTCATAATAGGAAGCAATAAATGCGGCGGTTTTTTCTTCTTGAAATga                     <  2:267988/137‑1 (MQ=255)                                                                                                                     gCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTGGGACTCCTAACGATTCATAATAGGAAGCAATAAATGCGGCGGTTTTTTCTTCTTGAAATGAGAGCt                <  8:26319/139‑1 (MQ=255)                                                                                                                     gCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTGGGACTCCTAACGATTCATAATAGGAAGCAATAAATGCGGCGGTTTTTTCTTCTTGAAATGAGAGCt                <  5:97056/139‑1 (MQ=255)                                                                                                                            tAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTGGGACTTCTAACGATTCATAATAGGAAGCAATAAATGCGGCGGTTTTTTCTTCTTGAAATGAGAGCTCAGGATa         <  3:175566/139‑1 (MQ=255)                                                                                                                            tAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTGGGACTCCTAACGATTCATAATAGGAAGCAATAAATGCGGCGGTTTTTTCTTCTTGAAATGAGAGCTCAGGa           <  7:369808/137‑1 (MQ=255)                                                                                                                            tAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTGGGACTCCTAACGATTCATAATAGGAAGCAATAAATGCGGCGGTTTTTTCTTCTTGAAATGAGAGCTCAGGa           <  1:272679/137‑1 (MQ=255)                                                                                                                            tAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTGGAACTCCTAACGATTCATAATAGGAAGCAATAAATGCGGCGGTTTTTTCTTCTTGAAATGAGAGCTCAGGa           <  3:332827/137‑1 (MQ=255)                                                                                                                            tAGGGCCGGGTTCGCTTCCTTCTAAATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTGGGACTCCTAACGATTCATAATAGGAAGCAATAAATGCGGCGGTTTTTTCTTCTTGAAATGAGAGCTCAGGATa         <  4:12051/139‑1 (MQ=255)                                                                                                                            tAGGGCCGGGTTAGCGTCCTTCTATATTTGCTAAAACCCCGCTACCGCCAACGTTTGTGCGGATTGGGACTCCTAACGATTCATAATAGGAAGCAATAAATGCGGCGGTTTTTTCTTCTTGAAATGAGAGCTCAGGATa         <  6:96249/139‑1 (MQ=255)                                                                                                                                    ggTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTGGGACTCCTAACGCTTCATAATAGGAAGCAATAAATGCGGCGGTTTTTTGTTTTTGAAATGTGCGGGCAGGGTACAAATgt  >  8:36847/1‑137 (MQ=255)                                                                                                                                    ggTTCGCGTCCTTCTATATTTGCTAAAACCCCTCTACCCCCAACGTTTGTTCGGATTGGGACTACTAACGATTCATAATAGGAAGCAATAAATGCGGCGGTTTTTTCTTAATTAAATAACTGCTCATCATACctatta  >  6:89627/1‑136 (MQ=255)                                                                                                                                           |                                                                                                                                   GAAGAGCTGCGGTGTGGCCGTCGTGGCCGCATGCATGCATGACACCAGGCACTTTGGAGGCGTAAGGGACATCTTTTTCATCTTGAATGGGGAGAGCGTCAAAGTCGGCCCTCAAAGCGACTGTAGGGCCGGGTTCGCTTCCTTCTATATTTGCTAAAACCCCTCTACCGCCAACGTTTGTGCGGATTGGGACTCCTAACGATTCATAATAGGAAGCAATAAATGCGGCGGTTTTTTCTTCTTGAAATGAGAGCTCAGGATACATATGA  >  CP006881/1044023‑1044291

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