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.0 / NA	28	F85F (TTT→TTC)	BW25113_RS24960	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (21/7);  total (21/7)

TGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATG  >  NZ_CP009273/524267‑524437                                                                                    |                                                                                        tGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAAc                                                                                  <  2:156861/91‑1 (MQ=255)            tttttttCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAatgatgatg                                                                       <  1:123134/91‑1 (MQ=255)            tttttttCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAatgatgatg                                                                       <  2:5279/91‑1 (MQ=255)                           aGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGa                                                        <  1:112370/91‑1 (MQ=255)                                tGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACa                                                   >  1:946626/1‑91 (MQ=255)                                 ggATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAatgatga                                                                         >  2:340746/1‑68 (MQ=255)                                 ggATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAatgatga                                                                         <  1:340746/68‑1 (MQ=255)                                  gATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACAtt                                                 >  2:797609/1‑91 (MQ=255)                                     atatCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGAc                                              <  2:709553/91‑1 (MQ=255)                                      tatCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACa                                             >  1:378550/1‑91 (MQ=255)                                      tatCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACa                                             >  1:162314/1‑91 (MQ=255)                                      tatCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACa                                             >  1:1116883/1‑91 (MQ=255)                                             ataCACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAggg                                      >  1:55956/1‑91 (MQ=255)                                             ataCACAAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAggg                                      >  1:30205/1‑91 (MQ=255)                                                  caAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAgg                                 >  2:108187/1‑91 (MQ=255)                                                  caAATGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAgg                                 >  2:883014/1‑91 (MQ=255)                                                   aaaTGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAggg                                >  2:75081/1‑91 (MQ=255)                                                    aaTGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGa                               >  2:890160/1‑91 (MQ=255)                                                     aTGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGaa                              >  1:92305/1‑91 (MQ=255)                                                     aTGTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGaa                              >  2:810719/1‑91 (MQ=255)                                                       gTAGAAGATGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACa                            >  2:517036/1‑91 (MQ=255)                                                         agaagaTGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAg                          >  1:766326/1‑91 (MQ=255)                                                            agaTGATTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTgg                       >  2:91030/1‑91 (MQ=255)                                                               tgatTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGa                    >  1:63799/1‑91 (MQ=255)                                                                gatTATTTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGaa                   >  1:428487/1‑91 (MQ=255)                                                                 attattTTTTTCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGaaa                  >  1:380579/1‑91 (MQ=255)                                                                       tttttCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTa            >  1:649954/1‑91 (MQ=255)                                                                           tCTTATTTTCCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCaa        <  1:959557/91‑1 (MQ=255)                                                                                 tttcCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATg  <  2:1049766/91‑1 (MQ=255)                                                                                    |                                                                                        TGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATG  >  NZ_CP009273/524267‑524437

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