Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,495,050	C→T	*81* (TAG→TAA)	yfdY ←	DUF2545 family putative inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,495,050	0	C	T	100.0%	230.9 / NA	70	*81* (TAG→TAA)	yfdY	DUF2545 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/36);  total (34/36)

GAAATGCCGCCCCCTTTTCCGACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAA  >  NC_000913/2494917‑2495164                                                                                                                                      |                                                                                                                   gaaATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCg                                                                                                               <  2:1283092/139‑1 (MQ=255)     tGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTg                                                                                                           >  2:1189895/1‑139 (MQ=255)       ccgccCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGcc                                                                                                         <  1:219230/139‑1 (MQ=255)           ccccTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATAc                                                                                                     >  1:810521/1‑139 (MQ=255)                          aGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGa                                                                                      >  2:326134/1‑139 (MQ=255)                                      ttATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCca                                                                          >  2:1352422/1‑139 (MQ=255)                                              ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCgg                                                                  >  2:1041334/1‑139 (MQ=255)                                              ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCgg                                                                  >  2:579147/1‑139 (MQ=255)                                              ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCgg                                                                  >  1:287557/1‑139 (MQ=255)                                                  tgCGATAACGTCTTTCACATGGCAGGTGAGCAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCa                                                              >  1:1128882/1‑139 (MQ=255)                                                  tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCa                                                              >  1:33291/1‑139 (MQ=255)                                                  tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCa                                                              >  2:523385/1‑139 (MQ=255)                                                  tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCa                                                              <  1:1189895/139‑1 (MQ=255)                                                    cGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAg                                                            <  2:6507/139‑1 (MQ=255)                                                      aTAACGTCTTTCACATGGCAGGTGAGCAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGaa                                                          <  2:1128882/139‑1 (MQ=255)                                                       taacGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCa                                                                                                        >  2:391149/4‑92 (MQ=255)                                                       taacGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCa                                                                                                        <  1:391149/89‑1 (MQ=255)                                                       taacGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGaaa                                                         <  1:445374/136‑1 (MQ=255)                                                           gTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATa                                                     <  2:33291/139‑1 (MQ=255)                                                            tCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATaa                                                    >  2:232033/1‑139 (MQ=255)                                                               ttCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCca                                                                          >  1:1205824/1‑114 (MQ=255)                                                               ttCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCca                                                                          <  2:1205824/114‑1 (MQ=255)                                                                 cacaTGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAg                                                            >  2:777655/1‑126 (MQ=255)                                                                 cacaTGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAg                                                            <  1:777655/126‑1 (MQ=255)                                                                              aaCAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGtt                                  >  1:603175/1‑139 (MQ=255)                                                                                  aTCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAAtcatca                                                                               <  2:325896/90‑1 (MQ=255)                                                                                  aTCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAAtcatca                                                                               >  1:325896/1‑90 (MQ=255)                                                                                     aTGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTAcc                           >  1:1229715/1‑139 (MQ=255)                                                                                       gTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGc                         <  1:506464/139‑1 (MQ=255)                                                                                        tCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGcc                        >  2:1127411/1‑139 (MQ=255)                                                                                           gatgctgatgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGt                     <  1:881881/139‑1 (MQ=255)                                                                                           gatgctgatgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGt                     >  2:1014323/1‑139 (MQ=255)                                                                                             tgctgatgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATcccc                                            >  1:592176/1‑114 (MQ=255)                                                                                             tgctgatgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTa                   <  2:603175/139‑1 (MQ=255)                                                                                             tgctgatgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAAGCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATcccc                                            <  2:592176/114‑1 (MQ=255)                                                                                              gctgatgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTaa                  <  2:287557/139‑1 (MQ=255)                                                                                                 gatgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACca               >  1:1065964/1‑139 (MQ=255)                                                                                                  atgTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACcac              >  2:1042949/1‑139 (MQ=255)                                                                                                   tgtgAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACcaca             <  1:1042949/139‑1 (MQ=255)                                                                                                   tgtgAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACcaca             <  2:743049/139‑1 (MQ=255)                                                                                                       aaCACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCca                                                                          <  2:800698/74‑1 (MQ=255)                                                                                                       aaCACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCca                                                                          >  1:800698/1‑74 (MQ=255)                                                                                                       aaCACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCAc                                       <  2:760026/109‑1 (MQ=255)                                                                                                       aaCACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCAc                                       >  1:760026/1‑109 (MQ=255)                                                                                                       aaCACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGt         >  2:1130200/1‑139 (MQ=255)                                                                                                         cacCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGa       <  2:1065964/139‑1 (MQ=255)                                                                                                         cacCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGa       <  1:1352422/139‑1 (MQ=255)                                                                                                         cacCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGa       <  1:1130200/139‑1 (MQ=255)                                                                                                         cacCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGa       <  2:700116/139‑1 (MQ=255)                                                                                                         cacCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGa       <  1:523385/139‑1 (MQ=255)                                                                                                           ccTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAg     <  1:326134/139‑1 (MQ=255)                                                                                                            cTGCCCAAATGAGCCATCCTGACCCATTACCGCCTg                                                                                                           <  2:372013/36‑1 (MQ=38)                                                                                                            cTGCCCAAATGAGCCATCCTGACCCATTACCGCCTg                                                                                                           >  1:372013/1‑36 (MQ=38)                                                                                                            cTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAatcatc                                                                                >  1:271923/1‑63 (MQ=255)                                                                                                            cTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAatcatc                                                                                <  2:271923/63‑1 (MQ=255)                                                                                                              gCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGtaa  <  1:232033/139‑1 (MQ=255)                                                                                                                    aTGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAAcc                <  1:1352662/119‑1 (MQ=255)                                                                                                                    aTGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAAcc                >  2:1352662/1‑119 (MQ=255)                                                                                                                         ccATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAg                                                            >  1:383459/1‑70 (MQ=255)                                                                                                                         ccATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAg                                                            <  2:383459/70‑1 (MQ=255)                                                                                                                         ccATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGta   >  2:387098/1‑127 (MQ=255)                                                                                                                         ccATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGta   <  1:387098/127‑1 (MQ=255)                                                                                                                             ccTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCg                          <  1:1368178/100‑1 (MQ=255)                                                                                                                             ccTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCg                          >  2:1368178/1‑100 (MQ=255)                                                                                                                                 aCCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCca                                                                          >  1:825503/1‑48 (MQ=255)                                                                                                                                 aCCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCca                                                                          <  2:825503/48‑1 (MQ=255)                                                                                                                                 aCCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGa                                                           >  1:1250087/1‑63 (MQ=255)                                                                                                                                 aCCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGa                                                           <  2:1250087/63‑1 (MQ=255)                                                                                                                                 aCCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGt                     <  1:984702/101‑1 (MQ=255)                                                                                                                                 aCCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGt                     >  2:984702/1‑101 (MQ=255)                                                                                                                                      |                                                                                                                   GAAATGCCGCCCCCTTTTCCGACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAAGAAACATAAATATCCCCGTCACTAGTTCATTACCGCCAGTTAACCACATGGTGAAGTAA  >  NC_000913/2494917‑2495164

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