Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	1,558,439	A→G	Y518H (TAT→CAT)	dosP ←	oxygen‑sensing cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	1,558,439	0	A	G	100.0%	197.2 / NA	52	Y518H (TAT→CAT)	dosP	oxygen‑sensing cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (26/26);  total (26/26)

GCGCCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAA  >  NZ_CP009273/1558354‑1558526                                                                                      |                                                                                        gcgcCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATcc                                                                                     >  2:2895723/1‑90 (MQ=255)      tAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcattgc                                                                                >  1:1968446/1‑90 (MQ=255)          cccAACCGCTTTTTTCCCTTTTCGTTCATCGCGGGGTTGAAAACCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcattgcattg                                                                            <  1:1785185/89‑1 (MQ=255)          acaaaaCCCTTTTTTCCCTTTTCTTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcattgcattg                                                                            <  1:851950/87‑1 (MQ=255)             aaaCGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGa                                                                         >  2:1262331/1‑90 (MQ=255)             aaaCGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGa                                                                         >  1:173220/1‑90 (MQ=255)              aaCGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcattgcattgctgt                                                                        >  1:2831619/1‑86 (MQ=255)              aaCGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcattg                                                                                 <  2:1917009/81‑1 (MQ=255)              aaCGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcattg                                                                                 >  1:1917009/1‑81 (MQ=255)               ccGCTTTTTTCCCTTTTCGTTCACCGCCGGGTTGAAGAACCGCCAGCCTTTCCCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGc                                                                       <  1:1804086/89‑1 (MQ=255)               aCGCTTTTTTCCCTTTTCTTTCATCGCCGGGCTAAAAAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGc                                                                       <  1:640203/90‑1 (MQ=255)               aCGCTTTTTTACCTTTTCGTTCATCCCGGGGTTGAAAAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGc                                                                       <  2:2957015/90‑1 (MQ=255)               aCGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAAAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGc                                                                       <  2:2303458/90‑1 (MQ=255)                    ctTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTgg                                                                  >  1:763912/1‑90 (MQ=255)                            aTTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAg                                                          >  1:2001577/1‑90 (MQ=255)                            aTTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAg                                                          >  2:1451720/1‑90 (MQ=255)                              ttCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATtgtg                                                                          <  1:644567/72‑1 (MQ=255)                              ttCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATtgtg                                                                          >  2:644567/1‑72 (MQ=255)                              ttCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTa                                                        >  2:2213600/1‑90 (MQ=255)                              ttCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTa                                                        >  2:988327/1‑90 (MQ=255)                              ttCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTa                                                        >  1:427443/1‑90 (MQ=255)                              ttCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTa                                                        >  2:2786014/1‑90 (MQ=255)                                 gTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATc                                                     >  2:2187981/1‑90 (MQ=255)                                    cATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcgcg                                                  >  2:550366/1‑90 (MQ=255)                                      tCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGt                                                >  1:1500509/1‑90 (MQ=255)                                              cTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTAccc                                        <  2:772109/90‑1 (MQ=255)                                                gaagaaCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAg                                      >  1:2380177/1‑90 (MQ=255)                                                gaagaaCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAg                                      <  1:2187981/90‑1 (MQ=255)                                                   gaaCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTc                                   <  2:1500509/90‑1 (MQ=255)                                                     aCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGt                                 >  2:2607235/1‑90 (MQ=255)                                                     aCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGt                                 >  2:515395/1‑90 (MQ=255)                                                     aCTCCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGt                                 <  2:2900776/90‑1 (MQ=255)                                                      cTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTa                                >  2:1714575/1‑90 (MQ=255)                                                        gCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTAccc                                        >  1:1102103/1‑80 (MQ=255)                                                        gCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTAccc                                        <  2:1102103/80‑1 (MQ=255)                                                        gCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGc                              >  1:2866756/1‑90 (MQ=255)                                                            gCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGAt                          <  2:1047326/90‑1 (MQ=255)                                                                ttACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCa                      <  2:2652976/90‑1 (MQ=255)                                                                ttACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCa                      <  2:2812875/90‑1 (MQ=255)                                                                 tACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCaa                     <  2:2164408/90‑1 (MQ=255)                                                                 tACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCaa                     <  1:550366/90‑1 (MQ=255)                                                                  aCCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAAt                    >  2:872376/1‑90 (MQ=255)                                                                    cGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATAc                  <  1:2042358/90‑1 (MQ=255)                                                                      ccATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTAccc                                        >  1:102118/1‑66 (MQ=255)                                                                      ccATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTAccc                                        <  2:102118/66‑1 (MQ=255)                                                                      ccATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTc                <  1:2436868/90‑1 (MQ=255)                                                                      ccATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCCCGGTTTTTACCCAGGTCGTACCTGATGCCAATACTc                <  1:2452789/90‑1 (MQ=255)                                                                       cATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCa               <  2:427443/90‑1 (MQ=255)                                                                        aTTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCaa              <  2:966362/90‑1 (MQ=255)                                                                             tGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTaa         <  1:1356158/90‑1 (MQ=255)                                                                               cGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGc       <  1:872376/90‑1 (MQ=255)                                                                                  atatGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCggg    >  1:1234859/1‑90 (MQ=255)                                                                                    ataATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGaa  <  1:2150897/90‑1 (MQ=255)                                                                                    ataATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGaa  <  1:1146447/90‑1 (MQ=255)                                                                                      |                                                                                        GCGCCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGGAA  >  NZ_CP009273/1558354‑1558526

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