Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	2,982,682	(A)8→7	coding (20/483 nt)	BW25113_RS14815 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	2,982,675	0	A	.	100.0%	146.3 / NA	33	coding (13/483 nt)	BW25113_RS14815	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (19/14);  total (19/14)

ATCCGTGGTTATACTTCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTATAAAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATAACGAAACATGCACTGCAAATTGGGTAATTTTTAAT  >  NZ_CP009273/2982540‑2982797                                                                                                                                        |                                                                                                                           aTCCGTGGTTATACTTCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTATaaaa                                                                                                                         <  2:655843/139‑1 (MQ=255)                 cccTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATtatttatt                                                                                                        >  2:683171/1‑139 (MQ=255)                   ctctTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTatttattt                                                                                                       >  2:243509/1‑138 (MQ=255)                           aaaTCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTa                                                                                              >  2:437526/1‑139 (MQ=255)                                cTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAg                                                                                          <  1:435851/138‑1 (MQ=255)                                          atgatTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATtc                                                                               >  2:545185/1‑139 (MQ=255)                                              ttGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCtt                                                                           <  1:683171/139‑1 (MQ=255)                                                 ttATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCa                                                                         <  1:594052/138‑1 (MQ=255)                                                       cATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTa                                                                  >  1:846662/1‑139 (MQ=255)                                                         tGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTAtt                                                                >  1:481075/1‑139 (MQ=255)                                                           cAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCt                                                              >  1:177708/1‑139 (MQ=255)                                                                aaaaCTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTaca                                                          <  1:243509/138‑1 (MQ=255)                                                                                tGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGata                                         <  2:177708/139‑1 (MQ=255)                                                                                        cAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATAACg                                  >  1:396091/1‑138 (MQ=255)                                                                                               tttAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATAACGAAACATg                           >  2:302166/1‑138 (MQ=255)                                                                                               tttAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATAACGAAACATg                           >  2:790985/1‑138 (MQ=255)                                                                                                 tAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATAACGAAACATGCAc                        >  2:349320/1‑139 (MQ=255)                                                                                                       tGTCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATAACGAAACATGCACTGCaaa                  >  2:125651/1‑139 (MQ=255)                                                                                                         tCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTc                                                                    >  2:314774/1‑87 (MQ=255)                                                                                                         tCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTc                                                                    <  1:314774/87‑1 (MQ=255)                                                                                                         tCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATAACGAAACATGCACTGCaaa                  <  1:349320/137‑1 (MQ=255)                                                                                                         tCTTACCGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATAACGAAACATGCACTGCAAAtt                >  1:439903/1‑139 (MQ=255)                                                                                                              ccGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATAACGAAACATGCACTGCAAATTGGGt            <  1:378223/138‑1 (MQ=255)                                                                                                               cGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTTTTCACTGTc                                                                    >  1:387985/1‑81 (MQ=255)                                                                                                               cGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTc                                                                    <  2:387985/81‑1 (MQ=255)                                                                                                               cGTGAGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATAACGAAACATGCACTGCAAATTGGGTaa          <  2:113056/139‑1 (MQ=255)                                                                                                                   aGGTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATAACGAAACATGCACTGCAAATTGGGTAAtttt      >  2:182885/1‑139 (MQ=255)                                                                                                                     gTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATAGCGAAACATGCACTGCAAATTGGGTAATTTTTa    <  1:545185/139‑1 (MQ=255)                                                                                                                     gTTAATCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATAACGAAACATGCACTGCAAATTGGGTAATTTTTa    <  2:821736/139‑1 (MQ=255)                                                                                                                        aaTCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATAACGAAACATGCACTGCAAATTGGGTAATTTTTAAt  >  2:192993/1‑138 (MQ=255)                                                                                                                         aTCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAg                                                                                          >  1:46413/1‑49 (MQ=255)                                                                                                                         aTCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAg                                                                                          <  2:46413/49‑1 (MQ=255)                                                                                                                         aTCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTc                                                                          >  1:91611/1‑65 (MQ=255)                                                                                                                         aTCGTGATTGATTAT‑AAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTc                                                                          <  2:91611/65‑1 (MQ=255)                                                                                                                                        |                                                                                                                           ATCCGTGGTTATACTTCCCTCTTGCAAAATCCTCACCGGGGATGATTGTTATGGCATGCAATAAAAACTATAACCAACATGAAGTGGCAAACTGTTTAACCCTGTCTTACCGTGAGGTTAATCGTGATTGATTATAAAAAAAACTTATTATTTATTTTAGTTTTTATCAGTGGGTTTATTCTCTTCACTGTCTATTCTTACACTGCAGAAAAAATGATATATAACGAAACATGCACTGCAAATTGGGTAATTTTTAAT  >  NZ_CP009273/2982540‑2982797

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