Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,062,912	G→A	W222* (TGG→TGA)	yihN →	putative transporter YihN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,062,912	0	G	A	100.0%	219.2 / NA	63	W222* (TGG→TGA)	yihN	putative transporter YihN
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/31);  total (32/31)

GTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACATCTATCTATTGCCCGTTGCGCTGGTGGGGGCTTACGGCATCATT  >  NC_000913/4062778‑4063032                                                                                                                                       |                                                                                                                         gtgATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTAt                                                                                                                     <  1:753889/140‑1 (MQ=255)       gCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGGCCTCGGTGCTGAAGGACAAAACTATCTGACTTATGGCttt                                                                                                               >  2:871220/1‑140 (MQ=255)         cGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTa                                                                                                             >  1:429545/1‑140 (MQ=255)         cGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTa                                                                                                             >  1:329509/1‑140 (MQ=255)           tAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAAc                                                                                                           >  1:531967/1‑140 (MQ=255)             ggCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCG‑AAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGtc                                                                                                        >  1:402703/1‑140 (MQ=255)                  tattattTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGtcttct                                                                                                    >  2:576884/1‑140 (MQ=255)                          tCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTAc                                                                                            <  1:576884/140‑1 (MQ=255)                          tCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTAc                                                                                            <  1:971520/140‑1 (MQ=255)                           cTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTAcg                                                                                           <  1:154526/140‑1 (MQ=255)                              tGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCgg                                                                                        >  1:643672/1‑140 (MQ=255)                              tGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCgg                                                                                        >  1:614272/1‑140 (MQ=255)                              tGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCG‑AAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTAcg                                                                                           >  1:118376/1‑136 (MQ=255)                              tGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCG‑AAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTAcg                                                                                           <  2:118376/136‑1 (MQ=255)                               gTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGt                                                                                       >  1:304011/1‑140 (MQ=255)                                          cAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTAc                                                                                            <  1:721823/124‑1 (MQ=255)                                          cAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTAc                                                                                            >  2:721823/1‑124 (MQ=255)                                          cAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGcc                                                                            <  2:329509/140‑1 (MQ=255)                                          cAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGcc                                                                            <  2:429545/140‑1 (MQ=255)                                              gaagaGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCG‑AAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACa                                                                       <  2:402703/140‑1 (MQ=255)                                               aagaGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACa                                                                       <  2:531967/140‑1 (MQ=255)                                                    gCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGTTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACAttctt                                                                  >  1:738253/1‑140 (MQ=255)                                                     cACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACt                                                                                  >  2:102194/1‑123 (MQ=255)                                                     cACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACt                                                                                  <  1:102194/123‑1 (MQ=255)                                                         gtccgtTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTcattc                                                             >  2:356147/1‑140 (MQ=255)                                                         gtccgtTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTcattc                                                             >  2:99450/1‑140 (MQ=255)                                                         gtccgtTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTcattc                                                             >  1:624996/1‑140 (MQ=255)                                                         gtccgtTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTcattc                                                             >  1:456009/1‑140 (MQ=255)                                                         gtccgtTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTcattc                                                             >  1:333421/1‑140 (MQ=255)                                                         gtccgtTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTcattc                                                             >  1:377868/1‑140 (MQ=255)                                                         gtccgtTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACCATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTcattc                                                             >  1:213256/1‑140 (MQ=255)                                                          tccgtTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACa                                                                       <  2:938829/129‑1 (MQ=255)                                                          tccgtTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACa                                                                       >  1:938829/1‑129 (MQ=255)                                                               tGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCattccattcc                                                       >  1:228376/1‑140 (MQ=255)                                                                   gTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGaa                                                   >  1:439858/1‑140 (MQ=255)                                                                     gAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTcattcca                                                           >  1:230907/1‑130 (MQ=255)                                                                     gAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTcattcca                                                           <  2:230907/130‑1 (MQ=255)                                                                      aaaaaaGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGaaaaa                                                <  2:72989/140‑1 (MQ=255)                                                                      aaaaaaGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGaaaaa                                                <  2:643672/140‑1 (MQ=255)                                                                       aaaaaGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAAc                                               >  1:210525/1‑140 (MQ=255)                                                                        aaaaGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGTTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACa                                              <  2:738253/140‑1 (MQ=255)                                                                          aagaagACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACatc                                            <  1:356147/140‑1 (MQ=255)                                                                          aagaagACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACatc                                            <  2:356296/140‑1 (MQ=255)                                                                          aagaagACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACatc                                            <  1:842105/140‑1 (MQ=255)                                                                             aagaCGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACatctat                                         >  2:353490/1‑140 (MQ=255)                                                                                aCGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACAtctatcta                                      <  2:304011/140‑1 (MQ=255)                                                                                aCGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACAtctatcta                                      <  1:353490/140‑1 (MQ=255)                                                                                                  ccAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACATCTATCTATTGCCCGTTGCGCtggtg                    <  2:228376/140‑1 (MQ=255)                                                                                                         gACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACATCTATCTATTGCCCg                               <  1:264502/122‑1 (MQ=255)                                                                                                         gACCTCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACATCTATCTATTGCCCg                               >  2:264502/1‑122 (MQ=255)                                                                                                             tCGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACATCTATCTATTGCCCGTTGCGCTGGTGGGGGCTTACgg         <  2:592525/140‑1 (MQ=255)                                                                                                              cGGTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACATCTATCTATTGCCCGTTGCGCTGGTGGGGGCTTACGGc        >  1:312583/1‑140 (MQ=255)                                                                                                                gTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACATCTATCTATTGCCCGTTGCGCTGGTGGGGGCTTACGGcat      >  1:264611/1‑140 (MQ=255)                                                                                                                gTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACATCTATCTATTGCCCGTTGCGCTGGTGGGGGCTTACGGc        >  1:766771/1‑138 (MQ=255)                                                                                                                gTGCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACATCTATCTATTGCCCGTTGCGCTGGTGGGGGCTTACGGc        <  2:766771/138‑1 (MQ=255)                                                                                                                  gCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACATCTATCTATTGCCCGTTGCGCTGGTGGGGGCTTACGGcatca    <  1:55727/140‑1 (MQ=255)                                                                                                                  gCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACATCTATCTATTGCCCGTTGCGCTGGTGGGGGCTTACGGcatca    <  2:952310/140‑1 (MQ=255)                                                                                                                  gCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACATCTATCTATTGCCCGTTGCGCTGGTGGGGGCTTACGGcatca    <  2:624996/140‑1 (MQ=255)                                                                                                                  gCTGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACATCTATCTATTGCCCGTTGCGCTGGTGGGGGCTTACGGcatca    <  2:439858/140‑1 (MQ=255)                                                                                                                  gCTGAAAGACAAAACCATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACATCTATCTATTGCCCGTTGCGCTGGTGGGGGCTTACGGcatca    <  2:213256/140‑1 (MQ=255)                                                                                                                    tGAAAGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACATCTATCTATTGCCCGTTGCGCTGGTGGGGGCTTACGGCATCAtt  >  2:214093/1‑140 (MQ=255)                                                                                                                      aaaGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACATCTATCTATTGCCCGTTGCGCTGGTGGGGGCTTACGGcatca    >  2:212463/1‑136 (MQ=255)                                                                                                                      aaaGACAAAACTATCTGACTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACATCTATCTATTGCCCGTTGCGCTGGTGGGGGCTTACGGcatca    <  1:212463/136‑1 (MQ=255)                                                                                                                                       |                                                                                                                         GTGATTGCCGTAGGCATTATTATTTTCTTTGTCCTGAATGACAAAGAAGAGGCACCGTCCGTTGAGGTGAAAAAAGAAGACGGAGCATCGAAAAACACCAGTATGACCTCGGTGCTGAAAGACAAAACTATCTGGCTTATCGCTTTTAACGTCTTCTTCGTTTACGCGGTTTACTGTGGCCTGACATTCTTCATTCCATTCCTGAAAAACATCTATCTATTGCCCGTTGCGCTGGTGGGGGCTTACGGCATCATT  >  NC_000913/4062778‑4063032

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