Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	325,250	G→A	10.0%	L178L (CTG→CTA)	oprQ →	outer‑membrane porin D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	325,250	0	G	A	10.0%	102.3 / 3.4	40	L178L (CTG→CTA)	oprQ	outer‑membrane porin D
Reads supporting (aligned to +/- strand):  ref base G (19/17);  new base A (2/2);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

GCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGAAGTCGATCAACGTCTTAGGTGGCAGCTACAAGTTCACCGACAACTTCACCGCTTCGCTGTACACCGCCGACAACGAAGATGTGAT  >  NC_002947/325101‑325399                                                                                                                                                      |                                                                                                                                                      gctgTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGCCGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTg                                                                                                                                                       <  1:588109/150‑1 (MQ=255)          ttCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACgccggc                                                                                                                                              >  1:119364/1‑150 (MQ=255)                    ccAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCAccc                                                                                                                                    <  1:233432/150‑1 (MQ=255)                           cGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAgc                                                                                                                             <  1:217361/150‑1 (MQ=255)                           cGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAgc                                                                                                                             <  1:423967/150‑1 (MQ=255)                                     tACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAgagc                                                                                                                   >  1:678713/1‑150 (MQ=255)                                        ggTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGc                                                                                                                                           <  2:449182/123‑1 (MQ=255)                                        ggTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGc                                                                                                                                           >  1:449182/1‑123 (MQ=255)                                        ggTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAgagcgcc                                                                                                                <  2:678713/150‑1 (MQ=255)                                                gATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCCCCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGagcgccgagcgccg                                                                                                        >  1:270401/1‑150 (MQ=255)                                                gATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCAATTGACAGCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGagcgccgagcgccg                                                                                                        >  1:251158/1‑150 (MQ=255)                                                    ccAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAt                                                                                                    >  2:416275/1‑150 (MQ=255)                                                    ccAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTAAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAt                                                                                                    >  2:549667/1‑150 (MQ=255)                                                    ccAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTAAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGAt                                                                                                    >  2:549642/1‑150 (MQ=255)                                                             ccAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAgagcgccgagcgcc                                                                                                         >  1:367720/1‑136 (MQ=255)                                                             ccAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAgagcgccgagcgcc                                                                                                         <  2:367720/136‑1 (MQ=255)                                                             ccAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGc                                                                                           >  1:511907/1‑150 (MQ=255)                                                                  gctgcAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTAAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGaa                                                                                      <  1:549642/150‑1 (MQ=255)                                                                  gctgcAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTAAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGaa                                                                                      <  1:549667/150‑1 (MQ=255)                                                                               gaTGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGAAGTCGATCAACGTc                                                                         >  1:719614/1‑150 (MQ=255)                                                                                  gCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGAAGTCGATCAACGTCTTa                                                                      >  1:659915/1‑150 (MQ=255)                                                                                     cGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGAAGTCGATCAACGTCTTAGGt                                                                   <  1:299884/150‑1 (MQ=255)                                                                                        ctgctgCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGAAGTCGATCAACGTCTTAGGTGGc                                                                <  2:251158/150‑1 (MQ=255)                                                                                         tgctgcCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGAAGTCGATCAACGTCTTAGGTGGCa                                                               <  2:659915/150‑1 (MQ=255)                                                                                          gctgcCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGAAGTCGATCAACGTCTTAGGTGGCAg                                                              <  2:719614/150‑1 (MQ=255)                                                                                                aGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGAAGTCGATCAACGTCTTAGGTGGCAGCTACaa                                                        <  1:115628/150‑1 (MQ=255)                                                                                                     gCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGAAGTCGATCAACGTCTTAGGTGGCAGCTACAAGTTCa                                                   >  2:524188/1‑150 (MQ=255)                                                                                                      cTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAgagcgccgagcg                                                                                                           <  2:305033/93‑1 (MQ=255)                                                                                                      cTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAgagcgccgagcg                                                                                                           >  1:305033/1‑93 (MQ=255)                                                                                                                      aTCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGAAGTCGATCAACGTCTTAGGTGGCAGCTACAAGTTCACCGACAACTTCACCGCt                                  >  2:273108/1‑150 (MQ=255)                                                                                                                       tCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCt                                                                                         <  1:321776/94‑1 (MQ=255)                                                                                                                       tCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCt                                                                                         >  2:321776/1‑94 (MQ=255)                                                                                                                           cTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGAAGTCGATCAACGTCTTAGGTGACAGCTACAAGTTCACCGACAACTTCACCGCTTCGCt                             >  2:295260/1‑150 (MQ=255)                                                                                                                            tCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGAAGTCGATCAACGTCTTAGGTGGCAGCTACAAGTTCACCGACAACTTCACCGCTTCGCTg                            >  2:467596/1‑150 (MQ=255)                                                                                                                              cAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGAAGTCGATCAACGTCTTAGGTGGCAGCTACAAg                                                       >  2:505515/1‑121 (MQ=255)                                                                                                                              cAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGAAGTCGATCAACGTCTTAGGTGGCAGCTACAAg                                                       <  1:505515/121‑1 (MQ=255)                                                                                                                                             ccTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGAAGTCGATCAACGTCTTAGGTGGCAGCTACAAGTTCACCGACAACTTCACCGCTTCGCTGTACACCGCCGACAACGa           <  1:273108/150‑1 (MQ=255)                                                                                                                                               tGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGAAGTCGATCAACGTCTTAGGTGGCAGCTACAAGTTCACCGACAACTTCACCGCTTCGCTGTACACCGCCGACAACGAAg         >  1:524854/1‑150 (MQ=255)                                                                                                                                                     tGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGAAGTCGATCAACGTCTTAGGTGGCAGCTACAAGTTCACCGACAACTTCACCGCTTCGCTGTACACCGCCGACAACGAAGa        <  1:434624/145‑1 (MQ=255)                                                                                                                                                     tGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGAAGTCGATCAACGTCTTAGGTGGCAGCTACAAGTTCACCGACAACTTCACCGCTTCGCTGTACACCGCCGACAACGAAGa        >  2:434624/1‑145 (MQ=255)                                                                                                                                                     tGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGAAGTCGATCAACGTCTTAGGTGGCAGCTACAAGTTCACCGACAACTTCACCGCTTCGCTGTACACCGCCGACAACGAAGATGtga   <  1:295260/150‑1 (MQ=255)                                                                                                                                                      gAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGAAGTCGATCAACGTCTTAGGTGGCAGCTACAAGTTCACCGACAACTTCACCGCTTCGCTGTACACCGCCGACAACGAAGATGtgat  <  2:119364/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GCTGTGAAGTTCCGCATTTCCAACACCGTGCTGAAGTACGGTGACCAGATGCCAGCACTGCCAGTGCTGCAGTACGACGATGCACGTCTGCTGCCAGAAAGCTTCACCGGCACCCTGATCACCTCCAAAGAGATCCAGGGCCTGGAACTGAACGCCGGCCGCTTCACCCAGGAAGCGCGCAAGAGCGCCGAGCGCCGTGATGGTGGTGGCCTGAAGTCGATCAACGTCTTAGGTGGCAGCTACAAGTTCACCGACAACTTCACCGCTTCGCTGTACACCGCCGACAACGAAGATGTGAT  >  NC_002947/325101‑325399

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