Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	2,305,946	G→T	77.8%	R95S (CGT→AGT)	ompC ←	outer membrane porin protein C

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	2,305,946	0	G	T	77.8%	48.0 / 11.5	31	R95S (CGT→AGT)	ompC	outer membrane porin protein C
Reads supporting (aligned to +/- strand):  ref base G (3/4);  new base T (10/14);  total (13/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 = 5.58e-01

GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCAAGACGCATGTAGGTCTGGTCGCCATCT  >  CP009273/2305802‑2306079                                                                                                                                                 |                                                                                                                                      gcATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGCAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACGCACATCCTGGAATTTCAGGCCTGGGTGTGCCGCACTGGTc                                                                                                                                   >  2:253086/1‑149 (MQ=255)        gTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACTGGTCCAGGAgt                                                                                                                            <  2:367483/149‑1 (MQ=255)         ttGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCGGGCgtt                                                                                                                           >  2:47413/1‑149 (MQ=255)          tcttcgacccGGAGGTGGCACCACCGAATTCGGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAg                                                                                                                                <  1:57930/135‑1 (MQ=255)                aaCCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACTGGTCCAGGAgttg                                                                                                                          >  1:196609/1‑143 (MQ=255)                    gTAGGTGTCGCCACCGAATTCTGGGAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACTGGTCCAGGAGTTgtttt                                                                                                                      <  1:90857/143‑1 (MQ=255)                               caccGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACTGGTCCAGGAGTTGTTTTCGTTTTCAgcg                                                                                                           >  1:290368/1‑143 (MQ=255)                                    aaTTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATATACAACGCCATACATACGTCCGTAGCCGAATTAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACTGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTg                                                                                                <  2:190556/149‑1 (MQ=255)                                         tGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTTGAATTTCAGACCTGCGAATGCCACACTGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCt                                                                                           <  2:35825/149‑1 (MQ=255)                                           gaaGTACGTCGGTCCAGGAAGTTACGCCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGa                                                                                         <  2:301153/147‑1 (MQ=255)                                           gCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACTGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGa                                                                                         <  2:284822/149‑1 (MQ=255)                                               taggtcgtgccATGAAGGTACGGCATAAACAACGGCGTAGTTACGTCCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACTGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCt                                                                                           <  1:365377/134‑1 (MQ=255)                                               tACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACTGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCt                                                                                           >  1:318125/1‑143 (MQ=255)                                               tACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACTGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCt                                                                                           >  1:105623/1‑143 (MQ=255)                                                         cAGGCAGTTCCGTTACCATCAACGCCGCAGTTACGCCCGCAGTCGCAACAACCCCCCTCCTGGAATTTCAGACCTGCGAATGCCACACTGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGcc                                                                           <  2:198305/149‑1 (MQ=255)                                                          aGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACTGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATCTTCCCACTGGCCg                                                                          >  2:32704/1‑149 (MQ=255)                                                             aaGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACTGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTgg                                                                             <  1:120900/143‑1 (MQ=255)                                                                ctaggTCATAAACAACGCCGTAGTGACGACCGCAGTCGCAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACTGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATCTTCCCACTGGCCg                                                                          <  1:366794/139‑1 (MQ=255)                                                                         aaaCAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACTGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCCGCTgg                                                           >  2:402584/1‑149 (MQ=255)                                                                                gCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGt                                                          <  1:385382/143‑1 (MQ=255)                                                                                       ttACGACCGTAGTCGAAAGAACCAATATCTTGACTGATGAGACCTTCGAATCCCACCCTAGCCCAAGACTTTTTTTCGTTTTCAGCTCTGCGGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAAc                                                   <  1:222337/143‑1 (MQ=255)                                                                                           ggccGTAGTCGAAAGAACCCACATCCTGGAAATTCAGCCCTGCGAATGCCACACTGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCa                                         <  2:119409/147‑1 (MQ=255)                                                                                           gACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACTGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCa                                         <  2:356423/149‑1 (MQ=255)                                                                                           gACCGTAGCCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGa                                               <  1:255902/143‑1 (MQ=255)                                                                                                   tCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGaa                                 >  2:149154/1‑149 (MQ=255)                                                                                                     gAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACTGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGGCAGTAACCTGAGTTTAACCtt                                     >  1:313014/1‑143 (MQ=255)                                                                                                      aaaGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACTGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGcc                              >  2:361328/1‑149 (MQ=255)                                                                                                             ccACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCa                             >  1:4877/1‑143 (MQ=255)                                                                                                                       ggaTTTCAGACCTGCGAATGCCACACTGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCAAGACGCATGt                   <  1:304469/141‑1 (MQ=255)                                                                                                                             cAGACCTGCGAATGCCACACTGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCAAGACGCATGTAGGTCTGGTCGc       <  2:326254/149‑1 (MQ=255)                                                                                                                                  cTGCGAATGCCACACTGGTCCAGGAGTTGATTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCGCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCAAAACGCATGTAGGTCTGGTCGCCATCt  >  2:263805/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                      GCATGAAGTTGTCAGAACCGTAGGTGTCACCACCGAATTCTGGCAGTACGTCGGTCCAGGAAGTTACGTCATAAACAACGCCGTAGTTACGACCGTAGTCGAAAGAACCCACATCCTGGAATTTCAGACCTGCGAATGCCACACGGGTCCAGGAGTTGTTTTCGTTTTCAGCGCTGTTGCCCTGGATCTGATATTCCCACTGGCCGTAACCGGTCAGCTGGTCAGTAACCTGAGTTTCACCTTTGAAGCCAAGACGCATGTAGGTCTGGTCGCCATCT  >  CP009273/2305802‑2306079

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