Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	888,480	A→C	F447C (TTC→TGC)	ybjL ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	888,480	0	A	C	100.0%	185.9 / NA	51	F447C (TTC→TGC)	ybjL	putative transporter
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (22/29);  total (22/29)

AACAAGAAACAGATAACCACGGGCACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGAACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGCGTTCCCCATGCCGAAACTGAATGTGCTGAACTGGAA  >  NC_000913/888342‑888619                                                                                                                                           |                                                                                                                                            aaCAAGAAACAGATACCCACGGGCCCCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCAc                                                                                                                                   <  2:267273/149‑1 (MQ=255) aaCAAGAAACAGATACCCACGGGCACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCa                                                                                                                                 <  2:905560/151‑1 (MQ=255)    agaaaaaaaaTACCCACGGCCCCCAGCTGGACAATCAACCCGTAAATCAACACCTTGCCGCCAATCGCGCCCAGCCCGTTATTAATACCGCTCCCGTCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCa                                                                                                                                 <  1:327773/139‑1 (MQ=255)    aaaaatcAGATACCCGCGGGCTCCAGACTGACAATCAATCCGGCAATCACCATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCAGCCATAAACACCATCAAGCCGCACTCTTTCACCATGc                                                                                                                              <  2:786039/145‑1 (MQ=255)       aacaaaaaaacccggggCCCCAGCCGCACAATCAATCCGCCAATCAACTTCGGCCCGCCAATCGCCCCCAGGCCGTTATTAATACCGCTCCCGGCCCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTa                                                                                                                           <  1:476996/137‑1 (MQ=255)         aCAGATAACCACGGGCACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAAt                                                                                                                         <  2:663332/151‑1 (MQ=255)               acccACGGGCACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCAcc                                                                                                                    <  2:346319/148‑1 (MQ=255)                  cACGGGCACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGc                                                                                                                <  2:474122/151‑1 (MQ=255)                        cACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATg                                                                                                          >  1:158318/1‑151 (MQ=255)                          cctgacgaaCAATCACTCCGCCATTCACCCTCGCCCCGCCAATCCCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTAAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGaa                                                                                                            <  1:836180/139‑1 (MQ=255)                             gACTGACAACAATTACGACACACAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCACCGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAAcc                                                                                                     <  1:91420/151‑1 (MQ=255)                                      tcaatcCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCgg                                                                                            >  2:372345/1‑151 (MQ=255)                                        aatcCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGggt                                                                                          <  2:824186/151‑1 (MQ=255)                                         acccGGAAATCAACACCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGggtg                                                                                         <  2:364849/149‑1 (MQ=255)                                           cgggCAACCAACAGCTGCCCGCCAATCGCCCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGggt                                                                                          <  1:441217/146‑1 (MQ=255)                                                aaTCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCa                                                                                  <  2:910879/151‑1 (MQ=255)                                                    aaCACCTGCCCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTCCCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCAc                                                                                 <  1:468265/148‑1 (MQ=255)                                                             ccGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGccc                                                                      <  1:926925/150‑1 (MQ=255)                                                              cGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGccc                                                                      >  1:151693/1‑149 (MQ=255)                                                                ccAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGc                                                                   >  2:683234/1‑150 (MQ=255)                                                                 cAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCAt                                                                 >  1:590540/1‑151 (MQ=255)                                                                       cgcCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTcc                                                           <  2:322666/151‑1 (MQ=255)                                                                        gcCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCg                                                          >  1:658542/1‑151 (MQ=255)                                                                           cAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCg                                                       <  2:283855/151‑1 (MQ=255)                                                                                cGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATaa                                                  <  2:416519/151‑1 (MQ=255)                                                                                  ttattaATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATaac                                                 >  1:254253/1‑150 (MQ=255)                                                                                     ttaATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAAcc                                             >  1:579887/1‑151 (MQ=255)                                                                                       aaTACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCg                                           <  2:300728/151‑1 (MQ=255)                                                                                          aCCGCTCCCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAAccc                                            <  2:822775/147‑1 (MQ=255)                                                                                           ccGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGc                                         <  2:913163/149‑1 (MQ=255)                                                                                            cGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGc                                      >  2:19629/1‑151 (MQ=255)                                                                                              cTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGCGt                                    >  1:287401/1‑151 (MQ=255)                                                                                                 ccGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGCGt                                    <  1:457493/148‑1 (MQ=255)                                                                                                 ccGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGCGTTcc                                 >  2:676565/1‑151 (MQ=255)                                                                                                     cgcTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGCGTTCCCCAt                             <  2:171209/151‑1 (MQ=255)                                                                                                      gcTCAGCCCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGCGTTCCCCAt                             <  2:554443/150‑1 (MQ=255)                                                                                                      gcTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGCGTTCCCCATg                            >  1:336742/1‑151 (MQ=255)                                                                                                       cTCAGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGCGTTCCCCATGc                           >  1:824262/1‑151 (MQ=255)                                                                                                          aGACCAACGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGCGTTCCCCATGCCg                         >  1:433407/1‑150 (MQ=255)                                                                                                             ccAACCCCTGCCTTAACCTCCATCAAGCCGCACTCTTTCGCCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACATCCCGGCAGCTTTCCCCATGc                           <  1:435113/145‑1 (MQ=255)                                                                                                               aaCGCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGCGTTCCCCATGCCGAAACt                    >  1:535977/1‑150 (MQ=255)                                                                                                                  gCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGCGTTCCCCATGCCGAAACTg                   <  2:526955/148‑1 (MQ=255)                                                                                                                  gCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGCGTTCCCCATGCCGAAACTGAAt                >  2:742255/1‑151 (MQ=255)                                                                                                                  gCCTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGCGTTCCCCATGCCGAAACTGAAt                >  1:211572/1‑151 (MQ=255)                                                                                                                    cTGCCATAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGCGTTCCCCATGCCGAAACTGAAtgt              >  2:469824/1‑151 (MQ=255)                                                                                                                          tAAACACCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGCGTTCCCCATGCCGAAACTGAATGTGCTGa         >  2:33398/1‑150 (MQ=255)                                                                                                                             acacCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGCGTTCCCCATGCCGAAACTGAATGTGCTGAACTg     <  2:569212/151‑1 (MQ=255)                                                                                                                              cacCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGCGTTCCCCATGCCGAAACTGAATGTGCTGAACTgg    >  2:548907/1‑151 (MQ=255)                                                                                                                               acCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGCGTTCCCCATGCCGAAACTGAATGTGCTGAACTgg    >  1:347385/1‑150 (MQ=255)                                                                                                                               acCATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGCGTTCCCCATGCCGAAACTGAATGTGCTGAACTgg    >  2:354697/1‑150 (MQ=255)                                                                                                                                ccATCAAGCCGCACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGCGTTCCCCATGCCGAAACTGAATGTGCTGAACTGGaa  <  1:928278/151‑1 (MQ=255)                                                                                                                                           |                                                                                                                                            AACAAGAAACAGATAACCACGGGCACCAGACTGACAATCAATCCGGCAATCAACATCTGACCGCCAATCGCGCCCAGGCCGTTATTAATACCGCTACCGGCGCTCAGACCAACGCCTGCCATAAACACCATCAAGCCGAACTCTTTCACCATGCTTAATGCACCCTGCGGAATGTAACCGAAGGTCGGGTGGTTAGCACGCATAAAGCCCAGCATAATTCCGGCGAATAACAACCCGGCAGCGTTCCCCATGCCGAAACTGAATGTGCTGAACTGGAA  >  NC_000913/888342‑888619

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