Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	477,025	G→A	*118* (TAG→TAA)	ybaA →	DUF1428 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	477,025	0	G	A	100.0%	219.7 / NA	67	*118* (TAG→TAA)	ybaA	DUF1428 family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/32);  total (35/32)

TTCAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAGGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGCACGCCATGCT  >  NC_000913/476893‑477155                                                                                                                                     |                                                                                                                                   ttCAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCg                                                                                                                              >  2:861559/1‑139 (MQ=255)  tCAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGc                                                                                                                               <  2:97706/137‑1 (MQ=255)  tCAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGc                                                                                                                               >  1:97706/1‑137 (MQ=255)  tCAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGt                                                                                                                             >  1:986799/1‑139 (MQ=255)    aaaaGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGa                                                                                                                           >  2:841448/1‑139 (MQ=255)        gagGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACggg                                                                                                                       <  2:2040/139‑1 (MQ=255)             ccGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCtgtg                                                                                                                  <  1:847677/139‑1 (MQ=255)                cgACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTgcgg                                                                                                               >  1:1480433/1‑139 (MQ=255)                   cGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTa                                                                                                                                 >  2:1137627/1‑118 (MQ=255)                   cGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTa                                                                                                                                 <  1:1137627/118‑1 (MQ=255)                   cGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGt                                                                                                            <  2:594467/139‑1 (MQ=255)                                  gatgatGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCaaaa                                                                                             >  1:1021134/1‑139 (MQ=255)                                   atgatgTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCaaaaa                                                                                            >  2:1434750/1‑139 (MQ=255)                                       tgTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACggg                                                                                                                       <  2:1168727/108‑1 (MQ=255)                                       tgTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACggg                                                                                                                       >  1:1168727/1‑108 (MQ=255)                                          cGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACaa                                                                                     >  2:641129/1‑139 (MQ=255)                                            gACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTca                                                                                        <  2:870982/134‑1 (MQ=255)                                            gACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTca                                                                                        >  1:870982/1‑134 (MQ=255)                                            gACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAAt                                                                                   >  1:892074/1‑139 (MQ=255)                                            gACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGAC‑‑GCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGt                                                                                 >  2:1597793/1‑139 (MQ=255)                                                 aCGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGc                                                                              <  2:1480433/139‑1 (MQ=255)                                                  cGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCt                                                                             <  1:1434750/139‑1 (MQ=255)                                                  cGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCt                                                                             <  2:892074/139‑1 (MQ=255)                                                  cGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCt                                                                             >  1:290645/1‑139 (MQ=255)                                                   ggATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTc                                                                            <  2:491571/139‑1 (MQ=255)                                                       gAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                        <  2:290645/139‑1 (MQ=255)                                                       gAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                        <  1:275421/139‑1 (MQ=255)                                                            agTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCg                                                                   >  2:194225/1‑139 (MQ=255)                                                              ttCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTACTCGGCCCAGCGTa                                                                 >  2:569426/1‑139 (MQ=255)                                                               tCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGt                                                                                                       >  2:13323/1‑100 (MQ=255)                                                               tCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGt                                                                                                       <  1:13323/100‑1 (MQ=255)                                                                cGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGTCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTgcggcg                                                                                                             <  2:1524803/93‑1 (MQ=255)                                                                cGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTgcggcg                                                                                                             >  1:1524803/1‑93 (MQ=255)                                                                      gTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTATGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAAtct                                                         >  1:1197135/1‑139 (MQ=255)                                                                       tCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGAC‑‑GCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAAtcttct                                                      <  1:1597793/139‑1 (MQ=255)                                                                         cATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAAtcttct                                                      <  2:986799/139‑1 (MQ=255)                                                                            gCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTc                                                   >  1:1309301/1‑139 (MQ=255)                                                                                                 gATCTATGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGagc                              <  1:841448/139‑1 (MQ=255)                                                                                                      aTGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTgc                                                                                                                 >  2:147082/1‑51 (MQ=255)                                                                                                      aTGGCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTgc                                                                                                                 <  1:147082/51‑1 (MQ=255)                                                                                                        ggCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                        >  1:726822/1‑90 (MQ=255)                                                                                                        ggCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                        <  2:726822/90‑1 (MQ=255)                                                                                                         gCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACg                                                                                                                         >  1:476958/1‑40 (MQ=255)                                                                                                         gCGGATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTg                      >  1:1179864/1‑139 (MQ=255)                                                                                                         gCGGATACGAGTCAATCATCGACGAATAAGTAGCGTGACg                                                                                                                         <  2:476958/40‑1 (MQ=255)                                                                                                            gATTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGAcc                   >  1:43024/1‑139 (MQ=255)                                                                                                             aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCaa                                                                                               <  1:1009439/62‑1 (MQ=255)                                                                                                             aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCaa                                                                                               >  2:1009439/1‑62 (MQ=255)                                                                                                             aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCt                                                                             >  2:466095/1‑80 (MQ=255)                                                                                                             aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCt                                                                             >  1:1003131/1‑80 (MQ=255)                                                                                                             aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCt                                                                             >  2:974778/1‑80 (MQ=255)                                                                                                             aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCt                                                                             <  1:466095/80‑1 (MQ=255)                                                                                                             aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCt                                                                             <  1:974778/80‑1 (MQ=255)                                                                                                             aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                        <  1:267236/85‑1 (MQ=255)                                                                                                             aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGcc                                                                        >  2:267236/1‑85 (MQ=255)                                                                                                             aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAAtct                                                         >  2:447059/1‑100 (MQ=255)                                                                                                             aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAAtct                                                         <  1:447059/100‑1 (MQ=255)                                                                                                             aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCt                                         >  1:207096/1‑116 (MQ=255)                                                                                                             aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCt                                         <  2:207096/116‑1 (MQ=255)                                                                                                             aTTCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCg                  >  1:718574/1‑139 (MQ=255)                                                                                                             aTTCGAGTCAAGCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCt                                                                             <  2:1003131/80‑1 (MQ=255)                                                                                                              ttCGAGTCAATCATCGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGt                 >  2:586602/1‑139 (MQ=255)                                                                                                                           tcGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGCACGCCATg    <  2:718574/139‑1 (MQ=255)                                                                                                                           tcGACGAATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGCACGCCATg    <  1:586602/139‑1 (MQ=255)                                                                                                                             gacgaATAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGCACGCCATGCt  <  1:861559/139‑1 (MQ=255)                                                                                                                                  aTAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCttt                                                    <  2:324243/84‑1 (MQ=255)                                                                                                                                  aTAAGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCttt                                                    >  1:324243/1‑84 (MQ=255)                                                                                                                                     |                                                                                                                                   TTCAAAAGAGGTCCGCGACGCTGCTAATCAAAAGATGATGTCGGACCCACGGATGAAAGAGTTCGGCGAGTCCATGCCGTTTGACGGCAAGCGAATGATCTATGGCGGATTCGAGTCAATCATCGACGAATAGGTAGCGTGACGGGCTGTGCGGCGTGCCGTCAGCCCAAAAATTCACAAATGTTGCTCGGCCCAGCGTAAGAAATCTTCTTTCGGTAATGCCTTGCTGTAGAGCCAGCCCTGACCGTAGTGCACGCCATGCT  >  NC_000913/476893‑477155

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