Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,047,787	C→T	T950T (ACC→ACT)	yeeJ →	inverse autotransporter adhesin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,047,787	0	C	T	100.0%	172.6 / NA	48	T950T (ACC→ACT)	yeeJ	inverse autotransporter adhesin
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (28/20);  total (28/20)

CTGACCAGTTTGAAAAATGGTGATTATAGGGTTACGGCCTCTGTGAGCTCTGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACCGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATAAAGAAATCACCTTCTCTGTGCCAAACGACGTCGCAAGTAAGTTCTCGATTAGCAACGGAGGAAA  >  NC_000913/2047650‑2047921                                                                                                                                          |                                                                                                                                       cTGACCAGTTTGAAAAATGGTGATTATAGGGTTACGGCCTCTGTGAGCTCTGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGt                                                                                                                                      >  1:299230/1‑140 (MQ=255)      cAGTTTGAAAAATGGTGATTATAGGGTTACGGCCTCTGTGAGCTCTGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCa                                                                                                                                 >  1:616464/1‑140 (MQ=255)            gAAAAATGGTGATTATAGGGTTACGGCCTCTGTGAGCTCTGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAg                                                                                                                           <  1:274391/140‑1 (MQ=255)            gAAAAATGGTGATTATAGGGTTACGGCCTCTGTGAGCTCTGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAg                                                                                                                           <  2:724436/140‑1 (MQ=255)               aaaTGGTGATTATAGGGTTACGGCCTCTGTGAGCTCTGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTc                                                                                                                        <  2:303724/140‑1 (MQ=255)                 aTGGTGATTATAGGGTTACGGCCTCTGTGAGCTCTGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCg                                                                                                                      <  2:253387/140‑1 (MQ=255)                                 tACGGCCTCTGTGAGCTCTGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAatatat                                                                                                              >  1:321462/1‑132 (MQ=255)                                 tACGGCCTCTGTGAGCTCTGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAatatat                                                                                                              <  2:321462/132‑1 (MQ=255)                                 tACGGCCTCTGTGAGCTCTGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCaa                                                                                                      >  1:411401/1‑140 (MQ=255)                                  aCGGCCTCTGTGAGCTCTGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAAc                                                                                                     >  2:858805/1‑140 (MQ=255)                                  aCGGCCTCTGTGAGCTCTGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAAc                                                                                                     >  2:404165/1‑140 (MQ=255)                                   cGGCCTCTGTGAGCTCTGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCg                                                                                                                      <  1:383245/122‑1 (MQ=255)                                   cGGCCTCTGTGAGCTCTGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCg                                                                                                                      >  2:383245/1‑122 (MQ=255)                                   cGGCCTCTGTGAGCTCTGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAAcc                                                                                                    >  2:230152/1‑140 (MQ=255)                                   cGGCCTCTGTGAGCTCTGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAAcc                                                                                                    >  2:584320/1‑140 (MQ=255)                                   cGGCCTCTGTGAGCTCTGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAAcc                                                                                                    >  1:763687/1‑140 (MQ=255)                                        tctGTGAGCTCTGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCa                                                                                               >  2:297648/1‑140 (MQ=255)                                               gCTCTGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATaaa                                                                                        <  2:411401/140‑1 (MQ=255)                                                  ctGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTg                                                                                                 >  1:71243/1‑128 (MQ=255)                                                  ctGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTg                                                                                                 <  2:71243/128‑1 (MQ=255)                                                     gTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGc                                                                                  <  1:584320/140‑1 (MQ=255)                                                     gTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGc                                                                                  <  1:858805/140‑1 (MQ=255)                                                     gTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGc                                                                                  <  1:230152/140‑1 (MQ=255)                                                         ccAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAAcc                                                                              <  1:297648/140‑1 (MQ=255)                                                             gCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACt                                                                          <  2:616464/140‑1 (MQ=255)                                                               tAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTaa                                                                        >  2:486287/1‑140 (MQ=255)                                                                     aCAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATa                                                                  <  2:299230/140‑1 (MQ=255)                                                                       aGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATaaa                                                                >  1:768272/1‑140 (MQ=255)                                                                          tGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATAAAGaa                                                             <  2:768272/140‑1 (MQ=255)                                                                          tGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATAAAGaa                                                             <  1:404165/140‑1 (MQ=255)                                                                                   tCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATAAAGAAATCACCTtc                                                    >  1:77823/1‑140 (MQ=255)                                                                                         aTCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATAAAGAAATCACCTTCTCtgtg                                              >  2:762906/1‑140 (MQ=255)                                                                                                   tgctgcCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATAAAGAAATCACCTTCTCTGTGCCAAacgacg                                    >  2:229358/1‑140 (MQ=255)                                                                                                   tgctgcCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATAAAGAAATCACCTTCTCTGTGCCAAacgacg                                    >  2:853742/1‑140 (MQ=255)                                                                                                     ctgcCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATAAAGAAATCACCTTCTCTGTGCCAAACTACGTc                                  >  1:249611/1‑140 (MQ=255)                                                                                                     ctgcCCTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATAAAGAAATCACCTTCTCTGTGCCAAACGACGTc                                  <  1:761464/140‑1 (MQ=255)                                                                                                          cTGACCCTCAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATAAAGAAATCACCTTCTCTGTGCCAAACGACGTCGCaag                             >  1:459802/1‑140 (MQ=255)                                                                                                                  cAGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATAAAGAAATCACCTTCTCTGTGCCAAACGACGTCGCAAGTAAGTtct                     <  1:853742/140‑1 (MQ=255)                                                                                                                   aGTGTGCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTGAAAGATAAAGAAATCACCTTCTCTGTGCCAAACGACGTCGCAAGTAAGTtctc                    >  2:753779/1‑140 (MQ=255)                                                                                                                       tgCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATAAAGAAATCACCTTCTCTGTGCCAAACGACGTCGCAAGTAAGTTCTCGAtt                >  1:38083/1‑140 (MQ=255)                                                                                                                        gCCTTCAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATAAAGAAATCACCTTCTCTGTGCCAAACGACGTCGCAAGTAAGTTCTCGATTa               >  2:482680/1‑140 (MQ=255)                                                                                                                             cAGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATAAAGAAATCACCTTCTCTGTGCCAAACGACGTCGCAAGTAAGTTCTCGATTAGCAAc          <  1:482680/140‑1 (MQ=255)                                                                                                                              aGGTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATAAAGAAATCACCTTCTCTGTGCCAAACGACGTCGCAAGTAAGTTCTCGATTAGCAACg         <  2:77823/140‑1 (MQ=255)                                                                                                                               ggTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATAAAGAAATCACCTTCTCTGTGCCAAACGACGTCGCAAGTAAGTTCTCGATTAGCAACgg        >  1:372484/1‑140 (MQ=255)                                                                                                                               ggTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATAAAGAAATCACCTTCTCTGTGCCAAACGACGTCGCAAGTAAGTTCTCGATTAGCAACgg        >  2:733266/1‑140 (MQ=255)                                                                                                                               ggTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATAAAGAAATCACCTTCTCTGTGCCAAACGACGTCGCAAGTAAGTTCTCGATTAGCAACgg        >  1:518664/1‑140 (MQ=255)                                                                                                                               ggTGATATCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATAAAGAAATCACCTTCTCTGTGCCAAACGACGTCGCAAGTAAGTTCTCGATTAGCAACgg        >  2:811362/1‑140 (MQ=255)                                                                                                                                     atCACTGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATAAAGAAATCACCTTCTCTGTGCCAAACGACGTCGCAAGTAAGTTCTCGATTAGCAACGGAGGaaa  >  2:719503/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                       CTGACCAGTTTGAAAAATGGTGATTATAGGGTTACGGCCTCTGTGAGCTCTGGTTCCCAGGCTAATCAACAGGTGAATTTTATCGGTGATCAAAGTACTGCTGCCCTGACCCTCAGTGTGCCTTCAGGTGATATCACCGTCACCAACACAGCTCCGCAATATATGACTGCAACCTTGCAGGATAAAAATGGCAACCCACTAAAAGATAAAGAAATCACCTTCTCTGTGCCAAACGACGTCGCAAGTAAGTTCTCGATTAGCAACGGAGGAAA  >  NC_000913/2047650‑2047921

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