Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,564,565	A→G	L390P (CTG→CCG)	dosP ←	oxygen sensor, c‑di‑GMP phosphodiesterase, heme‑regulated; cold‑ and stationary phase‑induced bioflim regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,564,565	0	A	G	100.0%	170.5 / NA	48	L390P (CTG→CCG)	dosP	oxygen sensor, c‑di‑GMP phosphodiesterase, heme‑regulated; cold‑ and stationary phase‑induced bioflim regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (27/21);  total (27/21)

TCACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTT  >  NC_000913/1564418‑1564705                                                                                                                                                    |                                                                                                                                             tCACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGt                                                                                                                                           >  1:127660/1‑151 (MQ=255)  cACGACAGCGATTGACCACTTCCAGCAATGCCTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGCATATGGTCAACACCGATGAGATACACGACGGGAGAGACGGCTTTGTCGAGCAGGTCATCGAGGTAATTGTGCggggt                                                                                                                                          >  1:465920/1‑149 (MQ=255)  cACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCTCGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGtt                                                                                                                                          <  1:627725/151‑1 (MQ=255)                 ccACTTCCAGCAATGCCTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGc                                                                                                                             >  2:324880/1‑149 (MQ=255)                            aaTGCCTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCg                                                                                                                >  2:268488/1‑151 (MQ=255)                                ccTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTTTTGCGATTTGGCAGACCGGTCATCGGATc                                                                                                            <  2:244149/151‑1 (MQ=255)                                ccTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATc                                                                                                            >  2:354868/1‑151 (MQ=255)                                  tGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCa                                                                                                           >  2:50781/1‑150 (MQ=255)                                      cGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAAtt                                                                                                       <  2:693347/150‑1 (MQ=255)                                      cGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAAtt                                                                                                       <  2:651170/150‑1 (MQ=255)                                        gCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGa                                                                                                    <  1:621001/151‑1 (MQ=255)                                        gCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGa                                                                                                    <  1:487053/151‑1 (MQ=255)                                           cACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATg                                                                                                  >  1:283498/1‑150 (MQ=255)                                            aCGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGAAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGGTCAAATTGGAt                                                                                                   >  1:468967/1‑148 (MQ=255)                                                aTAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAgttgt                                                                                            >  2:203845/1‑151 (MQ=255)                                                 tAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAgttgtt                                                                                           >  2:388464/1‑151 (MQ=255)                                                  aGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTc                                                                                          <  1:555280/151‑1 (MQ=255)                                                          cTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCt                                                                                  >  2:501672/1‑151 (MQ=255)                                                               aatcaCATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACgg                                                                             <  1:718229/151‑1 (MQ=255)                                                                  cacaTCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCt                                                                           <  2:501421/150‑1 (MQ=255)                                                                     aTCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCttttt                                                                       <  1:108063/151‑1 (MQ=255)                                                                       ccTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTc                                                                     <  1:711958/151‑1 (MQ=255)                                                                         tGAATCTGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCATATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCt                                                                   >  1:456870/1‑151 (MQ=255)                                                                           aaTATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCCGt                                                                 <  2:461498/151‑1 (MQ=255)                                                                            atatGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGtt                                                                >  2:298272/1‑151 (MQ=255)                                                                             tatGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAGTTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTc                                                               >  1:748493/1‑151 (MQ=255)                                                                             tatGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTc                                                               >  2:726744/1‑151 (MQ=255)                                                                                   cAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcg                                                          <  2:351813/150‑1 (MQ=255)                                                                                    aaCACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgcc                                                        <  1:346434/151‑1 (MQ=255)                                                                                      cacCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgcc                                                        >  2:112993/1‑149 (MQ=255)                                                                                      cacCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgcc                                                        >  1:610625/1‑149 (MQ=255)                                                                                       acCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgccagc                                                     <  2:150913/151‑1 (MQ=255)                                                                                       acCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgccagc                                                     >  1:200980/1‑151 (MQ=255)                                                                                       acCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTccagcgccagc                                                     >  2:425087/1‑151 (MQ=255)                                                                                            tGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGc                                                >  2:743056/1‑151 (MQ=255)                                                                                              agaTACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGcc                                               >  1:317756/1‑150 (MQ=255)                                                                                                    accacGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATgct                                        <  1:374858/151‑1 (MQ=255)                                                                                                          ggAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATgctggctg                                   >  2:356636/1‑150 (MQ=255)                                                                                                            agagaCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGat                                 >  1:155308/1‑150 (MQ=255)                                                                                                              agaCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATc                              <  2:147246/151‑1 (MQ=255)                                                                                                                aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTg                            <  2:308873/151‑1 (MQ=255)                                                                                                                aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTg                            >  1:447542/1‑151 (MQ=255)                                                                                                                aCGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTg                            <  1:25916/151‑1 (MQ=255)                                                                                                                  ggCTTTGTCGACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGcc                          >  1:265503/1‑151 (MQ=255)                                                                                                                           gACCAGGTCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCg                 >  1:764771/1‑151 (MQ=255)                                                                                                                                  tCATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAgg          >  1:594568/1‑151 (MQ=255)                                                                                                                                   cATCGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGgc         <  2:238573/151‑1 (MQ=255)                                                                                                                                      cGAGGTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCt      <  2:742849/151‑1 (MQ=255)                                                                                                                                          gTAATTGTGCGGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGtt  >  2:337411/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                             TCACGAAAGCGATTGACCACTTCCAGCAATGCCTGATCGGCCCACGCATAGCCAAGGCTATCAATCACATCCTGAATATGGTCAACACCGATGAGATACACCACGGGAGAGACGGCTTTGTCGACCAGGTCATCGAGGTAATTGTGCAGGTTATTGCGATTTGGCAGACCGGTCATCGGATCAAATTGGATGAGTTGTTCAATATGCTGACGGCTTTTTTCCTGTTCCAGCGCCAGCGCGGCCATATGCTGGCTGATATCTGCCACGCGTTCGATAAAGGCGCTGGTT  >  NC_000913/1564418‑1564705

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