Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,022,584	T→C	T434A (ACG→GCG)	yccS ←	putative transporter, FUSC superfamily inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,022,584	0	T	C	100.0%	157.2 / NA	45	T434A (ACG→GCG)	yccS	putative transporter, FUSC superfamily inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/23);  total (22/23)

ACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAGCCGAAGCACAACACCAGCGACATTCTTACCGCATGACGGAAGAGGGCGGAT  >  NC_000913/1022439‑1022723                                                                                                                                                  |                                                                                                                                            aCGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGt                                                                                                                                        <  1:620551/151‑1 (MQ=255) aCGTGACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACAAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGt                                                                                                                                        <  2:282330/151‑1 (MQ=255)   gTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTta                                                                                                                                      <  1:526727/151‑1 (MQ=255)      aCGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTtat                                                                                                                                     >  2:242878/1‑149 (MQ=255)          aaGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTtg                                                                                                                               <  1:90048/151‑1 (MQ=255)           aGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTtgg                                                                                                                              <  2:229690/151‑1 (MQ=255)              cAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTtggctg                                                                                                                           <  1:633103/151‑1 (MQ=255)               aaaaaaGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTtggctgg                                                                                                                          >  2:160842/1‑151 (MQ=255)               aaaaaaGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTtggctgg                                                                                                                          >  1:605159/1‑151 (MQ=255)                aaaaaGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGTCGGTGGCGCGCGGCGTTATAGTTtggct                                                                                                                            >  1:282597/1‑148 (MQ=255)                 aaaaGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTtggctggc                                                                                                                         >  2:521870/1‑150 (MQ=255)                 aaaaGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTtggctggc                                                                                                                         >  2:433632/1‑150 (MQ=255)                       gCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGacaaa                                                                                                                  <  1:531814/151‑1 (MQ=255)                         aCGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGacaaaca                                                                                                                <  1:493702/151‑1 (MQ=255)                          cGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGacaaaca                                                                                                                >  2:285212/1‑150 (MQ=255)                            ccGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACaaacaaac                                                                                                             >  1:556668/1‑151 (MQ=255)                                 aataaCCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTgcc                                                                                                        >  2:530542/1‑149 (MQ=255)                                        agcagcACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCa                                                                                                    >  2:663082/1‑148 (MQ=255)                                                       ccTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGAt                                                                                  <  2:530868/151‑1 (MQ=255)                                                            cAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCAtt                                                                             >  1:187690/1‑151 (MQ=255)                                                                  tGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGtt                                                                       >  2:378481/1‑151 (MQ=255)                                                                    gCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTAt                                                                     <  1:230872/151‑1 (MQ=255)                                                                        aaaCCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTg                                                                  >  2:596839/1‑150 (MQ=255)                                                                         aaCCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGaa                                                                <  2:437780/151‑1 (MQ=255)                                                                           ccACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAatg                                                              <  1:465023/151‑1 (MQ=255)                                                                                    agGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAg                                                     <  1:147424/151‑1 (MQ=255)                                                                                     ggAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAGc                                                    >  1:341372/1‑151 (MQ=255)                                                                                     ggAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAGc                                                    >  2:428/1‑151 (MQ=255)                                                                                       aaTGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAGCCg                                                  >  2:412713/1‑151 (MQ=255)                                                                                       aaTGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAGCCg                                                  <  2:512263/151‑1 (MQ=255)                                                                                       aaTGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAGCCg                                                  <  1:400442/151‑1 (MQ=255)                                                                                                 gCGATACCTACCAGCGTACCAATAATCCTTAACTGCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAGCCGAAGCacaaca                                        <  2:203564/151‑1 (MQ=255)                                                                                                    aTACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAGCCGAAGCACAACACCa                                     <  2:568711/151‑1 (MQ=255)                                                                                                      acctaccAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAGCCGAAGCACAACACCAGc                                   >  2:241990/1‑151 (MQ=255)                                                                                                        ctaccAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAGCCGAAGCACAACACCAGCg                                  >  1:93123/1‑150 (MQ=255)                                                                                                          accAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAGCCGAAGCACAACACCAGCGACa                               <  2:666595/151‑1 (MQ=255)                                                                                                            cAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAGCCGAAGCACAACACCAGCGACAtt                             <  2:393912/151‑1 (MQ=255)                                                                                                                  aCCAATAATCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAGCCGAAGCACAACACCAGCGACATTCTTAcc                       <  1:532599/151‑1 (MQ=255)                                                                                                                     aataatCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAGCCGAAGCACAACACCAGCGACATTCTTACCGCa                    <  1:164856/151‑1 (MQ=255)                                                                                                                        aatCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAGCCGAAGCACAACACCAGCGACATTCTTACCGCATg                  <  2:450330/150‑1 (MQ=255)                                                                                                                        aatCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAGCCGAAGCACAACACCAGCGACATTCTTACCGCATg                  >  1:28881/1‑150 (MQ=255)                                                                                                                        aatCCTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAGCCGAAGCACAACACCAGCGACATTCTTACCGCATg                  >  1:283045/1‑150 (MQ=255)                                                                                                                           ccTTAACTTCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAGCCGAAGCACAACACCAGCGACATTCTTACCGCATGACg               <  2:259273/150‑1 (MQ=255)                                                                                                                                  ttCAGGCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAGCCGAAGCACAACACCAGCGACATTCTTACCGCATGACGGAAGAggg       >  1:143657/1‑151 (MQ=255)                                                                                                                                       gCGGTGGCGCGCGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAGCCGAAGCACAACACCAGCGACATTCTTACCGCATGACGGAAGAGGGCGGAt  >  1:635902/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                            ACGTTACGGAAGGCAAAAAAGAGCACGCCGGTAATAACCAGCAGCACCAGCTGCCCTTCCAGTGATGGCACAAACCACAGCACAGGAATGCCAATGGCGATACCTACCAGCGTACCAATAATCCTTAACTTCAGGCGGTGGCGCGTGGCGTTATAGTTTGGCTGGCAGACAAACAAACTTGTCAGCAAGATCCAATACCCGTGATGCATTCCGGTTATCTGAATGATGGCGTAGCCGAAGCACAACACCAGCGACATTCTTACCGCATGACGGAAGAGGGCGGAT  >  NC_000913/1022439‑1022723

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