Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,304,953	A→G	S549P (TCG→CCG)	yjcS ←	alkyl sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,304,953	0	A	G	100.0%	178.0 / NA	51	S549P (TCG→CCG)	yjcS	alkyl sulfatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (25/26);  total (25/26)

AAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCA  >  NC_000913/4304822‑4305086                                                                                                                                    |                                                                                                                                      aaGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCcgc                                                                                                                               >  2:549867/1‑140 (MQ=255)        cTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGgt                                                                                                                        >  1:494041/1‑140 (MQ=255)         tATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCAt                                                                                                                                    <  1:580743/127‑1 (MQ=255)         tATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCAt                                                                                                                                    >  2:580743/1‑127 (MQ=255)                cAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCgggg                                                                                                                <  2:494041/140‑1 (MQ=255)                cAGCGTCAGGTTGAGGTTATCGCCGTTGCTCACATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCgggg                                                                                                                <  2:95192/140‑1 (MQ=255)                 aGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGa                                                                                                               >  2:243647/1‑140 (MQ=255)                   cGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAAc                                                                                                             >  1:368504/1‑140 (MQ=255)                     tCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGaa                                                                                                                            >  2:450363/1‑123 (MQ=255)                     tCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGaa                                                                                                                            <  1:450363/123‑1 (MQ=255)                     tCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAAc                                                                                                             <  1:519311/138‑1 (MQ=255)                     tCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAAc                                                                                                             >  2:519311/1‑138 (MQ=255)                     tCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCg                                                                                                           >  2:575951/1‑140 (MQ=255)                     tCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCg                                                                                                           >  1:307884/1‑140 (MQ=255)                     tCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCg                                                                                                           >  1:772286/1‑140 (MQ=255)                       aGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCggt                                                                                                         >  2:442094/1‑140 (MQ=255)                        ggTTGAGGTTATCGCCGTTGCTCATATT‑AAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCggtgg                                                                                                       >  2:829383/1‑140 (MQ=255)                             aGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCGTCCCGCGAATGGTGTCCGGGGAACCGGTGGTGcc                                                                                                   >  2:308444/1‑140 (MQ=255)                             aGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGcc                                                                                                   >  2:67110/1‑140 (MQ=255)                             aGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGcc                                                                                                   >  2:519094/1‑140 (MQ=255)                                    cGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTg                                                                                            <  1:549867/140‑1 (MQ=255)                                       cGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAAc                                                                                         <  1:342169/140‑1 (MQ=255)                                       cGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAAc                                                                                         <  1:778733/140‑1 (MQ=255)                                                 atTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCAccc                                                                               >  2:140245/1‑140 (MQ=255)                                                   tGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCt                                                                             <  2:665819/140‑1 (MQ=255)                                                   tGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCt                                                                             >  1:665819/1‑140 (MQ=255)                                                     aaGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTc                                                                           >  1:60531/1‑140 (MQ=255)                                                          cAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGAGGCTGAACTTATGCACCCCTTCGCGCa                                                                      <  1:829383/140‑1 (MQ=255)                                                           aGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAg                                                                     >  2:92884/1‑140 (MQ=255)                                                           aGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAg                                                                     <  1:519551/140‑1 (MQ=255)                                                                gATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCtt                                                                <  1:243647/140‑1 (MQ=255)                                                                          cccGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGt                                                      <  1:92884/140‑1 (MQ=255)                                                                          cccGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGt                                                      <  1:427679/140‑1 (MQ=255)                                                                          cccGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGt                                                      <  1:442094/140‑1 (MQ=255)                                                                          cccGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGt                                                      <  1:519094/140‑1 (MQ=255)                                                                          cccGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGt                                                      <  2:431187/140‑1 (MQ=255)                                                                              cAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAgg                                                  >  2:701896/1‑140 (MQ=255)                                                                                   ttCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGaa                                             <  2:307884/140‑1 (MQ=255)                                                                                      gcgcTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAAcc                                          <  2:772286/140‑1 (MQ=255)                                                                                        cttATCGAGCCAACCGGCCATAAGGTCGAACGGAATTCCGCCCGGAATCCCGCGAATGGTGCCCGGGGAACCGGTGGTGCCGTGCCTGACCTTATGCCCCCCTTCGCGCACCTCTTCCGCGCCGGTCAGGTAGAACCcgc                                        <  2:774382/138‑1 (MQ=255)                                                                                         cTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACcgcg                                       <  1:680002/140‑1 (MQ=255)                                                                                         cTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACcgcg                                       <  2:60531/140‑1 (MQ=255)                                                                                                     aCGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGa                           <  1:656675/140‑1 (MQ=255)                                                                                                      cGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGcc                                     >  2:462375/1‑129 (MQ=255)                                                                                                      cGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGcc                                     <  1:713086/129‑1 (MQ=255)                                                                                                      cGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGcc                                     <  1:462375/129‑1 (MQ=255)                                                                                                      cGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGcc                                     >  2:713086/1‑129 (MQ=255)                                                                                                      cGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGAc                          >  1:269845/1‑140 (MQ=255)                                                                                                      cGGCCATAAAGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGAc                          >  2:254450/1‑140 (MQ=255)                                                                                                             aaaGTCGAACAGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCt                   <  1:67110/140‑1 (MQ=255)                                                                                                                       aGCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAg         <  2:269845/140‑1 (MQ=255)                                                                                                                        gCATTTCGACCGGCATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGc        >  1:398789/1‑140 (MQ=255)                                                                                                                              cgaccggcATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCa  <  2:398789/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                      AAGCACGCTATCGTTCAGCGTCAGGTTGAGGTTATCGCCGTTGCTCATATTGAAGTTCAGGCTGATATTTTTACCCGCAGCTTTCGCGCTATCGAGGCGAACGGCCATAAAGTCGAACAGCATTTCGACCGACATCCCGCGAATGGTGTCCGGGGAACCGGTGGTGCCGTGGCTGAACTTATGCACCCCTTCGCGCAGCTCTTTCGCGCCGGTCAGGTAGAAACCGCGCCATGTGGCGGACTCGGCCTGATAGCCCAGCTGTTCA  >  NC_000913/4304822‑4305086

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: ‑