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%	141.4 / NA	44	*135* (TAG→TAA)	ycgX	DUF1398 family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/19);  total (25/19)

GTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGCTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGCCAGTGTTGAACAGTAT  >  NC_000913/1212576‑1212838                                                                                                                                |                                                                                                                                        gTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATATTCTCTGCACTTTTATTTCCATTATGCTATTGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATgg                                                                                                                              >  2:716612/1‑139 (MQ=255) gTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATgg                                                                                                                              >  1:1632658/1‑139 (MQ=255) gTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATgg                                                                                                                              >  1:702887/1‑139 (MQ=255)           atGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTAttt                                                                                                                    >  1:1677735/1‑139 (MQ=255)              ggTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGATTTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCt                                                                                                                 <  1:716612/139‑1 (MQ=255)               gTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTa                                                                                                                <  2:1632658/139‑1 (MQ=255)                     aTAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTACACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCtt                                                                                                          <  1:1486399/139‑1 (MQ=255)                              cAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTAGTATTTTCTACATCTTCAAtatata                                                                                                 >  2:1050648/1‑139 (MQ=255)                                gTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATaaa                                                                                               <  2:702887/139‑1 (MQ=255)                                gTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATaaa                                                                                               <  1:1275476/139‑1 (MQ=255)                                 tATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCtt                                                                                                          <  2:462404/127‑1 (MQ=255)                                 tATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCtt                                                                                                          >  1:462404/1‑127 (MQ=255)                                    tAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGc                                                                                            <  1:1210007/138‑1 (MQ=255)                                    tAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGc                                                                                            >  2:1210007/1‑138 (MQ=255)                                           ccATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATc                                                                                    >  2:1445687/1‑139 (MQ=255)                                                     aGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATaaaa                                                                                              >  2:1882787/1‑119 (MQ=255)                                                     aGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATaaaa                                                                                              <  1:1882787/119‑1 (MQ=255)                                                     aGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCgtag                                                                          >  1:779959/1‑139 (MQ=255)                                                      gTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCgtagt                                                                         <  2:1606377/139‑1 (MQ=255)                                                        ctctGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTaa                                                                       <  1:1050648/139‑1 (MQ=255)                                                           tGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAata                                                                                                     <  2:1871834/106‑1 (MQ=255)                                                           tGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAata                                                                                                     >  1:1871834/1‑106 (MQ=255)                                                                           ttATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTaa                                                    <  2:779959/139‑1 (MQ=255)                                                                              tGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACAt                                                 >  1:390491/1‑139 (MQ=255)                                                                                   aTGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGta                                                                                         <  2:671795/94‑1 (MQ=255)                                                                                   aTGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGta                                                                                         >  1:671795/1‑94 (MQ=255)                                                                                   aTGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTa                                            >  2:1365964/1‑139 (MQ=255)                                                                                    tGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAAc                                                   >  2:762929/1‑131 (MQ=255)                                                                                    tGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAAc                                                   <  1:762929/131‑1 (MQ=255)                                                                                         tACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTtta                                                                                                                          >  1:124325/1‑55 (MQ=255)                                                                                         tACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTtta                                                                                                                          <  2:124325/55‑1 (MQ=255)                                                                                          aCTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAATCCAcc                                     >  1:711814/1‑139 (MQ=255)                                                                                          aCTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATATTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCAcc                                     >  1:1098141/1‑139 (MQ=255)                                                                                              aaTTATGAAGCATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGaa                                 >  1:62910/1‑139 (MQ=255)                                                                                                        cATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGta                                                                                         <  2:976009/73‑1 (MQ=255)                                                                                                        cATTTCTTAAGTATCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGta                                                                                         >  1:976009/1‑73 (MQ=255)                                                                                                                   tATCCAACTTTAGTTAAATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCgtagt                                                                         <  1:925630/78‑1 (MQ=255)                                                                                                                   tATCCAACTTTAGTTAAATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCgtagt                                                                         >  2:925630/1‑78 (MQ=255)                                                                                                                    aTCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGCCAGTGtt           <  2:390491/139‑1 (MQ=255)                                                                                                                     tCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGCCAGTGTTg          >  2:1649862/1‑139 (MQ=255)                                                                                                                     tCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGCCAGTGTTg          >  2:1224416/1‑139 (MQ=255)                                                                                                                     tCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGCCAGTGTTg          <  2:711814/139‑1 (MQ=255)                                                                                                                     tCCAACTTTAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATCAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGCCAGTGTTg          >  1:48507/1‑139 (MQ=255)                                                                                                                             tAGTTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGTCAGTGTTGAACAGTAt  >  1:593437/1‑139 (MQ=255)                                                                                                                                |                                                                                                                                        GTGATTATTTATGGGTGACAATAAAAAGACAGTATTAATCATCCATAGAGATAGTCTCTGCACTTTTATTTCCATTATGCTAATGCCTTACTGAATTATGAAGCATTTCTTAAGTATCCAACTTTAGCTAGATTAATGGTTTATTATTTTCTACATCTTCAATATATAAAAGCGTATTATCAATGGCGTAGTAACTGCGTTTGTTATGATTAACATCAGTAACCCACCGGAAAACGCCCGCGCCTGCCAGTGTTGAACAGTAT  >  NC_000913/1212576‑1212838

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