Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,290,243	0	T	G	49.4%	97.4 / 21.6	86	L87F (TTA→TTC)	ccmB	heme exporter protein CcmB
Reads supporting (aligned to +/- strand):  ref base T (5/38);  new base G (42/0);  total (48/38)
Fisher's exact test for biased strand distribution p-value = 5.16e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATTGCTACCAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAAGCAGCGC  >  NZ_CP009273/2290154‑2290332                                                                                          |                                                                                          aTTGCTACCAGTGGCGAAAGGATGAGAAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGt                                                                                           >  1:272429/1‑90 (MQ=255)   tGCTACCAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTaa                                                                                         <  2:2163627/90‑1 (MQ=255)        ccAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGca                                                                                    <  2:2646285/90‑1 (MQ=255)         cAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGcaa                                                                                   <  2:1057634/90‑1 (MQ=255)         cAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGcaa                                                                                   <  2:125128/90‑1 (MQ=255)         cAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGcaa                                                                                   <  1:1376632/90‑1 (MQ=255)         cAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGcaa                                                                                   <  1:2346740/90‑1 (MQ=255)         cAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGcaa                                                                                   >  1:1172982/1‑90 (MQ=255)          aGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGcaac                                                                                  <  1:2765292/90‑1 (MQ=255)          aGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGcaac                                                                                  <  2:2186313/90‑1 (MQ=255)           gTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGcaaca                                                                                 <  1:262856/90‑1 (MQ=255)            tGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGGAAGGGGAAAGGcaacaa                                                                                >  2:1549667/1‑90 (MQ=255)            tGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGcaacaa                                                                                >  1:2487559/1‑90 (MQ=255)            tGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAAAGGcaacaa                                                                                >  1:133822/1‑90 (MQ=255)             ggCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGGACGGGaacaac                                                                               >  2:454989/1‑90 (MQ=255)              gCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAacaaca                                                                              <  1:1786719/90‑1 (MQ=255)              gCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAacaaca                                                                              >  1:3322084/1‑90 (MQ=255)              gCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAacaaca                                                                              >  1:1838724/1‑90 (MQ=255)                   aGGATGAGTAACGGCAGGCCGGTTTCCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAAt                                                                         <  2:69812/90‑1 (MQ=255)                   aGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAAt                                                                         <  2:392744/90‑1 (MQ=255)                   aGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAAt                                                                         <  1:3244904/90‑1 (MQ=255)                   aGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAAt                                                                         <  2:3168470/90‑1 (MQ=255)                   aGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAAt                                                                         <  1:1562177/90‑1 (MQ=255)                           tAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAg                                                                 <  2:410251/90‑1 (MQ=255)                           tAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAg                                                                 <  2:3144700/90‑1 (MQ=255)                           tAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAg                                                                 <  1:1365161/90‑1 (MQ=255)                           tAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAAGGGAAACGCAACAACATCAAATGttgtag                                                                 >  2:1242747/1‑86 (MQ=255)                              cGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTTAAGACt                                                              >  1:1231224/1‑90 (MQ=255)                              cGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAAAGGGAAAAAAAACAAAAGTTCAAGAcc                                                              >  1:3230831/1‑89 (MQ=255)                               ggCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTg                                                             >  1:1593161/1‑90 (MQ=255)                               ggCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACAACAATTTTTCAAAACTg                                                             >  2:2159193/1‑90 (MQ=255)                                  aGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAAAGGCAACAACATCAATTTTTTAAAACTGCCg                                                          >  1:2508879/1‑90 (MQ=255)                                   gACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGGAAGGGGAAAGGCCACAACAACAAATTTTCAAAACTGCCgg                                                         >  2:633427/1‑89 (MQ=255)                                   gACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGc                                                            <  2:2320478/87‑1 (MQ=255)                                   gACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGc                                                            >  1:2320478/1‑87 (MQ=255)                                   gACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGt                                                         >  1:1284542/1‑90 (MQ=255)                                   gACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAAAATCAAATTTTCAAAACTGCCGt                                                         >  1:2182381/1‑90 (MQ=255)                                   gACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACAGCAACAACATCAAATGTTCAAGAATGGCgg                                                         >  2:2926414/1‑89 (MQ=255)                                   gACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAAAGGGCACAACACCAATTTTTCCAGACTGCCgg                                                         >  2:1863570/1‑89 (MQ=255)                                    aCCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGGGAGGGGAAACGCAACAACAACAATTTTTTAAGAATGCCGTc                                                        >  2:1549863/1‑90 (MQ=255)                                    aCCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAAAGGCAAAAACATCAATTTTTTAAGACTGCCggc                                                        >  2:1793771/1‑88 (MQ=255)                                     caggTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTcc                                                       <  1:347334/88‑1 (MQ=255)                                       ggTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTg                                                     >  1:1838728/1‑90 (MQ=255)                                        gTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGc                                                    <  2:796989/90‑1 (MQ=255)                                        gTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTTTTCAAGACTGCCGTCCTGc                                                    >  1:3290531/1‑90 (MQ=255)                                          tACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCaa                                                  >  1:3008542/1‑90 (MQ=255)                                             cATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTc                                                                    >  1:1847219/1‑69 (MQ=255)                                             cATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTc                                                                    <  2:1847219/69‑1 (MQ=255)                                               tCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGCAACGGCAACAACATCAATTGTTCAAAACTGCCGGCCTGCAAAACgg                                             >  1:1303314/1‑89 (MQ=255)                                                ccAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGGCCTGCAAAACGTc                                            >  1:3172618/1‑90 (MQ=255)                                                 cAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAAGGGAAAGGCAAAAAAAGCAAATTTTCAAGACTTGCGTCCCGCAAAACGGCa                                           >  1:2312086/1‑90 (MQ=255)                                                  aGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAAt                                                                         >  2:3094770/1‑59 (MQ=255)                                                  aGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAAt                                                                         >  1:2957534/1‑59 (MQ=255)                                                  aGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAAt                                                                         <  1:3094770/59‑1 (MQ=255)                                                  aGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAAt                                                                         <  2:2957534/59‑1 (MQ=255)                                                  aGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCAc                                          <  2:2741614/90‑1 (MQ=255)                                                  aGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTTTTCAAGACTGCCGTCCTGCAAAACGTCAc                                          >  1:926692/1‑90 (MQ=255)                                                  aGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAAACTGCCGGCCTTCAAAACGTctg                                          >  2:347334/1‑88 (MQ=255)                                                   gTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACAACAATTTTTCAAAAATGCCGCCCTGCAAAACCTCACg                                         >  2:1015712/1‑90 (MQ=255)                                                   gTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACAACAATTGTTCAAGACTCCCGTCCTGCAAACCGGCACg                                         >  1:2602976/1‑90 (MQ=255)                                                       gcCATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTTTTCAAAACTGCCGGCCTGCAAAACGTCACGGAAc                                     >  2:2412793/1‑90 (MQ=255)                                                        ccATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACa                                    <  2:2692791/90‑1 (MQ=255)                                                        ccATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTTAAGACTGCCGTCCTGCAAATCGTCACGGAACa                                    >  1:3276597/1‑90 (MQ=255)                                                         cATCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGGCCTGCAAATCGTCACGGGACAg                                   >  1:1513378/1‑90 (MQ=255)                                                          aTCACCTTCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGGCACGGAACAGt                                  >  1:690630/1‑90 (MQ=255)                                                                ttCGCCAGCACAACGGCGGGCAAGGGGGACGGCAACAACAACAATTGTTTAAAACAGCCGGCCCGCAAAACGTCACGGAACAATCGGTcc                            >  2:1531182/1‑90 (MQ=255)                                                                 tCGCCAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCa                           >  1:3291043/1‑90 (MQ=255)                                                                    ccAGCACAACGGCGGGGAAGGGGAACGGCAAAAAAATCAATTTTTCAAGACTGCCCGCCCGCAAAACGGCACCGGACAGGCGTTCCAgcg                        >  2:1930538/1‑90 (MQ=255)                                                                    ccAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAgcg                        <  2:1284542/90‑1 (MQ=255)                                                                    ccAGCACAACGGCGGGCAAGGGGAACGGCAACAACATCAATTTTTCAAAACTGCCGTCCTGCAAAACGTCACGGGACAATCGTTCCAgcg                        >  2:188556/1‑90 (MQ=255)                                                                        cacaACGGCGGGCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGGCGTTCCAGCGCCAg                    >  1:1153668/1‑90 (MQ=255)                                                                           aaCGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCaa                 <  2:1303314/90‑1 (MQ=255)                                                                           aaCGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCaa                 <  1:1454029/90‑1 (MQ=255)                                                                           aaCGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCaa                 <  1:1549863/90‑1 (MQ=255)                                                                           aaCGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCaa                 <  1:1806526/90‑1 (MQ=255)                                                                           aaCGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCaa                 <  1:1863570/90‑1 (MQ=255)                                                                           aaCGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCaa                 <  1:454989/90‑1 (MQ=255)                                                                           aaCGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCaa                 <  1:896517/90‑1 (MQ=255)                                                                                cggGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATg            <  2:994987/90‑1 (MQ=255)                                                                                 gggCAAGGGTAACGGCAACAACATCAATTGTTCAAg                                                                 <  2:1855638/36‑1 (MQ=255)                                                                                 gggCAAGGGGAACGGCAACAACATCAATTTTTTAAGACTGCCGTCCTGCAAATCGTCACGGAACAGGCGTTCCAGCGGCAGCAAGGGTGa           >  1:2380963/1‑90 (MQ=255)                                                                                 gggCAAGGGGAACGGCAACAACATCAATTGTTCCAGACTGCCGCCCTGTCAATTGTCACGGAACAGTCGTTCCAGCGCCCGAAAAGATGa           >  1:2430828/1‑90 (MQ=255)                                                                                 gggCAAGGGGAACGGCAACAACATCAATTGTTCAAg                                                                 >  1:1855638/1‑36 (MQ=255)                                                                                 gggCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAAACGTCACGGAACAGGCGTTCCAGCGCCAGCAAGGATGa           >  2:1132228/1‑90 (MQ=255)                                                                                 gggCAAGGGGAACGGCAACAACATCAATTGTTCAAGACTGCCGGCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGGTGa           >  2:3389752/1‑90 (MQ=255)                                                                                   gCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGaaa         <  1:1015712/90‑1 (MQ=255)                                                                                        ggTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAagcagc    <  2:487598/90‑1 (MQ=255)                                                                                          tAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAAgcagcgc  <  1:1276485/90‑1 (MQ=255)                                                                                          |                                                                                          ATTGCTACCAGTGGCGAAAGGATGAGTAACGGCAGACCGGTTACCATCCAGTGCGCCATCACCTTCGCCAGCACAACGGCGGGCAAGGGTAACGGCAACAACATCAATTGTTCAAGACTGCCGTCCTGCAAATCGTCACGGAACAGTCGTTCCAGCGCCAGCAAGGATGAAAGCAGCGC  >  NZ_CP009273/2290154‑2290332

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