Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,576,634	A→G	R79R (CGA→CGG)	BW25113_RS17825 →	YrhA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,576,634	0	A	G	100.0%	180.4 / NA	48	R79R (CGA→CGG)	BW25113_RS17825	YrhA family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/18);  total (30/18)

ATGAAGTTAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAA  >  NZ_CP009273/3576549‑3576719                                                                                      |                                                                                      aTGAAGTTAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGgttat                                                                                   >  1:111580/1‑90 (MQ=255)   gAAGTTAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTAttt                                                                                 >  1:1937535/1‑90 (MQ=255)   gAAGTTAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTAttt                                                                                 >  1:2877327/1‑90 (MQ=255)       ttAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGcc                                                                             <  2:547096/90‑1 (MQ=255)        tAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCct                                                                            >  1:5510/1‑90 (MQ=255)        tAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCct                                                                            >  1:1946559/1‑90 (MQ=255)            aTTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCg                                                                        <  1:318609/90‑1 (MQ=255)                tttttAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTc                                                                    >  2:2137084/1‑90 (MQ=255)                  tttAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCa                                                                  >  1:1902451/1‑90 (MQ=255)                   ttAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAg                                                                 <  1:2764848/90‑1 (MQ=255)                   ttAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAg                                                                 <  2:934625/90‑1 (MQ=255)                         aTCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCtt                                                           >  2:701725/1‑90 (MQ=255)                           cAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCa                                                         <  1:1850175/90‑1 (MQ=255)                           cAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCa                                                         <  2:1684862/90‑1 (MQ=255)                           cAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCa                                                         <  2:5510/90‑1 (MQ=255)                           cAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCa                                                         <  2:781727/90‑1 (MQ=255)                                 gATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTaaa                                                   >  1:1854758/1‑90 (MQ=255)                                      tCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCt                                              >  2:1347182/1‑90 (MQ=255)                                       cTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCtt                                             >  2:984779/1‑90 (MQ=255)                                        ttACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCttt                                            >  1:2782556/1‑90 (MQ=255)                                        ttACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCttt                                            >  2:1022793/1‑90 (MQ=255)                                         tACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCtttt                                           >  1:93860/1‑90 (MQ=255)                                          aCTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCttttt                                          >  1:1715965/1‑90 (MQ=255)                                           cTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTg                                         >  2:65312/1‑90 (MQ=255)                                           cTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTg                                         >  2:1312601/1‑90 (MQ=255)                                                   agagCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAAt                                 >  1:2443857/1‑90 (MQ=255)                                                    gagCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATg                                <  1:272667/90‑1 (MQ=255)                                                        aaTGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAAtt                            >  2:1223320/1‑90 (MQ=255)                                                        aaTGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAAtt                            <  2:930052/90‑1 (MQ=255)                                                             atggatTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTtata                       <  1:412938/90‑1 (MQ=255)                                                              tggatTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAg                                                                 <  1:2300806/47‑1 (MQ=255)                                                              tggatTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAg                                                                 >  2:57801/1‑47 (MQ=255)                                                              tggatTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAg                                                                 >  2:2300806/1‑47 (MQ=255)                                                              tggatTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAg                                                                 <  1:57801/47‑1 (MQ=255)                                                                gatTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTtata                       >  1:1392407/1‑87 (MQ=255)                                                                gatTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTtata                       <  2:1392407/87‑1 (MQ=255)                                                                 atTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGaaa                   >  1:2188912/1‑90 (MQ=255)                                                                   tCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAaata                 >  1:605098/1‑90 (MQ=255)                                                                     gAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAaataat               >  1:635817/1‑90 (MQ=255)                                                                          tAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTGAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAAtgatga          >  2:1182787/1‑90 (MQ=255)                                                                          tAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAAtgatga          >  2:90463/1‑90 (MQ=255)                                                                           aaTGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATg                                >  2:2432013/1‑67 (MQ=255)                                                                           aaTGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATg                                <  1:2432013/67‑1 (MQ=255)                                                                           aaTGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATgatgat         >  1:2158817/1‑90 (MQ=255)                                                                                 ttACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATa   <  2:2653444/90‑1 (MQ=255)                                                                                 ttACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATa   <  1:2276564/90‑1 (MQ=255)                                                                                  tACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATaa  <  2:1240217/90‑1 (MQ=255)                                                                                  tACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATaa  >  1:237991/1‑90 (MQ=255)                                                                                  tACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATaa  >  1:1531973/1‑90 (MQ=255)                                                                                      |                                                                                      ATGAAGTTAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAA  >  NZ_CP009273/3576549‑3576719

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