Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,862,057	T→C	Y210H (TAT→CAT)	ygbJ →	putative dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,862,057	0	T	C	100.0%	172.9 / NA	44	Y210H (TAT→CAT)	ygbJ	putative dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (15/29);  total (15/29)

CCGGAAAAGTTTATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGTATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCT  >  NC_000913/2861911‑2862202                                                                                                                                                   |                                                                                                                                                  tcGGATAAGTTTATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAACGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGCGCGGGGATCCCGCTGGATGTGATGCAt                                                                                                                                                 >  1:441469/2‑149 (MQ=255)      aaaGTTTATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGAcgtcg                                                                                                                                          <  1:488498/151‑1 (MQ=255)         gattATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGGAAAAATTATTCACCAGTTGTAAGCGGGCGTACATAATGCTGCCGGAACCGAAGCGAAGGCACTTGCAGCCCGTGCGGGTAACCAGCTGGATGTGATGCATGATGTCGTg                                                                                                                                        >  1:34945/3‑150 (MQ=255)             aTCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCa                                                                                                                                    >  1:198391/1‑150 (MQ=255)                 cATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATgcc                                                                                                                               <  1:562418/151‑1 (MQ=255)                    aGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATgccgc                                                                                                                             <  1:111126/150‑1 (MQ=255)                          aGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAAt                                                                                                                      <  2:197145/151‑1 (MQ=255)                            aaCCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTc                                                                                                                    >  2:282621/1‑151 (MQ=255)                              ccGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCt                                                                                                                  <  2:198394/151‑1 (MQ=255)                              ccGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCt                                                                                                                  <  2:659805/151‑1 (MQ=255)                                gggACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTgg                                                                                                                <  2:464458/151‑1 (MQ=255)                                    cTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGGTGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGt                                                                                                            >  2:675868/1‑151 (MQ=255)                                      aGGTTCGACCGTAAAAATTATTCACCAGTTGTTACCGGGCGTACATATTGCTGCCGGAGACGAAGCCATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGCAATTCCTGGATGtt                                                                                                           >  1:75048/1‑150 (MQ=255)                                        gTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGa                                                                                                        <  2:137610/151‑1 (MQ=255)                                         ttCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGa                                                                                                        >  1:346960/1‑150 (MQ=255)                                                     aaTTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGt                                                                                           <  2:538003/151‑1 (MQ=255)                                                        tattCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCa                                                                                         <  2:301624/150‑1 (MQ=255)                                                            cACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATgtgg                                                                                    >  2:397025/1‑151 (MQ=255)                                                             aCCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATgtggt                                                                                   <  2:322508/151‑1 (MQ=255)                                                                 gttgttAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGAt                                                                               >  1:686965/1‑151 (MQ=255)                                                                 gttgttAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGAt                                                                               <  2:656630/151‑1 (MQ=255)                                                                  ttgttAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATg                                                                              >  1:43120/1‑151 (MQ=255)                                                                   tgttAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATg                                                                              <  2:63405/150‑1 (MQ=255)                                                                    gttAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATgg                                                                             <  2:172239/150‑1 (MQ=255)                                                                    gttAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGc                                                                            <  2:291821/151‑1 (MQ=255)                                                                          gggCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTa                                                                       >  1:113081/1‑150 (MQ=255)                                                                                aCATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACAcccc                                                                 <  2:318774/150‑1 (MQ=255)                                                                                  atatTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGc                                                               >  1:231499/1‑150 (MQ=255)                                                                                   tatTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCAt                                                             >  1:336926/1‑151 (MQ=255)                                                                                            cGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTc                                                    <  1:364394/151‑1 (MQ=255)                                                                                                     aGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATAttttt                                           <  1:603454/151‑1 (MQ=255)                                                                                                      gCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTg                                          <  1:461128/151‑1 (MQ=255)                                                                                                         aTGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTa                                       <  1:663338/151‑1 (MQ=255)                                                                                                           ggCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTa                                       >  2:85960/1‑149 (MQ=255)                                                                                                             cACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGa                                   <  2:583789/151‑1 (MQ=255)                                                                                                                       ccGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTg                         <  1:176925/151‑1 (MQ=255)                                                                                                                       ccGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTg                         <  2:352729/151‑1 (MQ=255)                                                                                                                        cGTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTgg                        <  2:355375/151‑1 (MQ=255)                                                                                                                         gTGCGGGGATCGCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCTATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATTTTGGTCTg                         <  1:449450/149‑1 (MQ=255)                                                                                                                         gTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGt                       <  1:99032/151‑1 (MQ=255)                                                                                                                         gTGCGGGGATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGt                       >  2:20993/1‑151 (MQ=255)                                                                                                                              gggATCCCGCTGGATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCg                  <  2:238359/151‑1 (MQ=255)                                                                                                                                         ggATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCaaa       <  2:596041/151‑1 (MQ=255)                                                                                                                                         ggATGTGATGCATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCa         <  2:687000/149‑1 (MQ=255)                                                                                                                                             gtgatgcATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGcc    <  2:63874/150‑1 (MQ=255)                                                                                                                                               gatgcATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCt  <  2:618749/150‑1 (MQ=255)                                                                                                                                               gatgcATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCt  <  2:367223/150‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  CCGGAAAAGTTTATCGCATAGGTGCAGAACCGGGACTAGGTTCGACCGTAAAAATTATTCACCAGTTGTTAGCGGGCGTACATATTGCTGCCGGAGCCGAAGCGATGGCACTTGCAGCCCGTGCGGGGATCCCGCTGGATGTGATGTATGACGTCGTGACCAATGCCGCCGGAAATTCCTGGATGTTCGAAAACCGGATGCGTCATGTGGTGGATGGCGATTACACCCCGCATTCAGCCGTCGATATTTTTGTTAAGGATCTTGGTCTGGTTGCCGATACAGCCAAAGCCCT  >  NC_000913/2861911‑2862202

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