Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,579,335	G→A	T3I (ACC→ATC)	yddA ←	ABC transporter family protein YddA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,579,335	0	G	A	100.0%	202.3 / NA	60	T3I (ACC→ATC)	yddA	ABC transporter family protein YddA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/31);  total (29/31)

TAACCACGCCAAGGATCATGGCGAGGATAATAATGATTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGGTTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTATCATTTAACAATATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATG  >  NC_000913/1579204‑1579468                                                                                                                                    |                                                                                                                                      tAACCACGCCAAGGATCATGGCGAGGATAATAATGATTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCAtt                                                                                                                               >  1:147629/1‑140 (MQ=255) tAACCACGCCAAGGATCATGGCGAGGATAATAATGATTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCAtt                                                                                                                               >  2:163434/1‑140 (MQ=255)   aCCACGCCAAGGATCATGGCGAGGATAATAATGATTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAg                                                                                                                             <  2:311360/140‑1 (MQ=255)      aCGCCAAGGATCATGGCGAGGATAATAATGATTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGa                                                                                                                          >  2:430454/1‑140 (MQ=255)                  aTGGCGAGGATAATAATGATTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAattga                                                                                                                   >  2:551640/1‑135 (MQ=255)                  aTGGCGAGGATAATAATGATTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAattga                                                                                                                   <  1:551640/135‑1 (MQ=255)                  aTGGCGAGGATAATAATGATTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATtgtg                                                                                                              <  1:550825/140‑1 (MQ=255)                  aTGGCGAGGATAATAATGATTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATtgtg                                                                                                              <  2:266320/140‑1 (MQ=255)                     gCGAGGATAATAATGATTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCa                                                                                                           <  2:701134/140‑1 (MQ=255)                       gAGGATAATAATGATTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATc                                                                                                         <  2:541110/140‑1 (MQ=255)                        aGGATAATAATGATTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATcc                                                                                                        <  1:142313/140‑1 (MQ=255)                        aGGATAATAATGATTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATcc                                                                                                        <  2:523700/140‑1 (MQ=255)                        aGGATAATAATGATTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATc                                                                                                         >  2:511975/1‑139 (MQ=255)                        aGGATAATAATGATTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATc                                                                                                         <  1:511975/139‑1 (MQ=255)                             aataatGATTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGa                                                                                                   <  2:147629/140‑1 (MQ=255)                               taatGATTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGAtt                                                                                                 >  2:69138/1‑140 (MQ=255)                                  tGATTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATtgtg                                                                                                              <  2:260986/124‑1 (MQ=255)                                  tGATTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATtgtg                                                                                                              >  1:260986/1‑124 (MQ=255)                                  tGATTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATc                                                                                              <  1:163434/140‑1 (MQ=255)                                    aTTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCa                                                                                                       <  2:398940/129‑1 (MQ=255)                                    aTTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCa                                                                                                       >  1:398940/1‑129 (MQ=255)                                    aTTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTa                                                                                            >  2:470759/1‑140 (MQ=255)                                    aTTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTa                                                                                            >  1:55597/1‑140 (MQ=255)                                     ttAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTAt                                                                                           >  2:56346/1‑140 (MQ=255)                                     ttAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACCTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCa                                                                                                       <  2:559789/128‑1 (MQ=255)                                     ttAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACCTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCa                                                                                                       >  1:559789/1‑128 (MQ=255)                                        acaacaCTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTATGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGc                                                                                        <  1:99929/140‑1 (MQ=255)                                                      ttattaTTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACtgat                                                                          >  2:456309/1‑140 (MQ=255)                                                      ttattaTTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACtgat                                                                          <  2:480113/140‑1 (MQ=255)                                                         ttattCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGaa                                                                       >  2:629591/1‑140 (MQ=255)                                                           attCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACtgatg                                                                         <  1:525637/136‑1 (MQ=255)                                                           attCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACtgatg                                                                         >  2:525637/1‑136 (MQ=255)                                                              cTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAattta                                                                  >  1:443655/1‑140 (MQ=255)                                                                 aCGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTAc                                                                              >  1:506188/1‑125 (MQ=255)                                                                 aCGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTAc                                                                              <  2:506188/125‑1 (MQ=255)                                                                 aCGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAatttattt                                                               <  2:491151/140‑1 (MQ=255)                                                                    cAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAg                                                            <  2:604458/140‑1 (MQ=255)                                                                      gCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTa                                                          <  2:55597/140‑1 (MQ=255)                                                                       ccAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTAt                                                         >  2:123650/1‑140 (MQ=255)                                                                       ccAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTAt                                                         >  1:236348/1‑140 (MQ=255)                                                                       ccAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTAt                                                         >  1:479618/1‑140 (MQ=255)                                                                         aaaaGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTATCa                                                       <  1:456309/140‑1 (MQ=255)                                                                                cTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTATCATTTAACa                                                <  1:430454/140‑1 (MQ=255)                                                                                 ttCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTATCATTTAACaa                                               >  1:231543/1‑140 (MQ=255)                                                                                     acagacagAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTATCATTTAACAATATg                                           <  1:470759/140‑1 (MQ=255)                                                                                         acagaGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATGTATGGCAATGATTTACTGATGAAATTTATTTAAGTATCATTTAACAATATGAAat                                       >  2:351107/1‑140 (MQ=255)                                                                                         acagaGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTATCATTTAACAATATGAAat                                       >  1:243056/1‑140 (MQ=255)                                                                                           agagATACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTATCATTTAACAATATGAAatat                                     >  2:478397/1‑140 (MQ=255)                                                                                             agaTACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTATCATTTAACAATATGAAATATaa                                   >  1:694675/1‑140 (MQ=255)                                                                                               aTACTTTGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTATCATTTAACAATATGAAATATAAAg                                 >  2:507034/1‑140 (MQ=255)                                                                                                     tGCTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTATCATTTAACAATATGAAATATAAAGGAatat                           >  2:9373/1‑140 (MQ=255)                                                                                                       cTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATTAAATTTATTTAAGTATCATTTAACAATATGAAAtata                                    <  2:690281/129‑1 (MQ=255)                                                                                                       cTATCAACATACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATTAAATTTATTTAAGTATCATTTAACAATATGAAAtata                                    >  1:690281/1‑129 (MQ=255)                                                                                                               aTACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTATCATTTAACAATATGAAATATAAAGGAATATTTCATTTTCa                 <  1:9373/140‑1 (MQ=255)                                                                                                                tACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTATCATTTAACAATATGAAATATAAAGGAATATTTCATTTTCAg                <  2:243056/140‑1 (MQ=255)                                                                                                                tACGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTATCATTTAACAATATGAAATATAAAGGAATATTTCATTTTCAg                <  2:694675/140‑1 (MQ=255)                                                                                                                 aCGAAGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTATCATTTAACAATATGAAATATAAAGGAATATTTCATTTTCAGt               >  1:688721/1‑140 (MQ=255)                                                                                                                    aaGCGTAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTATCATTTAACAATATGAAATATAAAGGAATATTTCATTTTCAGtggt            <  2:634677/140‑1 (MQ=255)                                                                                                                         tAATGGGAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTATCATTTAACAATATGAAATATAAAGGAATATTTCATTTTCAGtggttatgg       <  1:478397/140‑1 (MQ=255)                                                                                                                              ggAATGATTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTATCATTTAACAATATGAAATATAAAGGAATATTTCATTTTCAGTGgttatggttatg  <  2:236348/140‑1 (MQ=255)                                                                                                                                    |                                                                                                                                      TAACCACGCCAAGGATCATGGCGAGGATAATAATGATTAACAACACTGAAGTTTTATTATTCTTACGCAGCCAAAAGGGCTTCAACAGACAGAGATACTTTGCTATCAACATACGAAGCGTAATGGGAATGGTTATCATTAGCGAAAATTGATTGTGCCATCCAATGATTATCTATGGCAATGATTTACTGATGAAATTTATTTAAGTATCATTTAACAATATGAAATATAAAGGAATATTTCATTTTCAGTGGTTATGGTTATG  >  NC_000913/1579204‑1579468

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