Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,976,582	T→G	Y147S (TAC→TCC)	araH ←	arabinose ABC transporter permease AraH

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,976,582	0	T	G	100.0%	124.7 / NA	33	Y147S (TAC→TCC)	araH	arabinose ABC transporter permease AraH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (15/18);  total (15/18)

AGCAGCAAACCAAAGATAATCAGACACGCGACGGTGAGCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGTACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGAACGCCCAGCAACAACCCCGCTGCCACGCCAA  >  NZ_CP009273/1976450‑1976711                                                                                                                                     |                                                                                                                                  agcagcAAACCAAAGATAATCAGACACGCGACGGTGAGCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGcc                                                                                                                              <  2:75544/138‑1 (MQ=255) agcagcAAACCAAAGATAATCAGACACGCGACGGTGAGCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGcc                                                                                                                              <  2:616655/138‑1 (MQ=255) agcagcAAACCAAAGATAATCAGACACGCGACGGTGAGCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGcc                                                                                                                              <  2:569100/138‑1 (MQ=255)   cagcaAACCAAAGATAATCAGACACGCGACGGTGAGCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGa                                                                                                                           >  1:342290/1‑139 (MQ=255)      caAACCAAAGATAATCAGACACGCGACGGTGAGCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCt                                                                                                                        >  1:546253/1‑139 (MQ=255)                   aTCAGACACGCGACGGTGAGCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCa                                                                                                           >  2:346746/1‑139 (MQ=255)                            gcgACGGTGAGCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCa                                                                                                  <  2:626544/139‑1 (MQ=255)                            gcgACGGTGAGCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCa                                                                                                  <  1:346746/139‑1 (MQ=255)                            gcgACGGTGAGCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCa                                                                                                  <  2:497756/139‑1 (MQ=255)                                      gCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACTGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGAt                                                                                        >  1:242856/1‑139 (MQ=255)                                      gCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGAt                                                                                        >  2:197096/1‑139 (MQ=255)                                                         aCCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGt                                                                     <  2:567460/139‑1 (MQ=255)                                                               ccAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGc                                                                >  2:249206/1‑138 (MQ=255)                                                                           aCCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTg                                                                                                     >  1:621114/1‑89 (MQ=255)                                                                           aCCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTg                                                                                                     <  2:621114/89‑1 (MQ=255)                                                                           aCCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCa                                                   <  1:207731/139‑1 (MQ=255)                                                                                gggCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGAc                                              <  2:342290/139‑1 (MQ=255)                                                                                 ggCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGAc                                              >  1:35768/1‑138 (MQ=255)                                                                                 ggCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGAc                                              >  1:203391/1‑138 (MQ=255)                                                                                       gaagCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCca                                       <  2:203391/139‑1 (MQ=255)                                                                                            tttCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCcaca                                     >  2:613012/1‑136 (MQ=255)                                                                                               cATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGaa                                <  1:613012/138‑1 (MQ=255)                                                                                               cATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGaa                                <  1:379015/138‑1 (MQ=255)                                                                                                          ccgaccgCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGAACGCCCAGcaac                     <  1:249206/138‑1 (MQ=255)                                                                                                             accgCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGt                                                                     >  1:81792/1‑87 (MQ=255)                                                                                                             accgCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGt                                                                     <  2:81792/87‑1 (MQ=255)                                                                                                             accgCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACa                                                         <  2:228562/99‑1 (MQ=255)                                                                                                             accgCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACa                                                         >  1:228562/1‑99 (MQ=255)                                                                                                             accgCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGAACGCCCAGCAACAAcc                 >  2:42846/1‑139 (MQ=255)                                                                                                               cgCTTTACCGTCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGAACGCCCAGCAACAAcccc               <  2:242856/139‑1 (MQ=255)                                                                                                                         tCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGAt                                                                                        <  2:389229/56‑1 (MQ=255)                                                                                                                         tCTGAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGAt                                                                                        >  1:389229/1‑56 (MQ=255)                                                                                                                            gAAATGATGGACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGAACGCCCAGCAACAACCCCGCTGCCACGCCaa  >  1:387291/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                  AGCAGCAAACCAAAGATAATCAGACACGCGACGGTGAGCCAGATTGGCGCAGGCAGACCGAACCAGTTGGCGTAACCAAGGGCAAAGAAGCTTTCATCTTCGATACCGACCGCTTTACCGTCTGAAATGATGTACGCCAGACCTCGAACAATCTGCATCGTTGCCAATGTCGTGATCAGAGCATTTATTTTCAGTTTGGCGATAACAAAGCCATTGACCAGGCCACAGAGAACGCCCAGCAACAACCCCGCTGCCACGCCAA  >  NZ_CP009273/1976450‑1976711

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