Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,179,271	G→A	A167T (GCC→ACC)	ygiL →	fimbrial‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,179,271	0	G	A	96.9%	105.6 / NA	33	A167T (GCC→ACC)	ygiL	fimbrial‑like protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (12/19);  minor base C (0/1);  total (13/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCCGGTAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCA  >  NZ_CP009273/3179188‑3179351                                                                                    |                                                                                 cTCCGGTAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCAccc                                                                            >  2:385255/1‑90 (MQ=255)       tagtaATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCa                                                                      <  1:55085/90‑1 (MQ=255)       tagtaATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCa                                                                      <  1:561309/90‑1 (MQ=255)        agtaATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCGCGGGCaa                                                                     >  1:431672/1‑90 (MQ=255)         gtaATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAAc                                                                    <  2:261511/90‑1 (MQ=255)             ttACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGttt                                                                <  1:285493/90‑1 (MQ=255)             ttACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGttt                                                                <  1:425020/90‑1 (MQ=255)                   tCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTc                                                          <  2:48995/90‑1 (MQ=255)                    caacaaTACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCa                                                         <  1:385255/90‑1 (MQ=255)                       caaTACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGc                                                      <  1:9399/90‑1 (MQ=255)                         aTACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATCCCACCCCGGGCAACGTTTACGCTCAGGCTg                                                    <  1:411000/90‑1 (MQ=255)                         aTACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTg                                                    <  1:185371/90‑1 (MQ=255)                         aTACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTg                                                    <  2:412994/90‑1 (MQ=255)                            ccAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACt                                                 <  1:604225/90‑1 (MQ=255)                            ccAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACt                                                 >  1:53287/1‑90 (MQ=255)                              aaCAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTAc                                               >  2:532586/1‑90 (MQ=255)                                cAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACaa                                             >  1:484777/1‑90 (MQ=255)                                     aTCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTg                                        >  1:586013/1‑90 (MQ=255)                                     aTCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTg                                        >  1:417591/1‑90 (MQ=255)                                            aCTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATTCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGa                                 >  2:520449/1‑90 (MQ=255)                                            aCTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGa                                 >  1:493888/1‑90 (MQ=255)                                            aCTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGa                                 >  2:152655/1‑90 (MQ=255)                                                  tttCAAAGCTTGTATGGAATCCCCAAGCAAAGGTACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATc                           <  2:110395/90‑1 (MQ=255)                                                  tttCAAAGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATc                           <  2:390026/90‑1 (MQ=255)                                                       aaGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCAtt                      <  1:532586/90‑1 (MQ=255)                                                       aaGCTCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCAtt                      >  1:371861/1‑90 (MQ=255)                                                           tCGTATGGAATCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGc                  <  1:72664/90‑1 (MQ=255)                                                                    aTCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGc                                                      >  2:103162/1‑45 (MQ=255)                                                                    aTCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGc                                                      <  1:103162/45‑1 (MQ=255)                                                                     tCCCCAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGcc        <  2:553366/90‑1 (MQ=255)                                                                      ccccAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTg                                                    <  1:172521/45‑1 (MQ=255)                                                                      ccccAAGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTg                                                    >  2:172521/1‑45 (MQ=255)                                                                           aGCAAAGATACCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCACCCGGCca  <  2:605937/90‑1 (MQ=255)                                                                                    |                                                                                 CTCCGGTAGTAATTACGTTCAACAATACCAACAGCTATCAGGAACTGAATTTCAAAGCTCGTATGGAATCCCCAAGCAAAGATGCCACCCCGGGCAACGTTTACGCTCAGGCTGACTACAAGATTGCTTACGAGTAATCGCATTTGGCCGGGTTCGCCCGGCCA  >  NZ_CP009273/3179188‑3179351

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