Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,727,102	A→G	V87A (GTT→GCT)	yiaA ←	conserved inner membrane protein YiaA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,727,102	0	A	G	100.0%	209.2 / NA	59	V87A (GTT→GCT)	yiaA	conserved inner membrane protein YiaA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (31/28);  total (31/28)

CCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAACAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCTTTTGATAAGACGCCGCAGAAAACAGTCCTAATACCAGTACGGCAAAATAATATC  >  NC_000913/3726976‑3727237                                                                                                                               |                                                                                                                                        ccGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATaaa                                                                                                                            <  1:260843/140‑1 (MQ=255)   ggCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGa                                                                                                                          >  1:75054/1‑140 (MQ=255)    gCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGAc                                                                                                                         >  1:616990/1‑140 (MQ=255)    gCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGAc                                                                                                                         >  1:312470/1‑140 (MQ=255)     cATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGcagcag                                                                                                                                      >  1:590757/1‑126 (MQ=255)     cATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGcagcag                                                                                                                                      <  2:590757/126‑1 (MQ=255)     cATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACa                                                                                                                        >  1:213184/1‑140 (MQ=255)          cGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTcag                                                                                                                   <  1:456631/140‑1 (MQ=255)           gAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTcagg                                                                                                                  <  2:206996/140‑1 (MQ=255)            aaTATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTcaggc                                                                                                                 <  2:632678/140‑1 (MQ=255)             atatTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTcaggca                                                                                                                >  2:651174/1‑140 (MQ=255)               atTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTcaggcagg                                                                                                              <  2:60172/140‑1 (MQ=255)               atTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTcaggcagg                                                                                                              <  2:506709/140‑1 (MQ=255)               atTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTcaggcagg                                                                                                              <  2:616990/140‑1 (MQ=255)                      tGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCAtata                                                                                                       <  1:377848/140‑1 (MQ=255)                         aCCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATataata                                                                                                    >  2:60939/1‑140 (MQ=255)                           cGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGCAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATaaa                                                                                                  <  1:242886/140‑1 (MQ=255)                                   cACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTg                                                                                          >  1:406677/1‑140 (MQ=255)                                     ccAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTg                                                                                          <  2:484008/138‑1 (MQ=255)                                     ccAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTg                                                                                          >  1:484008/1‑138 (MQ=255)                                       aaaCAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCg                                                                                      >  2:485281/1‑140 (MQ=255)                                        aaCAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGt                                                                                        <  1:254068/137‑1 (MQ=255)                                        aaCAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGt                                                                                        >  2:254068/1‑137 (MQ=255)                                         aCAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGACAATAAAGACAGTCAGGCAGGTCATATAATAAACGGAAGTGGTCGGt                                                                                    <  1:317737/140‑1 (MQ=255)                                          cAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTa                                                                                   <  2:406677/140‑1 (MQ=255)                                          cAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTa                                                                                   <  2:213184/140‑1 (MQ=255)                                          cAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTa                                                                                   <  1:138624/140‑1 (MQ=255)                                                tAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGt                                                                                           <  1:195563/126‑1 (MQ=255)                                                tAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGt                                                                                           >  2:195563/1‑126 (MQ=255)                                                tAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCtt                                                                             >  2:622052/1‑140 (MQ=255)                                                tAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCtt                                                                             >  2:287011/1‑140 (MQ=255)                                                tAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCtt                                                                             <  1:1756/140‑1 (MQ=255)                                                       aaaGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACtt                                                                      >  2:610917/1‑140 (MQ=255)                                                        aaGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTc                                                                                       <  1:699295/122‑1 (MQ=255)                                                        aaGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTc                                                                                       >  2:699295/1‑122 (MQ=255)                                                        aaGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTg                                                                     >  2:632179/1‑140 (MQ=255)                                                              gTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGa                                                               >  1:630871/1‑140 (MQ=255)                                                               tCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATAc                                                                        >  1:35560/1‑130 (MQ=255)                                                               tCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATAc                                                                        <  2:35560/130‑1 (MQ=255)                                                               tCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGaa                                                              >  2:126979/1‑140 (MQ=255)                                                                ccATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAAc                                                             <  1:656946/140‑1 (MQ=255)                                                                        aCCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCtttt                                                     >  1:550062/1‑140 (MQ=255)                                                                                      gTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCTTTTGATAAGACGCCGCa                                       <  2:630871/140‑1 (MQ=255)                                                                                        aataatGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCTTTTGATAAGACGCCGCAGa                                     <  1:112478/140‑1 (MQ=255)                                                                                        aataatGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCTTTTGATAAGACGCCGCAGa                                     <  1:66643/140‑1 (MQ=255)                                                                                           aatGTCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCTTTTGATAAGACGCCGCAGaaaa                                  <  2:75054/140‑1 (MQ=255)                                                                                               tCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCTTTTGATAAGAcg                                            <  2:146843/126‑1 (MQ=255)                                                                                               tCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCTTTTGATAAGAcg                                            >  1:146843/1‑126 (MQ=255)                                                                                               tCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCTTTTGATAAGACGCCGCAGAAAACAGt                              >  1:644924/1‑140 (MQ=255)                                                                                               tCGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCTTTTGATAAGACGCCGCAGAAAACAGt                              >  1:713591/1‑140 (MQ=255)                                                                                                cGCATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCTTTTGATAAGACGCCGCAGAAAACAGTc                             <  1:651174/140‑1 (MQ=255)                                                                                                  cATTCCACAGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCTTTTGATAAGACGCCGCAGAAAACAGTCCt                           >  1:93062/1‑140 (MQ=255)                                                                                                       cacaGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCTTTTGATAAGACGCCGCAGAAAACAGTCCTAATAc                      <  1:746501/140‑1 (MQ=255)                                                                                                          aGACCTACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCTTTTGATAAGACGCCGCAGAAAACAGTCCTAATACCAg                   <  1:485281/140‑1 (MQ=255)                                                                                                               tACCATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCTTTTGATAAGACGCCGCAGAAAACAGTCCTAATACCAGTACgg              >  2:328195/1‑140 (MQ=255)                                                                                                                 ccATCAGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCTTTTGATAAGACGCCGCAGAAAACAGTCCTAATACCAGTACGGCa            >  1:286074/1‑140 (MQ=255)                                                                                                                      aGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCTTTTGATAAGACGCCGCAGAAAACAGTCCTAATACCAGTACGGCAAaata       >  1:527177/1‑140 (MQ=255)                                                                                                                      aGTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCTTTTGATAAGACGCCGCAGAAAACAGTCCTAATACCAGCACGGCAAaata       >  2:763423/1‑140 (MQ=255)                                                                                                                       gTAATGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCTTTTGATAAGACGCCGCAGAAAACAGTCCTAATACCAGTACGGCAAaataa      >  2:562228/1‑140 (MQ=255)                                                                                                                           tGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCTTTTGATAAGACGCCGCAGAAAACAGTCCTAATACCAGTACGGCAAAATAATATc  >  1:674521/1‑140 (MQ=255)                                                                                                                           tGCAGCAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCTTTTGATAAGACGCCGCAGAAAACAGTCCTAATACCAGTACGGCAAAATAATATc  >  2:644711/1‑140 (MQ=255)                                                                                                                               |                                                                                                                                        CCGGCATCACGAATATTTTTCTGCACCGCTACTGCACCAAACAAGCTTAAGAAGAAAGCCAGTCCATAAAAACCTTTTTCGCTGAGTAATAATGTCGCATTCCACAGACCTACCATCAGTAATGCAACAGAGATAATAAAGACAGTCAGGCAGGTCATATAATAAATGGAAGTGGTCGGTATGCCTTCATACTTGTCCCGAACGGTCTTTTGATAAGACGCCGCAGAAAACAGTCCTAATACCAGTACGGCAAAATAATATC  >  NC_000913/3726976‑3727237

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