Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,239,048	T→C	S384P (TCA→CCA)	lamB →	maltoporin LamB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,239,048	0	T	C	82.6%	99.1 / 20.4	46	S384P (TCA→CCA)	lamB	maltoporin LamB
Reads supporting (aligned to +/- strand):  ref base T (4/4);  new base C (20/18);  total (24/22)
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.48e-01

CAAGTGGACGCCAATCATGAGCACCGTGATGGAAATCGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGTCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAATGGGGTTACGACTACACCGGTAACGCTGATAACAACGCGAACTTCGGCAAAGCCGTTCCTGCTGATTTCAACGGCG  >  NZ_CP009273/4238915‑4239174                                                                                                                                      |                                                                                                                               cAAGTGGACGCCAATCATGAGCACCGTGATGGAAATCGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGc                                                                                                                           >  1:581819/1‑139 (MQ=255)     tGGACGCCAATCATGAGCACCGTGATGGAAATCGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCgg                                                                                                                       <  2:393465/139‑1 (MQ=255)       gACGCCAATCATGAGCACCGTGATGGAAATCGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGTCACGCCCGGc                                                                                                                      >  2:224760/1‑138 (MQ=255)       gACGCCAATCATGAGCACCGTGATGGAAATCGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACg                                                                                                                            >  2:97375/1‑132 (MQ=255)       gACGCCAATCATGAGCACCGTGATGGAAATCGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACg                                                                                                                            <  1:97375/132‑1 (MQ=255)       gACGCCAATCATGAGCACCGTGATGGAAATCGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGc                                                                                                                      >  2:195957/1‑138 (MQ=255)       gACGCCAATCATGAGCACCGTGATGGAAATCGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGc                                                                                                                      >  1:118959/1‑138 (MQ=255)              aTCATGAGCACCGTGATGGAAATCGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCgtg                                                                                                              <  2:151987/139‑1 (MQ=255)               tCATGAGCACCGTGATGGAAATCGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCgtgt                                                                                                             >  2:745167/1‑139 (MQ=255)                    aGCACCGTGATGGAAATCGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGTCACGCCCGGCTATTCGTGTCTTCg                                                                                                        <  2:524305/139‑1 (MQ=255)                    aGCACCGTGATGGAAATCGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCg                                                                                                        <  2:232649/139‑1 (MQ=255)                      cACCGTGATGGAAATCGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCg                                                                                                        <  2:771700/137‑1 (MQ=255)                      cACCGTGATGGAAATCGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCg                                                                                                        >  1:771700/1‑137 (MQ=255)                      cACCGTGATGGAAATCGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGc                                                                                                       >  2:7154/1‑138 (MQ=255)                           tGATGGAAATCGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTa                                                                                                 >  1:252383/1‑139 (MQ=255)                                  aaTCGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAg                                                                                          <  2:478689/139‑1 (MQ=255)                                   aTCGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGt                                                                                         >  1:141258/1‑139 (MQ=255)                                   aTCGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGt                                                                                         >  1:8208/1‑139 (MQ=255)                                     cGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTgg                                                                                       >  1:851169/1‑139 (MQ=255)                                     cGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTgg                                                                                       >  2:181496/1‑139 (MQ=255)                                        cTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGAt                                                                                    <  2:650651/139‑1 (MQ=255)                                           cGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGAt                                                                                    >  2:809741/1‑136 (MQ=255)                                               aaCGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAAt                                                                             <  2:851169/139‑1 (MQ=255)                                                     gAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAATGGGGt                                                                        <  1:312955/138‑1 (MQ=255)                                                     gAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAATGGGGt                                                                        <  2:581819/138‑1 (MQ=255)                                                     gAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAATGGGGt                                                                        <  2:789280/138‑1 (MQ=255)                                                          ccAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAATGGGGTTACGAc                                                                  <  1:421679/139‑1 (MQ=255)                                                            aGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGt                                                                                                                                <  1:844642/75‑1 (MQ=255)                                                            aGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGt                                                                                                                                >  2:844642/1‑75 (MQ=255)                                                                  ccGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAATGGGGTTACGACTACACCgg                                                          <  2:141258/139‑1 (MQ=255)                                                                   cGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAATGGGGTTACGACTACACCgg                                                          >  1:480501/1‑138 (MQ=255)                                                                    ggCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGTCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAATGGGGTTACGACTACACCGGTa                                                        <  2:637392/139‑1 (MQ=255)                                                                    ggCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGTCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAATGGGGTTACGACTACACCGGTa                                                        <  2:100372/139‑1 (MQ=255)                                                                    ggCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAATGGGGTTACGACTACACCGGTa                                                        <  1:236309/139‑1 (MQ=255)                                                                                    aGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAATGGGATTACGACTACACCGGTAACGCTGATAACAAcg                                         <  1:745167/138‑1 (MQ=255)                                                                                      tACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGt                                                                                         >  2:312120/1‑88 (MQ=255)                                                                                      tACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGt                                                                                         <  1:312120/88‑1 (MQ=255)                                                                                      tACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAATGGGGTTACGACTACACCGGTAACGCTGATAACAAcgcg                                       >  1:334350/1‑138 (MQ=255)                                                                                       aCAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAATGGGGTTACGACTACACCGGTAACGCTGATAACAACGCGaa                                     >  2:436351/1‑139 (MQ=255)                                                                                             ttACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAATGGGGTTACGACTACACCGGTAACGCTGATAACAACGCGAACTTCgg                               <  2:334350/139‑1 (MQ=255)                                                                                                     gCACAACAATGGCAGGCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAATGGGGTTACGACTACACCGGTAACGCTGATAACAACGCGAACTTCGGCAAAGCCg                       >  1:186558/1‑139 (MQ=255)                                                                                                              tGGCAGGCTGGCGACAGCATCTGGTCACGCCCGGCTATTCGTGTCTTCg                                                                                                        <  2:800176/49‑1 (MQ=255)                                                                                                              tGGCAGGCTGGCGACAGCATCTGGTCACGCCCGGCTATTCGTGTCTTCg                                                                                                        >  1:800176/1‑49 (MQ=255)                                                                                                              tGGCAGGCTGGCGACAGCATCTGGTCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAATGGGGTTACGACTACACCGGTAACGCTGATAACAACGCGAACTTCGGCAAAGCCGTTCctgct               >  2:779995/1‑138 (MQ=255)                                                                                                                   ggCTGGCGACAGCATCTGGTCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAATGGGGTTACGACTACACCGGTAACGCTGATAACAACGCGAACTTCGGCAAAGCCGTTCCTGCTGAttt          >  1:359326/1‑138 (MQ=255)                                                                                                                   ggCTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAATGGGGTTACGACTACACCGGTAACGCTGATAACAACGCGAACTTCGGCAAAGCCGTTCCTGCTGAttt          >  1:515622/1‑138 (MQ=255)                                                                                                                     cTGGCGACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAATGGGGTTACGACTACACCGGTAACGCTGATAACAACGCGAACTTCGGCAAAGCCGTTCCTGCTGATTTCaa       >  1:378909/1‑139 (MQ=255)                                                                                                                          gACAGCATCTGGCCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAATGGGGTTACGACTACACCGGTAACGCTGATAACAACGCGAACTTCGGCAAAGCCGTTCCTGCTGATTTCAAcggcg  <  2:135232/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                               CAAGTGGACGCCAATCATGAGCACCGTGATGGAAATCGGCTACGACAACGTCGAATCCCAGCGCACCGGCGACAAGAACAATCAGTACAAAATTACCCTCGCACAACAATGGCAGGCTGGCGACAGCATCTGGTCACGCCCGGCTATTCGTGTCTTCGCAACCTACGCCAAGTGGGATGAGAAATGGGGTTACGACTACACCGGTAACGCTGATAACAACGCGAACTTCGGCAAAGCCGTTCCTGCTGATTTCAACGGCG  >  NZ_CP009273/4238915‑4239174

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