Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,637,373	T→C	V226V (GTT→GTC)	creC →	sensory histidine kinase in two‑component regulatory system with CreB or PhoB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,637,373	0	T	C	100.0%	172.5 / NA	46	V226V (GTT→GTC)	creC	sensory histidine kinase in two‑component regulatory system with CreB or PhoB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/29);  total (17/29)

AGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAA  >  NC_000913/4637237‑4637498                                                                                                                                         |                                                                                                                              agcgGCGAATTTTATGGGCCAGCGCTATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTccc                                                                                                                             <  1:492052/139‑1 (MQ=255)       gAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTccc                                                                                                                       <  2:360254/139‑1 (MQ=255)        aaTTTTATGGGCCAGCGCCATTTTGTTGCGGCTTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCCACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTcccc                                                                                                                      <  1:547873/139‑1 (MQ=255)         aTTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCg                                                                                                                     <  2:1536/139‑1 (MQ=255)                    cAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGgtagt                                                                                                          >  2:668824/1‑139 (MQ=255)                          taTTTTGTAGGGGATTGCGCTGCTGCTTGGCGCGGGGATGGTCTGGTGCCTCAACTGCCCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTGGGTAGTAGCGAg                                                                                                    <  1:1108454/138‑1 (MQ=255)                          cATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAg                                                                                                    <  2:83840/139‑1 (MQ=255)                          cATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAg                                                                                                    <  1:610328/139‑1 (MQ=255)                                         tGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTcgcg                                                                                     >  1:290388/1‑139 (MQ=255)                                            aCTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAg                                                                                  <  1:668824/139‑1 (MQ=255)                                              tGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGgc                                                                                <  2:401587/139‑1 (MQ=255)                                              tGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGgc                                                                                <  2:190499/139‑1 (MQ=255)                                               ggTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGgcg                                                                               <  2:805825/139‑1 (MQ=255)                                                  gATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTg                                                                            <  1:1011202/139‑1 (MQ=255)                                                   aTTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTgg                                                                           >  2:97429/1‑139 (MQ=255)                                                     tGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTcccc                                                                                                                      >  2:590186/1‑94 (MQ=255)                                                     tGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTcccc                                                                                                                      <  1:590186/94‑1 (MQ=255)                                                                       gtggATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTcgctactaccg                                                                              >  2:887041/1‑108 (MQ=255)                                                                       gtggATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTcgctactaccg                                                                              <  1:887041/116‑9 (MQ=255)                                                                       gtggATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTg                                                       >  1:962475/1‑139 (MQ=255)                                                                            tCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAgg                                                  >  1:676514/1‑139 (MQ=255)                                                                             cAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAggg                                                 <  2:499065/139‑1 (MQ=255)                                                                                        ttGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAggg                                                 <  2:95029/128‑1 (MQ=255)                                                                                        ttGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAggg                                                 >  1:95029/1‑128 (MQ=255)                                                                                        ttGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTatat                                      <  2:652961/139‑1 (MQ=255)                                                                                        ttGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTatat                                      <  1:466261/139‑1 (MQ=255)                                                                                         tGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATAtt                                     <  1:45048/139‑1 (MQ=255)                                                                                            ccGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAg                                  <  2:1043826/139‑1 (MQ=255)                                                                                            ccGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAg                                  <  2:64078/139‑1 (MQ=255)                                                                                              ggCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCa                                <  2:641765/139‑1 (MQ=255)                                                                                               gCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAg                               <  1:97429/139‑1 (MQ=255)                                                                                                    cTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGt                          <  2:924309/139‑1 (MQ=255)                                                                                                           tGCTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCa                   >  2:814250/1‑139 (MQ=255)                                                                                                             cTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACt                                          <  1:167580/114‑1 (MQ=255)                                                                                                             cTGATTCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACt                                          >  2:167580/1‑114 (MQ=255)                                                                                                                  tCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTc            <  1:1008411/139‑1 (MQ=255)                                                                                                                  tCCGTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTc            <  1:814250/139‑1 (MQ=255)                                                                                                                     gTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGaaaaa                                            <  2:612190/104‑1 (MQ=255)                                                                                                                     gTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGaaaaa                                            >  1:612190/1‑104 (MQ=255)                                                                                                                     gTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAg                               <  2:723399/117‑1 (MQ=255)                                                                                                                     gTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAg                               >  1:723399/1‑117 (MQ=255)                                                                                                                     gTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATg         >  1:464157/1‑139 (MQ=255)                                                                                                                     gTCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATg         >  2:965747/1‑139 (MQ=255)                                                                                                                      tCACTGACAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGa        >  2:533452/1‑139 (MQ=255)                                                                                                                            aCAATAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTaaa  >  1:834276/1‑139 (MQ=255)                                                                                                                               aTAAGCCCGTCCCTCTCCCCGGTCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAggg                                                 <  2:816848/89‑1 (MQ=255)                                                                                                                               aTAAGCCCGTCCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAggg                                                 >  1:816848/1‑89 (MQ=255)                                                                                                                                         |                                                                                                                              AGCGGCGAATTTTATGGGCCAGCGCCATTTTGTTGGGGATTGCACTGGTGATTGGCGCAGGCATGGTTTGGTGGATCAACCGCTCTATTGCCCGGCTCACTCGCTATGCTGATTCCGTCACTGACAATAAGCCCGTTCCTCTCCCCGATCTCGGTAGTAGCGAGTTGCGTAAACTCGCGCAGGCGCTGGAAAGTATGCGCGTGAAGCTGGAAGGGAAAAACTATATTGAGCAGTATGTTTACGCATTAACTCATGAGCTAAA  >  NC_000913/4637237‑4637498

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