Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	524,350	T→C	F85F (TTT→TTC)	BW25113_RS24960 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	524,350	0	T	C	100.0%	104.6 / NA	30	F85F (TTT→TTC)	BW25113_RS24960	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (19/11);  total (19/11)

ACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAG  >  NZ_CP009273/524229‑524475                                                                                                                          |                                                                                                                              aCTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAatgatgat                                                                                                              >  2:193848/1‑139 (MQ=255)             tACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGtt                                                                                                   >  1:142737/1‑138 (MQ=255)             tACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGtt                                                                                                   >  2:255370/1‑138 (MQ=255)             tACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGtt                                                                                                   >  2:170737/1‑138 (MQ=255)             tACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGtt                                                                                                   >  1:278751/1‑138 (MQ=255)                  gAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTaa                                                                                                                         <  2:26138/111‑1 (MQ=255)                  gAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTaa                                                                                                                         >  1:26138/1‑111 (MQ=255)                          tgGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGa                                                                                     >  1:306353/1‑139 (MQ=255)                                 tttACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAg                                                                              >  1:316936/1‑139 (MQ=255)                                         gCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAgg                                                                       >  2:352080/1‑138 (MQ=255)                                         gCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAgg                                                                       >  1:19936/1‑138 (MQ=255)                                                                  gCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCaaa                                             <  1:51971/139‑1 (MQ=255)                                                                    ggTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGa                                           <  1:193848/139‑1 (MQ=255)                                                                         atatatCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGAt                                      >  1:61798/1‑139 (MQ=255)                                                                             atCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGa                                           <  2:214116/130‑1 (MQ=255)                                                                             atCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGa                                           >  1:214116/1‑130 (MQ=255)                                                                             atCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATc                                     >  1:260543/1‑136 (MQ=255)                                                                             atCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATc                                     >  2:402746/1‑136 (MQ=255)                                                                                    taCACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACaa                            <  2:278751/138‑1 (MQ=255)                                                                                           aTGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAgtg                    <  2:19936/139‑1 (MQ=255)                                                                                           aTGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAgtg                    >  2:59862/1‑139 (MQ=255)                                                                                            tGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAgtgt                   >  2:213000/1‑139 (MQ=255)                                                                                                    atgatTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTaaa           <  1:237060/139‑1 (MQ=255)                                                                                                     tgatTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTaaa           >  1:207859/1‑138 (MQ=255)                                                                                                       attattTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGa                                                                                     >  1:114085/1‑62 (MQ=255)                                                                                                       attattTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGa                                                                                     <  2:114085/62‑1 (MQ=255)                                                                                                       attattTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGa                                                                                     <  2:173882/62‑1 (MQ=255)                                                                                                       attattTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGa                                                                                     >  1:173882/1‑62 (MQ=255)                                                                                                           tttttttCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTa    <  1:402746/139‑1 (MQ=255)                                                                                                              ttttCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAg  <  2:61798/138‑1 (MQ=255)                                                                                                                          |                                                                                                                              ACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAG  >  NZ_CP009273/524229‑524475

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: ‑