Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,652,429	T→C	intergenic (‑142/‑65)	yfhM ← / → sseA	alpha2‑macroglobulin/3‑mercaptopyruvate sulfurtransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,652,429	0	T	C	100.0%	180.3 / NA	50	intergenic (‑142/‑65)	yfhM/sseA	alpha2‑macroglobulin/3‑mercaptopyruvate sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/24);  total (26/24)

CTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGCCGACTGGCTCGCCGAACATATTGATGACCCGGAAATTCAGA  >  NC_000913/2652297‑2652560                                                                                                                                     |                                                                                                                                    cTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGt                                                                                                                              >  2:203621/1‑140 (MQ=255)   gCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCg                                                                                                                            <  1:637245/140‑1 (MQ=255)    ccACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGc                                                                                                                           <  1:894255/140‑1 (MQ=255)        aaTGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCtttt                                                                                                                       <  1:831095/140‑1 (MQ=255)        aaTGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCtttt                                                                                                                       <  1:203621/140‑1 (MQ=255)               cGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTaaa                                                                                                                >  1:88193/1‑140 (MQ=255)                       aCGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGAc                                                                                                        >  2:281756/1‑140 (MQ=255)                           cGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAAtt                                                                                                    <  1:551340/140‑1 (MQ=255)                                       atGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAAcccc                                                                                        >  2:839908/1‑140 (MQ=255)                                          cGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGc                                                                                     >  2:406747/1‑140 (MQ=255)                                          cGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGc                                                                                     >  1:572278/1‑140 (MQ=255)                                           gTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCa                                                                                    >  1:305287/1‑140 (MQ=255)                                            tGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCaa                                                                                   >  1:247618/1‑140 (MQ=255)                                            tGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCaa                                                                                   >  1:552733/1‑140 (MQ=255)                                                    tACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGgag                                                                           <  1:566868/140‑1 (MQ=255)                                                         tATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGgag                                                                           <  1:862995/135‑1 (MQ=255)                                                         tATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGgag                                                                           >  2:862995/1‑135 (MQ=255)                                                          aTGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCt                                                                     <  1:918572/140‑1 (MQ=255)                                                          aTGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCt                                                                     >  2:918572/1‑140 (MQ=255)                                                                     tATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGAc                                                          >  1:927964/1‑140 (MQ=255)                                                                     tATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGAc                                                          >  1:198985/1‑140 (MQ=255)                                                                     tATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGAc                                                          >  1:85922/1‑140 (MQ=255)                                                                     tATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGAc                                                          >  1:244787/1‑140 (MQ=255)                                                                     tATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGAc                                                          >  1:69115/1‑140 (MQ=255)                                                                      aTCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACa                                                         <  2:552733/140‑1 (MQ=255)                                                                      aTCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACa                                                         <  2:247618/140‑1 (MQ=255)                                                                        cACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATg                                                       >  1:614266/1‑140 (MQ=255)                                                                          cGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACa                                                         <  2:38020/136‑1 (MQ=255)                                                                          cGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACa                                                         >  1:38020/1‑136 (MQ=255)                                                                          cGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGt                                                     >  2:729444/1‑140 (MQ=255)                                                                            cccATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGttt                                                   >  2:445332/1‑140 (MQ=255)                                                                                     tCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGcc                                          >  1:155365/1‑140 (MQ=255)                                                                                     tCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGcc                                          >  2:703233/1‑140 (MQ=255)                                                                                     tCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGcc                                          >  2:380215/1‑140 (MQ=255)                                                                                       cGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGCCGa                                        <  2:473016/140‑1 (MQ=255)                                                                                       cGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGCCGa                                        >  1:851288/1‑140 (MQ=255)                                                                                         aaTTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGCCGACt                                      >  2:915499/1‑140 (MQ=255)                                                                                           ttttAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGTCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGCCGACTgg                                    <  2:88193/140‑1 (MQ=255)                                                                                            tttAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGCCGACTGGc                                   <  1:445332/140‑1 (MQ=255)                                                                                               aGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGCCGACTGACTCg                                <  1:839908/140‑1 (MQ=255)                                                                                                 tGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGCCGACTGGCTCGcc                              <  2:244948/140‑1 (MQ=255)                                                                                                 tGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGCCGACTGGCTCGcc                              <  1:915499/140‑1 (MQ=255)                                                                                                  gagaATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGCCGACTGGCTCGCCg                             <  2:244787/140‑1 (MQ=255)                                                                                                  gagaATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGCCGACTGGCTCGCCg                             <  2:85922/140‑1 (MQ=255)                                                                                                   agaATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGCCGACTGGCTCGCCGa                            <  2:155365/140‑1 (MQ=255)                                                                                                   agaATTGTGGCTCAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGCCGACTGGCTCGCCGa                            <  1:703233/140‑1 (MQ=255)                                                                                                               cAGACAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGCCGACTGGCTCGCCGAACATATtgatga                <  1:199681/140‑1 (MQ=255)                                                                                                                  aCAGCCAAATCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGCCGACTGGCTCGCCGAACATATTGATGAccc             >  1:273181/1‑140 (MQ=255)                                                                                                                         aaTCAGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGCCGACTGGCTCGCCGAACATATTGATGACCCGGAAAtt      <  2:69115/140‑1 (MQ=255)                                                                                                                             aGCCTCTTCTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATATCCACGACATGGTTTGTAGGAGCCGACTGGCTCGCCGAACATATTGATGACCCGGAAATTCAGa  <  2:305287/140‑1 (MQ=255)                                                                                                                                     |                                                                                                                                    CTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCCTAACCCCGGCAAAAATGGAGATGCCTATGTCCACGACATGGTTTGTAGGAGCCGACTGGCTCGCCGAACATATTGATGACCCGGAAATTCAGA  >  NC_000913/2652297‑2652560

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