Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,649,476	A→G	intergenic (‑78/+128)	rspA ← / ← ynfA	starvation‑sensing protein RspA/YnfA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,649,476	0	A	G	95.8%	72.4 / ‑3.4	24	intergenic (‑78/+128)	rspA/ynfA	starvation‑sensing protein RspA/YnfA family protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (10/13);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-01

CGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTT  >  NZ_CP009273/1649389‑1649554                                                                                        |                                                                               cGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGgaa                                                                              <  2:973358/90‑1 (MQ=255)    tCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGgaaaaa                                                                           >  2:145016/1‑90 (MQ=255)             tCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTg                                                                  >  2:90973/1‑90 (MQ=255)                        cATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACaaaaaa                                                       <  2:721444/90‑1 (MQ=255)                           cGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAgg                                                    <  2:555851/90‑1 (MQ=255)                              ttGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGt                                                 <  2:317402/90‑1 (MQ=255)                              ttGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGt                                                 >  1:60480/1‑90 (MQ=255)                                         cATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCt                                      >  2:392458/1‑90 (MQ=255)                                          aTGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCtg                                     >  1:968697/1‑90 (MQ=255)                                                   aCGCAACTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGtt                            <  1:639689/90‑1 (MQ=255)                                                        aCTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACaa                       <  1:90973/90‑1 (MQ=255)                                                        aCTGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACaa                       <  2:317923/90‑1 (MQ=255)                                                          tGAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACaaaa                     <  1:277473/90‑1 (MQ=255)                                                           gAGCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAg                    >  1:642803/1‑90 (MQ=255)                                                             gCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGAt                  <  2:325517/90‑1 (MQ=255)                                                             gCTACTACCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGAt                  <  2:742457/90‑1 (MQ=255)                                                               tactacCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGa                <  1:426022/90‑1 (MQ=255)                                                                actacCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGaa               >  1:580612/1‑90 (MQ=255)                                                                 ctacCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGaaa              <  1:369967/90‑1 (MQ=255)                                                                  tacCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACaa                       <  1:852512/80‑1 (MQ=255)                                                                  tacCATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACaa                       >  2:852512/1‑80 (MQ=255)                                                                     cATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACaaaaaaa                                                      <  1:918296/46‑1 (MQ=255)                                                                     cATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACaaaaaaa                                                      >  2:918296/1‑46 (MQ=255)                                                                     cATACAAGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGa          >  1:261478/1‑90 (MQ=255)                                                                           aGTATAAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCt    >  2:452688/1‑90 (MQ=255)                                                                             tataAAGATCGGAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCttt  >  2:408116/1‑90 (MQ=255)                                                                                        |                                                                               CGATCTTCATGTTCGATTCCTTGCATCGCTTGTCGTGATGCATGAAATCTACGCAACTGAGCTACTACCATACAAGTATAAAGATCGAAAAAAGCCGGAGTGATCACAAAAAAAGGCGTATATTTGCGCTGTGAATGGTTGACAAAAGATGAAATAGAATACCTTT  >  NZ_CP009273/1649389‑1649554

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