Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	687,119	A→C	40.0%	I41S (ATC→AGC)	ybeX ←	putative ion transport

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	687,119	0	A	C	40.0%	8.2 / 39.6	35	I41S (ATC→AGC)	ybeX	putative ion transport
Reads supporting (aligned to +/- strand):  ref base A (10/11);  new base C (7/7);  total (17/18)
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.10e-01

AAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCATGATGTCGCGAACGCGTTGGTCTGCGATGTCCATCACCCCTTCGAGCATATCGCGCGTATCTTCGTCGATAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACG  >  CP009273/686933‑687352                                                                                                                                                                                           |                                                                                                                                                                                                                                          aaaCGTGTGTGGTAGGTGTCTATGTTGACATTAAGTCATGCGTCCAGGGTCTGGTTGCGTTTCGGGGTAATCATCTGGGAGCGGGGGATCAGGAGGTCGCGAACGCGTTGGGCTGCGATGTCCATCACCCCTTCGAACATAGCGCGCGTATCTTCGTCGATAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTAt                                                                                                                                                   <  1:89557/275‑1 (MQ=255)                                        tCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCATGATGTCGCGAACGCGTTGGTCTGCGATGTCCATCACCCCTTCGAGCATATCGCGCGTATCTTCGTCGATAAGGTCGTTCTGCCCGGAATCACGGCTCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCGGAGTAACAGGGAGAACAATCCCTGCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGGCGCTGCTGGtgtc                                                                                                            >  2:641134/1‑275 (MQ=255)                                        tCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCATGATGTCGCGAACGCGTTGGTCTGCGATGTCCATCACCCCTTCGAGCATATCGCGCGTATCTTCGTCGATAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTaa                                                                                                                                                                                 >  1:604182/1‑206 (MQ=255)                                            ccAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCATGATGTCGCGAACGCGTTGGTCTGCGATGTCCATCACCCCTTCGAGCATATCGCGCGTATCTTCGTCGATAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAgg                                                                                                                                                                             <  2:604175/206‑1 (MQ=255)                                              aGCGTCTGGTTGCGTTTCAGGTTAATCATCTGGGAGCGGGGGATCATGATGTCGCGAACGCTTTGGTCTGCGATGTCCATCACCCCTTCGAGCATATCGCGCGTATCTTCGTCGATNAGGTCGTTCTGTCCGGAATCACGGCTCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATgg                                                                                                      <  1:641142/275‑1 (MQ=255)                                                                            tGGGAGCGGGGGATCATGATGTCGCGAACGCGTTGGTCTGCGATGTCCATCACCCCTTCGAGCATATCGCGCGTATCTTCGTCGATAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCTTCGCTCATGTCGGCGGATGGGTTGTGATGTTGGTTNNANaaaagcgcgtt                                                                        >  2:579152/1‑268 (MQ=255)                                                                                   gggggATCATGATGTCGCGAACGCGTTGGTCTGCGATGTCCATCACCCCTTCGAGCATATCGCGCGTATCTTCGTCGATAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATc                                                                 >  1:143546/1‑275 (MQ=255)                                                                                               ggctgcgATCGCGTTTGTCTGCGATTTCCATCACTCCTTCGGGCATATCGCGCGTTTCTTTGTCGTTTAGGGCTTTCTGCCCGGCATCACGGATCAGCGCCAGCAGTTCGTCACGGTTGTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAAcggcgg                                                      <  1:579159/270‑1 (MQ=255)                                                                                                                 cTGCGATGTCCATCATCCCTTCGAGCATAGCGCGCGTATCTTCGTCGATAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATNTTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCACACCGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAAt                                   >  1:506637/1‑275 (MQ=255)                                                                                                                         tCCATCACCCCTTCGAGCATATCGCGCGTATCTTCGTCGATAAGGTCGTTCTGCCCGGAATCACGGCTCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCcttc                                                                                                                                                              >  2:377735/1‑144 (MQ=255)                                                                                                                             tCACCCCTTCGAGCATATCGCGCGTATCTTCGTCGATAAGGTCGTTCTGCCCGGAAGCACGGCTCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCttgt                                                                                                                                                          <  1:377740/144‑1 (MQ=255)                                                                                                                                   cTTCGAGCATATCGCGCGTATCTTCGTCGATAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATcc                                                                                                                                                                  >  2:535504/1‑130 (MQ=255)                                                                                                                                   cTTCGAGCATATCGCGCGTATCTTCGTCGATAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGc                                                                                                                                                       >  2:431122/1‑141 (MQ=255)                                                                                                                                       gAGCATATCGCGCGTATCTTCGTCGATAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCcttc                                                                                                                                                              <  1:535510/130‑1 (MQ=255)                                                                                                                                       gAGCATATCGCGCGTATCTTCGTCGATAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTAt                                                                                                                                                   <  1:431128/141‑1 (MQ=255)                                                                                                                                         gtatttCGAGCGTATCTTCGTCGATAAGGTCGTGCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCgg           <  1:139496/270‑1 (MQ=255)                                                                                                                                             atCGCGCGTATCTTCGTCGATAAGGTCGTTCTGCCCGGAATCACGGCTCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCttgtt                                                                                                                                                         >  2:490067/1‑129 (MQ=255)                                                                                                                                               cgcgcgTATCTTCGTCGATAAGGTCGTTCTGCCCAGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTc                                                                                                                                                                                                              >  1:486743/1‑74 (MQ=255)                                                                                                                                                 cgcgTATCTTCGTCGATAAGGTCGTTCTGCCCGGAATCACGGCTCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACGGGGAGAAAAATCCCTTCTTGTTGCtt                                                                                                                                                     <  1:490073/129‑1 (MQ=255)                                                                                                                                                   cgTATCTTCGTCGATAAGGTCGTTCTGCCCGGAATCACGGCTCAgc                                                                                                                                                                                                                                      <  1:252747/46‑1 (MQ=255)                                                                                                                                                   cgTATCTTCGTCGATAAGGTCGTTCTGCCCGGAATCACGGCTCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGCCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCACGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGCCAAACGCCTTGCAAATTAtcccttcctcgcacttccct  >  2:192119/1‑261 (MQ=255)                                                                                                                                                   cgTATCTTCGTCGATAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGtt                                                                                                                                                                                                          <  2:486737/74‑1 (MQ=255)                                                                                                                                                       tCTTCGTCGATAAGGTCGTTCTGCCCGGAATCACGGCTCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCa                                                                                                                   >  1:214248/1‑157 (MQ=255)                                                                                                                                                       tCTTCGTCGATAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCa                                                                                                                                            >  1:296018/1‑132 (MQ=255)                                                                                                                                                           cgtcgATAAGGTCGTTCTGCCCGGAATCACGGCTCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGc                                                                                                               <  2:214245/157‑1 (MQ=255)                                                                                                                                                           cgtcgATAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCactac                                                                                                                                        <  2:296014/132‑1 (MQ=255)                                                                                                                                                                aTAAGGTCGTTCTGCCCGGAATCACGGCTCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCgt                                                                                                                                                >  1:261531/1‑119 (MQ=255)                                                                                                                                                                    ggTCGTTCTGCCCGGAATCACGGCTCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCa                                                                                                                                            <  2:261527/119‑1 (MQ=255)                                                                                                                                                                              cccGGAATCACGGCTCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAAtt                                                                                                                              >  2:490599/1‑123 (MQ=255)                                                                                                                                                                              cccGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCgg                                                                                                                                                                                                            >  2:647118/1‑45 (MQ=255)                                                                                                                                                                              cccGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCtt                                                                                                                                                     >  2:597812/1‑100 (MQ=255)                                                                                                                                                                                  gAATCACGGCTCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTgtcg                                                                                                                          <  1:490605/123‑1 (MQ=255)                                                                                                                                                                                  gAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCa                                                                                                                                                                                                        <  1:647126/45‑1 (MQ=255)                                                                                                                                                                                  gAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCg                                                                                                                                                 <  1:597819/100‑1 (MQ=255)                                                                                                                                                                                       aCGGATCAGCGCCAGCAGTTCGTGACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTAttctttct              <  1:448703/226‑1 (MQ=255)                                                                                                                                                                                           |                                                                                                                                                                                                                                          AAACGTGAGTGGGCGGACTCGATGATGACATCAAGACATTCGTCCAGCGTCTGGTTGCGTTTCAGGGTAATCATCTGGGAGCGGGGGATCATGATGTCGCGAACGCGTTGGTCTGCGATGTCCATCACCCCTTCGAGCATATCGCGCGTATCTTCGTCGATAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGACGTCTGCCAGTCAAATGCCTGGCAAATTATTCTTTCTCGGCAATGTACG  >  CP009273/686933‑687352

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