Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	882,114	C→G	A228P (GCC→CCC)	ybjJ ←	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	882,114	0	C	G	98.3%	221.4 / ‑3.2	60	A228P (GCC→CCC)	ybjJ	MFS transporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base G (30/29);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.62e-01

AAATCAGCGAGCCGGAAGTAGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGCCATCGCCAGCACCACAACACCTATCAGCAGCAACTGGATATCGCGATAAAAAGGTACGCCTTTTTCGCCATGCTGGGTGCCATCG  >  NZ_CP009273/882027‑882199                                                                                        |                                                                                      aaaTCAGCGAGCCGGAAGTAGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCa                                                                                     <  2:792080/90‑1 (MQ=255)  aaTCAGCGAGCCGGAAGTAGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCAt                                                                                    <  2:2150737/90‑1 (MQ=255)     cAGCGAGCCGGAAGTAGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGc                                                                                 >  2:469607/1‑90 (MQ=255)       gCGAGCCGGAAGTAGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCa                                                                               >  1:670703/1‑90 (MQ=255)       gCGAGCCGGAAGTAGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCa                                                                               <  1:878707/90‑1 (MQ=255)       gCGAGCCGGAAGTAGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCa                                                                               >  1:2407491/1‑90 (MQ=255)       gCGAGCCGGAAGTAGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCa                                                                               >  1:1510207/1‑90 (MQ=255)       gCGAGCCGGAAGTAGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCa                                                                               >  2:2525886/1‑90 (MQ=255)         gAGCCGGAAGTAGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGc                                                                             <  1:2906082/90‑1 (MQ=255)         gAGCCGGAAGTAGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGc                                                                             <  1:810278/90‑1 (MQ=255)         gAGCCGGAAGTAGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGc                                                                             <  2:1876877/90‑1 (MQ=255)         gAGCCGGAAGTAGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGc                                                                             <  1:1472032/90‑1 (MQ=255)         gAGCCGGAAGTAGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGc                                                                             <  1:1260792/90‑1 (MQ=255)         gAGCAGAAAGTCCGGCTAAAACCGTGTCCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGc                                                                             <  1:2752064/90‑1 (MQ=255)            ccGGAAGTAGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGcacc                                                                          <  1:1155654/90‑1 (MQ=255)                 aGTAGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCacaac                                                                     <  1:2304601/90‑1 (MQ=255)                 aGTAGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCacaac                                                                     <  1:2732492/90‑1 (MQ=255)                 aGTAGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCacaac                                                                     <  1:1984893/90‑1 (MQ=255)                 aGTAGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCacaac                                                                     <  1:1035593/90‑1 (MQ=255)                    aGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGca                                                                            >  1:967586/1‑80 (MQ=255)                    aGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGca                                                                            <  2:967586/80‑1 (MQ=255)                      ggCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTa                                                                <  1:1562362/90‑1 (MQ=255)                            aaCCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTAt                                                               >  1:1051954/1‑85 (MQ=255)                            aaCCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTAt                                                               <  2:1051954/85‑1 (MQ=255)                            aaCCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATcagca                                                          <  1:1695047/90‑1 (MQ=255)                            aaCCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATcagca                                                          <  2:383724/90‑1 (MQ=255)                              ccGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCagcagc                                                        <  2:720869/90‑1 (MQ=255)                               cGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAgcagca                                                       >  2:908000/1‑90 (MQ=255)                               cGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAgcagca                                                       >  1:1724549/1‑90 (MQ=255)                                gTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCGAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCaa                                                      >  2:1020952/1‑90 (MQ=255)                                gTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATcagca                                                          >  1:304439/1‑86 (MQ=255)                                gTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATcagca                                                          <  2:304439/86‑1 (MQ=255)                                gTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCaa                                                      >  1:2242171/1‑90 (MQ=255)                                gTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCaa                                                      >  2:566357/1‑90 (MQ=255)                                gTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCaa                                                      >  1:155668/1‑90 (MQ=255)                                 tGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCAAc                                                     >  2:932998/1‑90 (MQ=255)                                  gACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCAACt                                                    >  1:2733594/1‑90 (MQ=255)                                           ccATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCAACTGGATATcgc                                           <  1:2636877/90‑1 (MQ=255)                                           ccATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCAACTGGATATcgc                                           <  1:1549603/90‑1 (MQ=255)                                           ccATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCAACTGGATATcgc                                           <  1:1986978/90‑1 (MQ=255)                                                     gggTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGc                                                                             >  2:1879534/1‑46 (MQ=255)                                                     gggTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGc                                                                             <  1:1879534/46‑1 (MQ=255)                                                      ggTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATcagca                                                          >  1:2652445/1‑64 (MQ=255)                                                      ggTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATcagca                                                          <  2:2652445/64‑1 (MQ=255)                                                        tAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCAACTGGATATCGCGATAAAAAGGTAc                              >  2:1853812/1‑90 (MQ=255)                                                        tAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCAACTGGATATCGCGATAAAAAGGTAc                              >  2:241933/1‑90 (MQ=255)                                                        tAACCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCAACTGGATATCGCGATAAAAAGGTAc                              >  2:2643766/1‑90 (MQ=255)                                                         aaCCAGTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCAACTGGATATCGCGATAAAAAGGTACg                             >  1:510621/1‑90 (MQ=255)                                                              gTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCAACTGGATATCGCGATAAAAAGGTACGCCttt                        >  1:998144/1‑90 (MQ=255)                                                              gTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCAACTGGATATCGCGATAAAAAGGTACGCCttt                        <  2:2422247/90‑1 (MQ=255)                                                              gTCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCAACTGGATATCGCGATAAAAAGGTACGCCttt                        >  1:1466306/1‑90 (MQ=255)                                                               tCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCAACTGGATATCGCGATAAAAAGGTACGCCtttt                       >  2:731548/1‑90 (MQ=255)                                                               tCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCAACTGGATATCGCGATAAAAAGGTACGCCtttt                       >  1:2091086/1‑90 (MQ=255)                                                               tCGTTGGCAGAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCAACTGGATATCGCGATAAAAAGGTACGCCtttt                       >  2:941754/1‑90 (MQ=255)                                                                   tGGCAGAACCTTCGGCAAAGGGCATCGCTAGCACCACAACACCTATCAGCAGCAACTGGATATCGCGATAAAAAGGTACGCCTTTTTCGc                   <  1:257588/90‑1 (MQ=255)                                                                        gAACCTTCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCAACTGGATATCGCGATAAAAAGGTACGCCTTTTTCGCCATGc              >  2:12451/1‑90 (MQ=255)                                                                          acattCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCAACTGGATATCGCGATAAAAAGGTACGCCTTTTTCGCCATGCTg            <  2:1705726/87‑1 (MQ=255)                                                                              tCGGCAAAGGGCATCGCCAGCACCACAACACCTATCAGCAGCAACTGGATATCGCGATAAAAAGGTACGCCTTTTTCGCCATGCTGGGTg        >  1:2890424/1‑90 (MQ=255)                                                                                    aaGGGCATCGCCAGCACCACAACACCTATCAGCAGCAACTGGATATCGCGATAAAAAGGTACGCCTTTTTCGCCATGCTGGGTGCCATCg  >  1:902056/1‑90 (MQ=255)                                                                                    aaGGCCATCGCCAGCACCACAACACCTATCAGCAGCAACTGGATATCGCGATAAAAAGGTACGCCTTTTTCGCCATGCTGGGTGCCATCg  >  2:2329538/1‑90 (MQ=255)                                                                                        |                                                                                      AAATCAGCGAGCCGGAAGTAGGGCTAAAACCGTGACCGTCAACCATTAATAAGGGTAACCAGTCGTTGGCAGAACCTTCGGCAAAGGCCATCGCCAGCACCACAACACCTATCAGCAGCAACTGGATATCGCGATAAAAAGGTACGCCTTTTTCGCCATGCTGGGTGCCATCG  >  NZ_CP009273/882027‑882199

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