Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,916,069	T→C	G20G (GGA→GGG)	msrC ←	free methionine‑(R)‑sulfoxide reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,916,069	0	T	C	100.0%	144.8 / NA	41	G20G (GGA→GGG)	msrC	free methionine‑(R)‑sulfoxide reductase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (14/27);  total (14/27)

TATCCGGACACAGGCAATTTTGCCCTGAAATGGTCCGAGTACCAGTGTATCGTCCTCAAGCAAATAAAAACCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCTCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCTTACCTGGATCTGCCAGCCTATTAAAATAAGCATTAAATGCGTTAATGCTCAAGATCATTCCCATCA  >  NC_000913/1915922‑1916215                                                                                                                                                    |                                                                                                                                                   tATCCGGACACAGGCAATTTTGCCCTGAAATGGTCCGAGTACCAGTGTATCGTCCTCAAGCAAATAAAAACCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCcccc                                                                                                                                                 >  2:106041/1‑151 (MQ=255)  aTCCGGACACAGGCAATTTTGCCCTGAAATGGTCCGAGTACCAGTGTATCGTCCTCAAGCAAATAAAAACCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCcccc                                                                                                                                                 >  1:453822/1‑150 (MQ=255)  aTCCGGACACAGGCAATTTTGCCCTGAAATGGTCCGAGTACCAGTGTATCGTCCTCAAGCAAATAAAAACCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCcccc                                                                                                                                                 <  2:9345/150‑1 (MQ=255)       gACACAGGCAATTTTGCCCTGAAATGGTCCGAGTACCAGTGTATCGTCCTCAAGCAAATAAAAACCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATc                                                                                                                                           <  1:269276/151‑1 (MQ=255)        acacAGGCAATTTTGCCCTGAAATGGTCCGAGTACCAGTGTATCGTCCTCAAGCAAATAAAACCCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCa                                                                                                                                          <  2:123357/151‑1 (MQ=255)                aaTTTTGCCCTGAAATGGTCCGAGTACCAGTGTATCGTCCTCAAGCAAATAAAAACCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTa                                                                                                                                  <  1:290674/151‑1 (MQ=255)                   tttGCCCTGAAATGGTCCGAGTACCAGTGTATCGTCCTCAAGCAAATAAAAACCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTaa                                                                                                                                 >  1:346409/1‑149 (MQ=255)                    ttGCCCTGAAATGGTCCGAGTACCAGTGTATCGTCCTCAAGCAAATAAAAACCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAg                                                                                                                               <  1:572241/150‑1 (MQ=255)                            aaaTGGTCCGAGTACCAGTGTATCGTCCTCAAGCAAATAAAAACCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGAttt                                                                                                                      >  1:546802/1‑151 (MQ=255)                                 gTCCGAGTACCAGTGTATCGTCCTCAAGCAAATAAAAACCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATc                                                                                                                 <  1:120903/151‑1 (MQ=255)                                    cGAGTACCAGTGTATCGTCCTCAAGCAAATAAAAACCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCcgc                                                                                                              >  1:489597/1‑151 (MQ=255)                                     gAGTACCAGTGTATCGTCCTCAAGCAAATAAAAACCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCcgcg                                                                                                             >  1:605826/1‑151 (MQ=255)                                        tACCAGTGTATCGTCCTCAAGCAAATAAAAACCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTa                                                                                                           >  2:52972/1‑150 (MQ=255)                                          ccAGTGTATCGTCCTCAAGCAAATAAAAACCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTaaa                                                                                                         >  1:275202/1‑150 (MQ=255)                                             gtgtATCGTCCTCAAGCAAATAAAAACCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAAtt                                                                                                      <  1:592797/150‑1 (MQ=255)                                                   cGTCCTCAAGCAAATAAAAACCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGtttt                                                                                               <  2:619973/151‑1 (MQ=255)                                                         cAAGCAAATAAAAACCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCAt                                                                                         <  1:84250/151‑1 (MQ=255)                                                           aGCAAATAAAAACCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATta                                                                                       <  2:571026/151‑1 (MQ=255)                                                           aGCAAATAAAAACCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATta                                                                                       <  1:51502/151‑1 (MQ=255)                                                                  aaaaaCCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAAt                                                                                 <  1:378182/150‑1 (MQ=255)                                                                    aaaCCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCAc                                                                              >  2:552830/1‑151 (MQ=255)                                                                          gCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTgg                                                                          <  2:67280/149‑1 (MQ=255)                                                                           cccAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTg                                                                       <  2:630968/151‑1 (MQ=255)                                                                           cccAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTg                                                                       <  1:667861/151‑1 (MQ=255)                                                                                tttATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCtt                                                                   >  2:404409/1‑150 (MQ=255)                                                                                 ttATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCTTAc                                                                 >  1:388550/1‑151 (MQ=255)                                                                                          tGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCTTACCTGGATCTg                                                        >  1:146440/1‑151 (MQ=255)                                                                                               cGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCTTACCTGGATCTGCCAg                                                    <  1:125421/150‑1 (MQ=255)                                                                                                 cTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCTTACCTGGATCTGCCAGCCt                                                 >  2:423977/1‑151 (MQ=255)                                                                                                   cATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCTTACCTGGATCTGCCAGCCTa                                                <  2:369700/150‑1 (MQ=255)                                                                                                   cATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCTTACCTGGATCTGCCAGCCTAt                                               <  1:649511/151‑1 (MQ=255)                                                                                                   cATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCTTACCTGGATCTGCCAGCCTAt                                               <  1:279897/151‑1 (MQ=255)                                                                                                    atatAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCTTACCTGGATCTGCCAGCCTAt                                               <  2:9188/150‑1 (MQ=255)                                                                                                      ataACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCTTACCTGGATCTGCCAGCCTATTa                                             <  1:409731/150‑1 (MQ=255)                                                                                                             cGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCTTACCTGGATCTGCCAGCCTATTAAAATAAGc                                     <  1:332344/151‑1 (MQ=255)                                                                                                              gCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCTTACCTGGATCTGCCAGCCTATTAAAATAAGCa                                    <  2:533143/151‑1 (MQ=255)                                                                                                               cACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCTTACCTGGATCTGCCAGCCTATTAAAATAAGCAt                                   <  2:557455/151‑1 (MQ=255)                                                                                                                           cAAGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCTTACCTGGATCTGCCAGCCTATTAAAATAAGCATTAAATGCGTTa                        >  2:321133/1‑150 (MQ=255)                                                                                                                            aaGCGTTGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCTTACCTGGATCTGCCAGCCTATTAAAATAAGCATTAAATGCGTTAAt                      >  1:321233/1‑151 (MQ=255)                                                                                                                                 ttGCCAGAAAACTGGTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCTTACCTGGATCTGCCAGCCTATTAAAATAAGCATTAAATGCGTTAATGCTc                  <  2:621703/150‑1 (MQ=255)                                                                                                                                     cAGAAAACTGGTTTCCCCCGCCATCAGCGCGTCAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCTTACCTGGATCTGCCAGCCTATTAAAATAAGCATTAAATGCGTTAATGCTCAAGAt             <  1:70259/151‑1 (MQ=255)                                                                                                                                              ggTTTCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCTTACCTGGATCTGCCAGCCTATTAAAATAAGCATTAAATGCGTTAATGCTCAAGATCATTccc      <  1:209745/149‑1 (MQ=255)                                                                                                                                                tttCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCTTACCTGGATCTGCCAGCCTATTAAAATAAGCATTAAATGCGTTAATGCTCAAGATCATTCCcatc   >  1:639767/1‑150 (MQ=255)                                                                                                                                                 ttCCCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCTTACCTGGATCTGCCAGCCTATTAAAATAAGCATTAAATGCGTTAATGCTCAAGATCATTCCcatca  <  1:523476/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   TATCCGGACACAGGCAATTTTGCCCTGAAATGGTCCGAGTACCAGTGTATCGTCCTCAAGCAAATAAAAACCTGCCCAGTTTATGTCAGTGAGACGCTCATATAACAACGCACTGGTGTTCGCAAGCGTTGCCAGAAAACTGGTTTCTCCCGCCATCAGCGCGTTAAAGTCGCGATTTAAATCCGCGTAAAATTCTGTTTTGTTCATTATATAATCACTTGGTTGTCTTACCTGGATCTGCCAGCCTATTAAAATAAGCATTAAATGCGTTAATGCTCAAGATCATTCCCATCA  >  NC_000913/1915922‑1916215

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