Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,894,259	G→C	12.3%	G161A (GGC→GCC)	PP_4302 →	urea transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,894,259	0	G	C	12.3%	83.9 / 3.3	31	G161A (GGC→GCC)	PP_4302	urea transporter
Reads supporting (aligned to +/- strand):  ref base G (12/15);  new base C (2/2);  total (14/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ACGCACCAGTGGCGCAAGCGTGGCGGCAAGTTGCTGATTGCCTACACCTCCCCTTTCGTATTGCTGGGCTGGGCCACACTGCTGGTCGCCAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCGGGTTGCTGATCATCGTCGGCATGTTCATCGTCAATCCCTATGCGGCGATGTGGGCGGTGATCGGCTCGGCCATTGGTGGCGGCGTGGCGTT  >  NC_002947/4894114‑4894388                                                                                                                                                   |                                                                                                                                  aCGCACCAGTGGCGCAAGCGTGGCGGCAAGTTGCTGATTGCCTACACCTCCCCTTTCGTATTGCTGGGCTGGGCCACACTGCTGGTCGCCAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTg                                                                                                                               <  2:540626/150‑1 (MQ=255)  cGCACCAGTGGCGCAAGCGTGGCGGCAAGTTGCTGATTGCCTACACCTCCCCTTTCGTATTGCTGGGCTGGGCCACACTGCTGGTCGCCAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTgg                                                                                                                              <  2:107967/150‑1 (MQ=255)       cAGTGGCGCAAGCGTGGCGGCAAGTTGCTGATTGCCTACACCTCCCCTTTCGTATTGCTGGGCTGGGCCACACTGCTGGTCGCCAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCAg                                                                                                                         >  1:206548/1‑150 (MQ=255)        aGTGGCGCAAGCGTGGCGGCAAGTTGCTGATTGCCTACACCTCCCCTTTCGTATTGCTGGGCTGGGCCACACTGCTGGTCGCCAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑ACGCCGCGGGCCaac                                                                                                                        <  1:224016/150‑3 (MQ=255)         gTGGCGCAAGCGTGGCGGCAAGTTGCTGATTGCCTACACCTCCCCTTTCGTATTGCTGGGCTGGGCCACACTGCTGGTCGCCAGCCCCTCGCCCAGTGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCAGAt                                                                                                                       <  2:55296/150‑1 (MQ=255)              gcAAGCGTGGCGGCCAGCTGCTGATTGCCTACACCTCCCCTTTCGTATTGCTGGGCTGGGCCACACTGCTGGTCGCCAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCAGATTTTcc                                                                                                                  >  2:315981/1‑150 (MQ=255)                               ttGCTGATTGCCTACACCTCCCCTTTCGTATTGCTGGGCTGGGCCACACTGCTGGTCGCCAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGCCCCGCTGTATGCGCTGGCCC‑ACGCCGCGGGCCaac                                                                                                                        <  2:331896/127‑3 (MQ=255)                               ttGCTGATTGCCTACACCTCCCCTTTCGTATTGCTGGGCTGGGCCACACTGCTGGTCGCCAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGCTGGCCC‑ACGCCGCGGGCCaac                                                                                                                        >  1:331896/1‑125 (MQ=255)                                    gATTGCCTACACCTCCCCTTTCGTATTGCTGGGCTGGGCCACACTGCTGGTCCCCAGCCCCTCGCCCAGCGGGTTCGTCGAATCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCgg                                                                                            <  1:224048/150‑1 (MQ=255)                                         ccTACACCTCCCCTCTCGTATTGCTGGGCTGGGCCACACTGCTGGTCGCCAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCa                                                                                                                          >  2:47508/1‑115 (MQ=255)                                         ccTACACCTCCCCTCTCGTATTGCTGGGCTGGGCCACACTGCTGGTCGCCAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCa                                                                                                                          <  1:47508/115‑1 (MQ=255)                                          cTACACCTCCCCTTTCGTATTGCTGGGCTGGGCCACACTGCTGGTCGCCAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCGGGtt                                                                                         <  1:348568/147‑1 (MQ=255)                                          cTACACCTCCCCTTTCGTATTGCTGGGCTGGGCCACACTGCTGGTCGCCAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCGGGtt                                                                                         >  2:348568/1‑147 (MQ=255)                                            acacCTCCCCTTTCGTATTGCTGGGCTGGGCCACACTGCTGGTCGCCAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCGGGTTGCTGa                                                                                    <  2:206548/150‑1 (MQ=255)                                             cacCTCCCCTTTCGTATTGCTGGGCTGGGCCACACTGCTGGTCGCCAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCGGGTTGCTGat                                                                                   <  1:315981/150‑1 (MQ=255)                                                 tCCCCTTTCGTATTGCTGGGCTGGGCCACACTGCTGGTCGCCAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCGGGTTGCTGatcatc                                                                               >  2:11518/1‑150 (MQ=255)                                                  ccccTTTCGTATTGCTGGGCTGGGCCACACTGCTGGTCGCCAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCGGGTTGCTGATCAtcg                                                                              >  1:55479/1‑150 (MQ=255)                                                        tCGTATTGCTGGGCTGGGCCACACTGCTGGTCGCCAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGAACAGCTCCCGCCCTTGGCCCGGCAGCGT‑GGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCGGGTTGCTGATCATCGTCGGCa                                                                        <  1:383196/150‑1 (MQ=255)                                                                     ctgggcCACACTGCTGGTCGCCAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCCCCCGGGGCCAGATTTTCCTGCTCGCCCAGCCCTTGCCCGGGTTGCTGATCATCGTCGGCATGTTCATCGTCaa                                                           >  1:291536/1‑150 (MQ=255)                                                                         gcCACACTGCTGGTCGCCAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCGGGTTGCTGATCATCGTCGGCATGTTCATCGTCAATccc                                                       >  2:179402/1‑150 (MQ=255)                                                                                tgctgGTCGCCAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCGGGTTGCTGATCATCGTCGGCATGTTCATCGTCAATCCCTATgcgg                                                >  2:462787/1‑150 (MQ=255)                                                                                          cAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCggg                                                                                           <  1:409192/97‑1 (MQ=255)                                                                                          cAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCggg                                                                                           >  2:409192/1‑97 (MQ=255)                                                                                               ccTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCGGGTTGCTGATCATCGTCGGCATGTTCATCGTCAATCCCTATGCGGCGATGTGGGCGGTGa                                 >  2:231244/1‑150 (MQ=255)                                                                                               ccTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCGGGTTGCTGATCATCGTCGGCATGTTCATCGTCAATCCCTATGCGGCGATGTGGGCGGTGa                                 <  2:391987/150‑1 (MQ=255)                                                                                                      cAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCGGGTTGCTGATCATCGTCGGCATGTTCATCGTCAATCCCTATGCGGCGATGTGGGCGGTGAtcggctc                          <  1:179402/150‑1 (MQ=255)                                                                                                       aGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCGGGTTGCTGATCATCGTCGGCATGTTCATCGTCAATCCCTATGCGGCGATGTGGGCGGTGAtcggctcg                         <  2:401913/150‑1 (MQ=255)                                                                                                       aGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCGGGTTGCTGATCATCGTCGGCATGTTCATCGTCAATCCCTATGCGGCGATGTGGGCGGTGAtcggctcg                         <  2:55479/150‑1 (MQ=255)                                                                                                                              gCTGTATGCGTTGGCCC‑TCGGCGTGGGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCGGGTTGCTGAGCATCGTCGGCATGTTCGTCGTCAATCCCTATGCGGCGATGTGGGCGGTGATCGGCTCGGCCATTGGTGGCGGCGTGGCGtt  >  1:339109/1‑150 (MQ=255)                                                                                                                                           gCCC‑GCGGCGTGGGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCGGGTTGCTGATCATCGTCGGCATGTTCATCGTCAATCCCTATGCGGCGATGTGGGCGGTGATCGGCTCGGCCATtggt                >  1:356357/1‑123 (MQ=255)                                                                                                                                           gCCC‑GCGGCGTGGGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCGGGTTGCTGATCATCGTCGGCATGTTCATCGTCAATCCCCATGCGGCGATGTGGGCGGTGATCGGCTCGGCCATtggt                <  2:356357/123‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                  ACGCACCAGTGGCGCAAGCGTGGCGGCAAGTTGCTGATTGCCTACACCTCCCCTTTCGTATTGCTGGGCTGGGCCACACTGCTGGTCGCCAGCCCCTCGCCCAGCGGGTTCGTCGAAGCCGACCCGCTGTATGCGTTGGCCC‑GCGGCGTGGGCCAGATTTTCCTGCTCGACCAGCCCTTGGCCGGGTTGCTGATCATCGTCGGCATGTTCATCGTCAATCCCTATGCGGCGATGTGGGCGGTGATCGGCTCGGCCATTGGTGGCGGCGTGGCGTT  >  NC_002947/4894114‑4894388

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