Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,258,591	(T)6→5	intergenic (+356/+146)	tyrB → / ← BW25113_RS26020	aromatic amino acid transaminase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,258,586	0	T	.	100.0%	171.9 / NA	40	intergenic (+351/+151)	tyrB/BW25113_RS26020	aromatic amino acid transaminase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (23/17);  total (23/17)

ATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTCTTTTTTCCGCTTTGTGTAAGTTTAACCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATC  >  NZ_CP009273/4258497‑4258662                                                                                          |                                                                             atatCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAActct                                                                              >  2:895308/1‑90 (MQ=255) atatCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAActct                                                                              >  2:69312/1‑90 (MQ=255) atatCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAActct                                                                              >  1:306933/1‑90 (MQ=255)     ccTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTCttttt                                                                          >  1:992910/1‑90 (MQ=255)      cTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTCtttttc                                                                         >  1:955483/1‑89 (MQ=255)        cAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTCtttttccg                                                                       >  1:160084/1‑87 (MQ=255)        cAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTCtttttccg                                                                       >  2:628444/1‑87 (MQ=255)        cAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTCtttttccg                                                                       >  2:364487/1‑87 (MQ=255)             atatGCTTTGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTtg                                                                 >  2:169956/1‑90 (MQ=255)             atatGCTTTGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTtg                                                                 >  1:308611/1‑90 (MQ=255)               atGCTTTGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTtgtg                                                               >  2:976673/1‑90 (MQ=255)                tGCTTTGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTgtgt                                                              <  2:832333/90‑1 (MQ=255)                 gCTTTGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTa                                                             <  1:976673/90‑1 (MQ=255)                    ttGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCtt                                                                   >  1:37477/1‑81 (MQ=255)                    ttGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCtt                                                                   <  2:652156/81‑1 (MQ=255)                    ttGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCtt                                                                   <  2:37477/81‑1 (MQ=255)                    ttGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCtt                                                                   >  1:652156/1‑81 (MQ=255)                    ttGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGt                                                          >  2:173393/1‑90 (MQ=255)                    ttGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGt                                                          >  1:236207/1‑90 (MQ=255)                    ttGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGt                                                          >  1:288037/1‑90 (MQ=255)                    ttGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGt                                                          >  1:435673/1‑90 (MQ=255)                      gCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGttt                                                        <  1:895308/90‑1 (MQ=255)                          tATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAAcc                                                    >  2:774216/1‑90 (MQ=255)                                cTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGGGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTatat                                              >  1:580093/1‑90 (MQ=255)                                cTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTatat                                              >  2:692271/1‑90 (MQ=255)                                cTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTatat                                              >  1:452580/1‑90 (MQ=255)                                cTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTatat                                              <  1:122919/90‑1 (MQ=255)                                cTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTatat                                              >  1:40023/1‑90 (MQ=255)                                   tgACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTATATTGa                                           <  2:777417/90‑1 (MQ=255)                                        ttCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTATATTGAtt                                         <  2:355873/87‑1 (MQ=255)                                        ttCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTATATTGAtt                                         >  1:355873/1‑87 (MQ=255)                                          cTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTATATTGATTACCaa                                    <  1:973383/90‑1 (MQ=255)                                                ttCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTATATTGATTACCAATATCaa                              <  1:169956/90‑1 (MQ=255)                                                  cAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTATATTGATTACCAATATCAAAt                            >  1:715830/1‑90 (MQ=255)                                                         aGTATCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTATATTGATTACCAATATCAAATAAAAAgg                     <  1:198236/90‑1 (MQ=255)                                                            aTCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTgtgt                                                              >  2:920634/1‑46 (MQ=39)                                                            aTCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTgtgt                                                              <  1:920634/46‑1 (MQ=39)                                                            aTCGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTATATTGATTACCAATATCAAATAAAAAGGTcc                  >  1:443878/1‑90 (MQ=255)                                                              cGTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTATATTGATTACCAATATCAAATAAAAAGGTCCTg                >  2:522625/1‑90 (MQ=255)                                                               gTTAAATTCTAATGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGt               >  2:794124/1‑90 (MQ=255)                                                                         aaTGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGaata     >  1:1029830/1‑90 (MQ=255)                                                                         aaTGTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGaata     >  2:789985/1‑90 (MQ=255)                                                                            gTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATc  <  2:443878/90‑1 (MQ=255)                                                                            gTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATc  <  1:364487/90‑1 (MQ=255)                                                                            gTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATc  <  1:628444/90‑1 (MQ=255)                                                                            gTTAAAGAGAACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATc  <  1:692271/90‑1 (MQ=255)                                                                                    gaACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATc  >  2:258623/1‑82 (MQ=255)                                                                                    gaACTC‑TTTTTCCGCTTTGTGTAAGTTTAACCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATc  <  1:258623/82‑1 (MQ=255)                                                                                          |                                                                             ATATCCTCAAGAATATGCTTTGCCTTATTGCCTGTGACTTTCTGAGATTCAATTATAGTATCGTTAAATTCTAATGTTAAAGAGAACTCTTTTTTCCGCTTTGTGTAAGTTTAACCTTATATTGATTACCAATATCAAATAAAAAGGTCCTGTAATGAGAATAATC  >  NZ_CP009273/4258497‑4258662

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