Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,212,703	C→T	*135* (TAG→TAA)	ycgX ←	DUF1398 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,212,703	0	C	T	100.0%	200.8 / NA	64	*135* (TAG→TAA)	ycgX	DUF1398 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (33/31);  total (33/31)

GATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGCTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGCCAGTGTTGAA  >  NC_000913/1212578‑1212832                                                                                                                              |                                                                                                                                  gATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGtt                                                                                                                      >  2:936452/1‑139 (MQ=255)  attattTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGttt                                                                                                                     >  1:405561/1‑139 (MQ=255)              tGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTAc                                                                                                         <  1:1160103/139‑1 (MQ=255)              atttgtataagaGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTAc                                                                                                         <  2:290707/128‑1 (MQ=255)                 cAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGt                                                                                                                       <  2:290133/122‑1 (MQ=255)                 cAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGt                                                                                                                       >  1:290133/1‑122 (MQ=255)                   aTAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCtt                                                                                                    <  2:405561/139‑1 (MQ=255)                            cAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTttatt                                                                                                                  <  2:496782/116‑1 (MQ=255)                            cAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTttatt                                                                                                                  >  1:496782/1‑116 (MQ=255)                              gTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCa                                                                                                  >  2:501439/1‑130 (MQ=255)                              gTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCa                                                                                                  <  1:501439/130‑1 (MQ=255)                                         ccATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTGAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATc                                                                              >  1:1282684/1‑139 (MQ=255)                                         ccATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTattatt                                                                                                               <  2:978580/106‑1 (MQ=255)                                         ccATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTattatt                                                                                                               >  1:978580/1‑106 (MQ=255)                                         ccATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATc                                                                              >  1:952061/1‑139 (MQ=255)                                         ccATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATc                                                                              >  2:1108830/1‑139 (MQ=255)                                         ccATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATc                                                                              >  2:306153/1‑139 (MQ=255)                                         ccATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATc                                                                              >  1:339128/1‑139 (MQ=255)                                         ccATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATc                                                                              >  1:175054/1‑139 (MQ=255)                                         ccATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATc                                                                              <  2:724810/139‑1 (MQ=255)                                          cATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCa                                                                             <  2:339128/139‑1 (MQ=255)                                            tAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAAt                                                                           <  1:306153/139‑1 (MQ=255)                                                 aTAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTc                                                                                                            <  1:237808/101‑1 (MQ=255)                                                 aTAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTc                                                                                                            >  2:237808/1‑101 (MQ=255)                                                   aGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTttatta                                                                                                                 >  2:1035104/1‑94 (MQ=255)                                                   aGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTttatta                                                                                                                 <  1:1035104/94‑1 (MQ=255)                                                   aGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCt                                                                                                           >  2:492740/1‑100 (MQ=255)                                                   aGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCt                                                                                                           <  1:492740/100‑1 (MQ=255)                                                   aGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCt                                                                                                           <  2:651885/100‑1 (MQ=255)                                                   aGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCt                                                                                                           >  1:651885/1‑100 (MQ=255)                                                   aGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCgtag                                                                    >  2:1079338/1‑139 (MQ=255)                                                           cACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCt                                                                                                           >  1:687536/1‑92 (MQ=255)                                                           cACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCt                                                                                                           <  2:687536/92‑1 (MQ=255)                                                                 tATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGtt                                                      <  1:1108830/139‑1 (MQ=255)                                                                          tATGCTAATGCCTTACTTAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATaaa                                                                                         <  1:1175945/95‑1 (MQ=255)                                                                          tATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATaaa                                                                                         >  2:1175945/1‑95 (MQ=255)                                                                                 aTGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCa                                                                             <  2:832455/100‑1 (MQ=255)                                                                                 aTGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCa                                                                             >  1:832455/1‑100 (MQ=255)                                                                                  tGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATAAATGGCGTAGTAAc                                                                <  1:1160915/112‑1 (MQ=255)                                                                                  tGCCTTACTGAATTATGAAGCATTTCTTAAGTATCAAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTaat                                                                >  2:1160915/1‑111 (MQ=255)                                                                                      ttACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCtt                                                                                                    <  2:562322/72‑1 (MQ=255)                                                                                      ttACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCtt                                                                                                    >  1:562322/1‑72 (MQ=255)                                                                                      ttACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCa                                                                                                  >  2:769777/1‑74 (MQ=255)                                                                                      ttACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCa                                                                                                  <  2:1171283/74‑1 (MQ=255)                                                                                      ttACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCa                                                                                                  >  1:1171283/1‑74 (MQ=255)                                                                                      ttACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCa                                                                                                  <  1:769777/74‑1 (MQ=255)                                                                                       tACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTa                                               >  1:1007218/1‑124 (MQ=255)                                                                                       tACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTa                                               <  2:1007218/124‑1 (MQ=255)                                                                                        aCTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCg                                                                       >  1:51896/1‑99 (MQ=255)                                                                                        aCTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCg                                                                       <  2:51896/99‑1 (MQ=255)                                                                                        aCTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCAcc                               >  1:638452/1‑139 (MQ=255)                                                                                                 tGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGc                      >  1:353487/1‑139 (MQ=255)                                                                                                      cATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCcgcg                 <  2:353487/139‑1 (MQ=255)                                                                                                           cTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTa                                      <  1:760976/113‑1 (MQ=255)                                                                                                           cTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTa                                      >  2:760976/1‑113 (MQ=255)                                                                                                                gTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTg                                                              >  2:1151028/1‑84 (MQ=255)                                                                                                                gTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTg                                                              <  1:1151028/84‑1 (MQ=255)                                                                                                                   tCCAACTTTAGTTAGATTAATTGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGc            >  2:809385/1‑131 (MQ=255)                                                                                                                   tCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTa                                               >  2:660232/1‑96 (MQ=255)                                                                                                                   tCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTa                                               <  1:660232/96‑1 (MQ=255)                                                                                                                   tCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGc            <  1:809385/131‑1 (MQ=255)                                                                                                                    ccAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGc            >  2:543292/1‑130 (MQ=255)                                                                                                                    ccAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGc            <  1:543292/130‑1 (MQ=255)                                                                                                                     cAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGCCAGTGTTGaa  <  2:1282684/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                  GATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGCTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGCCAGTGTTGAA  >  NC_000913/1212578‑1212832

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