Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,710,430	C→T	A175V (GCG→GTG)	yiaD →	OmpA family lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,710,430	0	C	T	88.2%	168.4 / NA	51	A175V (GCG→GTG)	yiaD	OmpA family lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (22/23);  minor base G (6/0);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 2.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAA  >  NZ_CP009273/3710296‑3710546                                                                                                                                       |                                                                                                                     ccTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGGGATTGGTTATACCGACAGCACCGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGtt                                                                                                                   >  1:71471/1‑138 (MQ=255)  cTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGG‑TAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGAt                                                                                                                <  2:510521/139‑1 (MQ=255)     aCCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGAt                                                                                                                <  1:395387/137‑1 (MQ=255)     aCCGGCGTGGCAATGGTACTGAAAGAGTATCCAAAAACGGCGGTTAACGTGATTGGGTATAACGACAGCACGGGGGGGCACGAACTGAACATGCGTCTCTCCCAGCAACGGGCGGGTTACGTTGCCAGCGGGTTGATCa                                                                                                              >  1:611909/1‑139 (MQ=255)       cGGCGTGGCAATGGTACTGAAAGAGTATTCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCAcc                                                                                                            <  1:675207/139‑1 (MQ=255)       cGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCAcc                                                                                                            <  1:260413/139‑1 (MQ=255)        ggCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCAccc                                                                                                           <  2:445346/139‑1 (MQ=255)         gCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAACACGGGTGGTCACGACCTGAACATGCGTCTCTCCCCGCAACGTGCGGATTCCGTTGGCCGCGTGCTGATCAtcc                                                                                                           >  1:515306/1‑138 (MQ=255)         gCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGGGATTGGGTATACCGACAGCACGGGGGGGCACGACCCGAACATGCGTCTCTCCCAGCAACGGGCGGATTCCGTTGCCAGCGTGTTGATCAcc                                                                                                            >  1:500877/1‑137 (MQ=255)                          aaaGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGc                                                                                          <  1:349181/138‑1 (MQ=255)                           aaGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCg                                                                                        <  2:706640/139‑1 (MQ=255)                              agTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCAt                                                                                     <  1:511737/139‑1 (MQ=255)                                 aTCCGAAAACGGCGGTTAACGTGATTGGTTTTACCCACAGCACGGGGGGGCACGACCTGAAAATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTTCCCGCGGGGTGATCACCCCGGGCGGGGACGCCAGCCGCCTCCg                                                                                  >  2:349181/1‑139 (MQ=255)                                 aTCCGAAAACGGCGGTTAAAGTGATTGGGTATACCGACAGGACGGGGGGGCACGACCTGAAAATGCGTCTCTCCCAGCAACGGGCGGGTTCCGTTGCCAGCGGGTTGAGCACCCAGGGCGTGGGCGCCAGCCGCATCCg                                                                                  >  2:67909/1‑139 (MQ=255)                                  tCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACAAGCGTCTCTCCCAGCAACGGGCGGATTTCGTTGCCAGCGGGTTGATCACCCCGGGCGTGGGCGCCCGCCGCATCCg                                                                                  >  1:210036/1‑138 (MQ=255)                                  tCCGAAAACGGCGGTTAACGTGATTGGGTATACCGACAGCACGGGGGGTCACGACCTGAACATGCGTCTCTCCCCGCAACGTGCGGATTCCGTTGCCAGCGGGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCg                                                                                  >  1:233986/1‑138 (MQ=255)                                       aaaCGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCa                                                                            <  2:825046/139‑1 (MQ=255)                                       aaaCGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCa                                                                            <  2:39335/139‑1 (MQ=255)                                       aaaCGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCa                                                                            <  2:501966/139‑1 (MQ=255)                                            gcggTTAAAGTGATTGGTTATACCGACAGCACGGGGGGGCCAGACCTGAACATGCGTCTTTCCCAGCAAAGGGCGGGTTTCGTTGTCAGCGGGTTGGTCACCCAGGGGGTGGGCGCCAGCCCCATCCCGACTCAGGGcc                                                                       >  2:741323/1‑139 (MQ=255)                                                        aTTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCgtgt                                                                                                                    <  2:317664/82‑1 (MQ=255)                                                        aTTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCgtgt                                                                                                                    >  1:317664/1‑82 (MQ=255)                                                         ttGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTaa                                                          >  1:634980/1‑139 (MQ=255)                                                                   cGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGcc                                                <  2:500877/139‑1 (MQ=255)                                                                     aCAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAg                                              >  1:821985/1‑139 (MQ=255)                                                                     aCAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAg                                              >  1:256442/1‑139 (MQ=255)                                                                      cAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCg                                                                                        >  1:692090/1‑96 (MQ=255)                                                                      cAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCg                                                                                        <  2:692090/96‑1 (MQ=255)                                                                      cAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAg                                              >  2:366680/1‑138 (MQ=255)                                                                       aGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGc                                             >  2:155878/1‑138 (MQ=255)                                                                       aGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGc                                             <  2:413165/138‑1 (MQ=255)                                                                       aGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGc                                             >  1:413165/1‑138 (MQ=255)                                                                       aGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGc                                             >  2:827751/1‑138 (MQ=255)                                                                             ggtggtCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCa                                      <  1:798104/139‑1 (MQ=255)                                                                                ggtCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCg                                   >  2:628990/1‑139 (MQ=255)                                                                                        ccTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCa                                                                            <  2:335023/90‑1 (MQ=255)                                                                                        ccTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCa                                                                            >  1:335023/1‑90 (MQ=255)                                                                                        ccTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGACgg                            >  2:395224/1‑138 (MQ=255)                                                                                            aaCATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAgg                       >  2:237059/1‑139 (MQ=255)                                                                                                  cGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCaa                  >  2:216466/1‑138 (MQ=255)                                                                                                    tctctcCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCg                                                  >  2:533382/1‑104 (MQ=255)                                                                                                    tctctcCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCg                                                  <  1:533382/104‑1 (MQ=255)                                                                                                        tcCCAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAAccg             <  2:821985/137‑1 (MQ=255)                                                                                                           cAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCgtg        <  2:611909/139‑1 (MQ=255)                                                                                                           cAGCAACGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCgtg        <  2:256442/139‑1 (MQ=255)                                                                                                                 cGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGaaa  >  2:306717/1‑139 (MQ=255)                                                                                                                 cGTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGaaa  >  1:565479/1‑139 (MQ=255)                                                                                                                  gTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCaa                                           >  2:57230/1‑97 (MQ=255)                                                                                                                  gTGCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCaa                                           <  1:57230/97‑1 (MQ=255)                                                                                                                    gCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAAccgcc           <  2:213105/127‑1 (MQ=255)                                                                                                                    gCGGATTCCGTTGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAAccgcc           >  1:213105/1‑127 (MQ=255)                                                                                                                              ttGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAAccg             <  2:591456/115‑1 (MQ=255)                                                                                                                              ttGCCAGCGTGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAAccg             >  1:591456/1‑115 (MQ=255)                                                                                                                                       |                                                                                                                     CCTGACCGGCGTGGCAATGGTACTGAAAGAGTATCCGAAAACGGCGGTTAACGTGATTGGTTATACCGACAGCACGGGTGGTCACGACCTGAACATGCGTCTCTCCCAGCAACGTGCGGATTCCGTTGCCAGCGCGTTGATCACCCAGGGCGTGGACGCCAGCCGCATCCGTACTCAGGGCCTTGGCCCGGCTAACCCAATCGCCAGCAACAGCACCGCAGAAGGTAAGGCGCAAAACCGCCGTGTAGAAA  >  NZ_CP009273/3710296‑3710546

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