Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,428,964	T→C	E39G (GAG→GGG)	ytfJ ←	YtfJ family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,428,964	0	T	C	100.0%	126.8 / NA	34	E39G (GAG→GGG)	ytfJ	YtfJ family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/10);  total (24/10)

GCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGCAGGCA  >  NZ_CP009273/4428877‑4429049                                                                                        |                                                                                      gCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTCGCTAAACTGATCTTTATCAAGCACCAACCCGc                                                                                    >  1:772773/1‑91 (MQ=255) gCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACCCGc                                                                                    <  2:861677/91‑1 (MQ=255) gCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACCCGc                                                                                    >  1:461529/1‑91 (MQ=255)   agtatGCTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTACCAATCACCAACCCGccc                                                                                  <  2:162018/89‑1 (MQ=255)      atGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACCCGCCCCGa                                                                               <  1:156164/91‑1 (MQ=255)         ttGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACCCGCCCCGATCg                                                                            >  1:420202/1‑91 (MQ=255)           gCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACCCGCCCCGATCGGt                                                                          >  1:723273/1‑91 (MQ=255)           gCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACCCGCCCCGATCGGt                                                                          <  1:861508/91‑1 (MQ=255)           gCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACCCGCCCCGATCGGt                                                                          >  1:686862/1‑91 (MQ=255)                  tCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACCCGCCCCGATCGGTAATGCCa                                                                   >  1:719616/1‑91 (MQ=255)                   cGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACCCGCCCCGATCGGTAATGCCaa                                                                  >  2:264384/1‑91 (MQ=255)                     cACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACCCGCCCCGATCGGTAATGCCAATc                                                                >  1:745792/1‑91 (MQ=255)                     cACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACCCGCCCCGATCGGTAATGCCAATc                                                                >  1:263582/1‑91 (MQ=255)                          ttCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACCCGCCCCGATCGGTAATGCCAATcggcgg                                                           >  1:205135/1‑91 (MQ=255)                          ttCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACCCGCCCCGATCGGTAATGCCAATcggcgg                                                           >  2:802075/1‑91 (MQ=255)                              cACTAACTGCGCGCTGTGCCGGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACCCGCCCCGATCGGTAATGCCAATCGGCGGCACt                                                       <  2:163359/91‑1 (MQ=255)                               aCTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACCCGCCCCGATCGGTAATGCCAATCGGCGGCACTc                                                      >  1:387792/1‑91 (MQ=255)                                      gcgcgcTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACCCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGAcc                                               >  1:273956/1‑91 (MQ=255)                                           cTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACCCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGttt                                          <  2:273956/91‑1 (MQ=255)                                                  aGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACCCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTg                                   >  2:104824/1‑91 (MQ=255)                                                     ttttGTAGCTAAACTGATCTTTATCAAGCACCAACCCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGAcc                                               >  1:708/1‑76 (MQ=255)                                                     ttttGTAGCTAAACTGATCTTTATCAAGCACCAACCCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGAcc                                               <  2:708/76‑1 (MQ=255)                                                       ttGTAGCTAAACTGATCTTTATCAAGCACCAACCCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGAtgtg                              >  1:859833/1‑91 (MQ=255)                                                       ttGTAGCTAAACTGATCTTTATCAAGCACCAACCCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGAtgtg                              >  1:21605/1‑91 (MQ=255)                                                       ttGTAGCTAAACTGATCTTTATCAAGCACCAACCCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGAtgtg                              >  1:288702/1‑91 (MQ=255)                                                              tAAACTGATCTTTATCAAGCACCAACCCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGc                       <  1:694859/91‑1 (MQ=255)                                                                   tGATCTTTATCAAGCACCAACCCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCcatca                  <  2:686862/91‑1 (MQ=255)                                                                      tCTTTATCAAGCACCAACCCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGAcc                                               <  1:663606/59‑1 (MQ=255)                                                                      tCTTTATCAAGCACCAACCCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGAcc                                               >  2:663606/1‑59 (MQ=255)                                                                      tCTTTATCAAGCACCAACCCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCg               >  1:521709/1‑91 (MQ=255)                                                                              aaGCACCAACCCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAAcagc       >  1:206613/1‑91 (MQ=255)                                                                                gCACCAACCCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCcatcat                 <  2:366223/79‑1 (MQ=255)                                                                                gCACCAACCCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCcatcat                 >  1:366223/1‑79 (MQ=255)                                                                                gCACCAACCCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAAcagcag     >  2:401259/1‑91 (MQ=255)                                                                                   ccAACCCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAgcaggca  >  2:202244/1‑91 (MQ=255)                                                                                        |                                                                                      GCAATATGTTGCAGTACTCGCACTTTTCCCACTAACTGCGCGCTGTTCCAGGTTTTGTAGCTAAACTGATCTTTATCAAGCACCAACTCGCCCCGATCGGTAATGCCAATCGGCGGCACTCGCTGACCGGTTTCGAACTGATGTGCGGAAGCCATCATCGGCAACAGCAGGCA  >  NZ_CP009273/4428877‑4429049

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