Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,997,454	C→A	A121S (GCC→TCC)	tcyN ←	cystine ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,997,454	0	C	A	100.0%	199.5 / NA	59	A121S (GCC→TCC)	tcyN	cystine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (31/28);  total (31/28)

CCTCAGGACGCATTGCCAGCGCACGCGCAATCGCAACACGCTGCTGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGCCTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCTGGAAGACAAACCCGACGTGCTGACGCAACTGGCGAATCAGAGATTTTTG  >  NC_000913/1997322‑1997589                                                                                                                                     |                                                                                                                                        ccTCAGGACGCATTGCCAGCGCACGCGCAATCGCAACACGCTGCTGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACtcttct                                                                                                                                  <  1:1068702/140‑1 (MQ=255) ccTCAGGACGCATTGCCAGCGCACGCGCAATCGCAACACGCTGCTGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACtcttct                                                                                                                                  >  2:244250/1‑140 (MQ=255) ccTCAGGACGCATTGCCAGCGCACGCGCAATCGCAACACGCTGCTGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACtcttct                                                                                                                                  <  1:196869/140‑1 (MQ=255)            aTTGCCAGCGCACGCGCAATCGCAACACGCTGCTGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCAcc                                                                                                                       >  1:451239/1‑140 (MQ=255)            aTTGCCAGCGCACGCGCAATCGCAACACGCTGCTGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCAcc                                                                                                                       >  1:167817/1‑140 (MQ=255)            aTTGCCAGCGCACGCGCAATCGCAACACGCTGCTGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCAcc                                                                                                                       >  1:14626/1‑140 (MQ=255)                  aGCGCACGCGCAATCGCAACACGCTGCTGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGGCGTGGACTCTTCTTTCGGTTCACCTTTGAc                                                                                                                 >  2:722280/1‑140 (MQ=255)                    cgcACGCGCAATCGCAACACGCTGCTGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGa                                                                                                               >  1:1342572/1‑140 (MQ=255)                        cgcgCAATCGCAACACGCTGCTGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCa                                                                                                                         >  1:66649/1‑126 (MQ=255)                        cgcgCAATCGCAACACGCTGCTGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCa                                                                                                                         <  2:66649/126‑1 (MQ=255)                        cgcgCAATCGCAACACGCTGCTGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCAc                                                                                                           >  1:1071566/1‑140 (MQ=255)                         gcgcAATCGCAACACGCTGCTGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCAcc                                                                                                          <  1:336862/140‑1 (MQ=255)                         gcgcAATCGCAACACGCTGCTGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCAcc                                                                                                          <  1:1047506/140‑1 (MQ=255)                           gcAATCGCAACACGCTGCTGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTGCCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCgg                                                                                                        <  1:1094122/140‑1 (MQ=255)                              aTCGCAACACGCTGCTGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGccc                                                                                                     >  2:717611/1‑140 (MQ=255)                              aTCGCAACACGCTGCTGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGccc                                                                                                     <  2:1202679/140‑1 (MQ=255)                              aTCGCAACACGCTGCTGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGccc                                                                                                     <  2:1071566/140‑1 (MQ=255)                                       cGCTGCTGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCaataat                                                                                            <  1:1228146/140‑1 (MQ=255)                                       cGCTGCTGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCaataat                                                                                            >  1:281379/1‑140 (MQ=255)                                        gctgctGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATg                                                                                           <  2:344697/140‑1 (MQ=255)                                        gctgctGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATg                                                                                           <  2:1224604/140‑1 (MQ=255)                                        gctgctGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATg                                                                                           <  1:244250/140‑1 (MQ=255)                                         ctgctgTTGACCGCCAGACAAACGACGTGGATATCTGTTTTCTTTACCTACCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGAGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGt                                                                                          >  1:1135454/1‑140 (MQ=255)                                            ctgTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGccc                                                                                                     <  2:971035/126‑1 (MQ=255)                                            ctgTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGccc                                                                                                     >  1:971035/1‑126 (MQ=255)                                             tgTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTtctc                                                                                      <  2:451239/140‑1 (MQ=255)                                             tgTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCACCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTtctc                                                                                      <  2:441822/140‑1 (MQ=255)                                                  aCCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCCGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCa                                                                                 >  1:442106/1‑140 (MQ=255)                                                    cGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCaa                                                                                                >  1:358692/1‑123 (MQ=255)                                                    cGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCaa                                                                                                <  2:358692/123‑1 (MQ=255)                                                      ccAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGt                                                                             >  1:922717/1‑140 (MQ=255)                                                         gACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTtctc                                                                                      <  2:743083/128‑1 (MQ=255)                                                         gACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTtctc                                                                                      >  1:743083/1‑128 (MQ=255)                                                         gACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACg                                                                          >  1:1089251/1‑140 (MQ=255)                                                                  cgTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGaaa                                                                 <  2:1314307/140‑1 (MQ=255)                                                                  cgTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGaaa                                                                 <  2:281379/140‑1 (MQ=255)                                                                  cgTGGATAGCTGGTTTCTTTACCCGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGaaa                                                                 <  2:442106/140‑1 (MQ=255)                                                                            tGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAg                                                       <  1:1179424/140‑1 (MQ=255)                                                                            tGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAg                                                       <  1:926635/140‑1 (MQ=255)                                                                                ttCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCt                                                   >  1:670721/1‑140 (MQ=255)                                                                                     tACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCTGGAAg                                              >  1:1429351/1‑140 (MQ=255)                                                                                     tACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCTGGAAg                                              >  2:79033/1‑140 (MQ=255)                                                                                     tACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCTGGAAg                                              >  2:386216/1‑140 (MQ=255)                                                                                        cTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCTGGAAGACa                                           >  2:944632/1‑140 (MQ=255)                                                                                          gCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCTGGAAGACaaa                                         >  1:548906/1‑140 (MQ=255)                                                                                            cAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCTGGAAGACAAAcc                                       >  2:935398/1‑140 (MQ=255)                                                                                            cAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCTGGAAGACAAAcc                                       >  2:1404282/1‑140 (MQ=255)                                                                                            cAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCTGGAAGACAAAcc                                       <  2:670721/140‑1 (MQ=255)                                                                                               cccAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCTGGAAGACAAACCCGa                                    <  1:935398/140‑1 (MQ=255)                                                                                               cccAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCTGGAAGACAAACCCGa                                    <  1:386216/140‑1 (MQ=255)                                                                                                         gCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCTGGAAGACAAACCCGACGTg                                >  1:977503/1‑134 (MQ=255)                                                                                                         gCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCTGGAAGACAAACCCGACGTg                                <  2:977503/134‑1 (MQ=255)                                                                                                         gCCAGCAGCTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCTGGAAGACAAACCCGACGTGCTGACg                          >  1:1097499/1‑140 (MQ=255)                                                                                                              cagcTCGCGAGCGCGCGCCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCTGGAAGACAAACCCGACGTGCTGACGCAACt                     >  1:580420/1‑140 (MQ=255)                                                                                                                        gcgcgcgcCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCTGGAAGACAAACCCGACGTGCTGACGCAACTGGCGAATCag           >  1:503099/1‑140 (MQ=255)                                                                                                                           cgcgcCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCTGGAAGACAAACCCGACGTGCTGACGCAACTGGCGAATCAgaga        >  2:1410397/1‑140 (MQ=255)                                                                                                                           cgcgcCGTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCTGGAAGACAAACCCGACGTGCTGACGCAACTGGCGAATCAgaga        >  2:1228018/1‑140 (MQ=255)                                                                                                                                 gTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCTGGAAGACAAACCCGACGTGCTGACGCAACTGGCGAATCAGAGATTTTTg  <  2:1429351/140‑1 (MQ=255)                                                                                                                                 gTGGACTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCTGGAAGACAAACCCGACGTGCTGACGCAACTGGCGAATCAGAGATTTTTg  <  2:1127514/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                        CCTCAGGACGCATTGCCAGCGCACGCGCAATCGCAACACGCTGCTGTTGACCGCCAGACAAACGACGTGGATAGCTGGTTTCTTTACCTGCCAGCCCAACTTTTGCCAGCAGCTCGCGAGCGCGCGCCGTGGCCTCTTCTTTCGGTTCACCTTTGACGATCACCGGCCCTTCAATAATGTTCTCCAGCACCGTACGATGCGGAAACAAATTAAAGTTCTGGAAGACAAACCCGACGTGCTGACGCAACTGGCGAATCAGAGATTTTTG  >  NC_000913/1997322‑1997589

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