Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,633,239	(T)8→7	intergenic (‑159/‑232)	uspB ← / → uspA	universal stress protein UspB/universal stress protein UspA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,633,232	0	T	.	100.0%	134.6 / NA	31	intergenic (‑152/‑239)	uspB/uspA	universal stress protein UspB/universal stress protein UspA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (20/11);  total (20/11)

CCTGCCCGCGGAGTTCCACCCCGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACTGAATCTTTA  >  NZ_CP009273/3633097‑3633364                                                                                                                                        |                                                                                                                                     ccTGCCCGCGGAGTTCCACCCCGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGtttt                                                                                                                                   >  2:41998/1‑139 (MQ=255)             gTTCCACCCCGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGt                                                                                                                      <  1:132708/139‑1 (MQ=255)                 cACCCCGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAAc                                                                                                                  <  1:137728/139‑1 (MQ=255)                    cccGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCAc                                                                                                               <  1:482856/139‑1 (MQ=255)                                    ccGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATg                                                                                               >  2:294174/1‑139 (MQ=255)                                              gCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAAttt                                                                                     >  1:299016/1‑139 (MQ=255)                                                              cAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGt                                                                      >  1:424727/1‑138 (MQ=255)                                                                ggTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAg                                                                   >  1:381601/1‑139 (MQ=255)                                                                 gTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAAc                                                                                                                  <  1:110125/91‑1 (MQ=255)                                                                 gTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAAc                                                                                                                  >  2:110125/1‑91 (MQ=255)                                                                 gTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTAt                                                                                                            <  1:188509/97‑1 (MQ=255)                                                                 gTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTAt                                                                                                            <  1:414066/97‑1 (MQ=255)                                                                 gTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTAt                                                                                                            >  2:414066/1‑97 (MQ=255)                                                                 gTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTAt                                                                                                            >  2:188509/1‑97 (MQ=255)                                                                 gTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAg                                                                   >  2:250784/1‑138 (MQ=255)                                                                            aCCGCTTTTGTGGTGACCAACTTACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCg                                                                                            >  1:167760/1‑102 (MQ=255)                                                                            aCCGCTTTTGTGGTGACCAACTTACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCg                                                                                            <  2:167760/102‑1 (MQ=255)                                                                                           aCCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGt                                        <  1:250784/139‑1 (MQ=255)                                                                                                 aTACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCtat                                  <  1:35954/139‑1 (MQ=255)                                                                                                       gCGGCTCTATAGATAGTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTAtaaggt                             >  2:467396/1‑138 (MQ=255)                                                                                                                     agTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGAt                       >  2:446560/1‑130 (MQ=255)                                                                                                                     agTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGAt                       <  1:446560/130‑1 (MQ=255)                                                                                                                     agTGTAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTaa               >  1:143150/1‑138 (MQ=255)                                                                                                                        gtAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCtat                                  <  2:470798/116‑1 (MQ=255)                                                                                                                        gtAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCtat                                  >  1:470798/1‑116 (MQ=255)                                                                                                                        gtAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACt           >  2:60678/1‑139 (MQ=255)                                                                                                                        gtAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACt           >  2:127050/1‑139 (MQ=255)                                                                                                                        gtAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACt           >  1:490065/1‑139 (MQ=255)                                                                                                                        gtAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACt           >  1:151124/1‑139 (MQ=255)                                                                                                                        gtAGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACt           >  1:130944/1‑139 (MQ=255)                                                                                                                          aGGAGATCAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAAATTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACTg          >  2:198490/1‑138 (MQ=255)                                                                                                                                 cAGGTTG‑TTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACTGAATCTTTa  >  2:109352/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                     CCTGCCCGCGGAGTTCCACCCCGGGGCTACCGCTCCCGATACGCTGCCAATCAGTTAACACCAGGTCCTGGAGAAACCGCTTTTGTGGTGACCAACATACGAGCGGCTCTATAGATAGTGTAGGAGATCAGGTTGTTTTTTTTCCAGAAGGTTAACCACTATCAATATATTCATGTCGAAAATTTGTTTATCTAACGAGTAAGCAAGGCGGATTGACGGATCATCCGGGTCGCTATAAGGTAAGGATGGTCTTAACACTGAATCTTTA  >  NZ_CP009273/3633097‑3633364

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