Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,441,538	C→T	*213* (TAG→TAA)	msrA ←	methionine sulfoxide reductase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,441,538	0	C	T	100.0%	219.4 / NA	67	*213* (TAG→TAA)	msrA	methionine sulfoxide reductase A
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (34/33);  total (34/33)

CCGATCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGT  >  NC_000913/4441408‑4441671                                                                                                                                   |                                                                                                                                      ccGATCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTc                                                                                                                               <  2:495275/139‑1 (MQ=255)    aTCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCcgg                                                                                                                            >  1:330000/1‑139 (MQ=255)    aTCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGAACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCcgg                                                                                                                            >  1:153475/1‑139 (MQ=255)         gCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGca                                                                                                                       >  1:245065/1‑139 (MQ=255)                   tCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcagaca                                                                                                                   <  1:860852/133‑1 (MQ=255)                   tCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcagaca                                                                                                                   >  2:860852/1‑133 (MQ=255)                   tCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgcc                                                                                                             >  1:1231341/1‑139 (MQ=255)                   tCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgcc                                                                                                             >  1:588943/1‑139 (MQ=255)                   tCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgcc                                                                                                             >  1:1456871/1‑139 (MQ=255)                           atatatTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                                     <  2:1456871/139‑1 (MQ=255)                             atatTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCATTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGAAGACAGAcgccaattccgcca                                                                                                   >  1:854592/1‑139 (MQ=255)                               atTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaat                                                                                                 <  2:588943/139‑1 (MQ=255)                                        ggACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgcc                                                                                                    <  1:762934/127‑1 (MQ=255)                                        ggACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgcc                                                                                                    >  2:762934/1‑127 (MQ=255)                                        ggACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTa                                                                                        >  1:683833/1‑139 (MQ=255)                                        ggACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTa                                                                                        >  2:92121/1‑139 (MQ=255)                                                ttGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACg                                                                                <  2:330000/139‑1 (MQ=255)                                                        cTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaatt                                                                                                <  2:1011850/115‑1 (MQ=255)                                                        cTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaatt                                                                                                >  1:1011850/1‑115 (MQ=255)                                                        cTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTa                                                                        >  1:258574/1‑139 (MQ=255)                                                         tGGTAAGTCGTTTACGCCACATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaatt                                                                                                <  1:343960/114‑1 (MQ=255)                                                         tGGTAAGTCGTTTACGCCACATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgccaatt                                                                                                >  2:343960/1‑114 (MQ=255)                                                          ggTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATg                                                                      <  2:140474/139‑1 (MQ=255)                                                          ggTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATg                                                                      <  1:1366857/139‑1 (MQ=255)                                                           gTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGc                                                                     >  2:189432/1‑139 (MQ=255)                                                           gTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGc                                                                     <  2:683833/139‑1 (MQ=255)                                                            tAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCa                                                                    <  1:1170484/139‑1 (MQ=255)                                                               gTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGat                                                                 <  2:1317995/139‑1 (MQ=255)                                                                  gTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATAtt                                                              <  2:245065/139‑1 (MQ=255)                                                                     tACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCt                                                           >  2:454401/1‑139 (MQ=255)                                                                     tACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCt                                                           >  2:566517/1‑139 (MQ=255)                                                                     tACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCt                                                           <  1:112641/139‑1 (MQ=255)                                                                     tACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCt                                                           >  1:1412256/1‑139 (MQ=255)                                                                                    aTAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCgg                                            <  1:92121/139‑1 (MQ=255)                                                                                    aTAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCgg                                            <  2:258574/139‑1 (MQ=255)                                                                                    aTAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCgg                                            <  1:454401/139‑1 (MQ=255)                                                                                            agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcaga                                                                                                                     <  1:901696/58‑1 (MQ=255)                                                                                            agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGcaga                                                                                                                     >  2:901696/1‑58 (MQ=255)                                                                                            agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                                     <  2:303935/74‑1 (MQ=255)                                                                                            agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgc                                                                                                     >  1:303935/1‑74 (MQ=255)                                                                                            agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTCACCATACGGGTTTTTATGCAGATATTGCCGGTGGTCATCTTCGGCAtaataa                                    >  1:1339181/1‑139 (MQ=255)                                                                                            agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCgg                                            >  1:245546/1‑131 (MQ=255)                                                                                            agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCgg                                            <  2:245546/131‑1 (MQ=255)                                                                                            agagCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCAtaataa                                    >  1:1297892/1‑139 (MQ=255)                                                                                              agCACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATaaaa                                  <  2:1339181/139‑1 (MQ=255)                                                                                                cACTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACg                                >  1:699686/1‑139 (MQ=255)                                                                                                  cTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATa                                                                                  >  1:847842/1‑87 (MQ=255)                                                                                                  cTTTATCAACAAACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATa                                                                                  <  2:847842/87‑1 (MQ=255)                                                                                                        caacaaACTACATTTGTACCCGTAACGTTATGCTTCcg                                                                                                                             >  1:503827/1‑38 (MQ=255)                                                                                                        caacaaACTACATTTGTACCCGTAACGTTATGCTTCcg                                                                                                                             <  2:503827/38‑1 (MQ=255)                                                                                                        caacaaACTACATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCg                        <  2:699686/139‑1 (MQ=255)                                                                                                                  cATTTGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGAttt              <  2:1231341/139‑1 (MQ=255)                                                                                                                     ttGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCg           <  1:566517/139‑1 (MQ=255)                                                                                                                     ttGTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCg           <  1:1318101/139‑1 (MQ=255)                                                                                                                       gTACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgtg         >  2:387995/1‑139 (MQ=255)                                                                                                                        tACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTAt                                                                       <  2:1405965/76‑1 (MQ=255)                                                                                                                        tACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTAt                                                                       >  1:1405965/1‑76 (MQ=255)                                                                                                                        tACCCGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCgtgg        <  1:387995/139‑1 (MQ=255)                                                                                                                          cccGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATaaaa                                  <  2:263336/111‑1 (MQ=255)                                                                                                                          cccGTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATaaaa                                  >  1:263336/1‑111 (MQ=255)                                                                                                                             gTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgcc                                                                                                    <  1:676104/42‑1 (MQ=255)                                                                                                                             gTAACGTTATGCTTCCGGCGGCAGACAGAcgccaattccgcc                                                                                                    >  2:676104/1‑42 (MQ=255)                                                                                                                             gTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGAtg   >  2:682421/1‑139 (MQ=255)                                                                                                                             gTAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGAtg   >  1:737268/1‑139 (MQ=255)                                                                                                                              tAACGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGAtgt  >  2:104063/1‑139 (MQ=255)                                                                                                                               aaCGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCt                                                           >  2:1199272/1‑81 (MQ=255)                                                                                                                               aaCGTTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCt                                                           <  1:1199272/81‑1 (MQ=255)                                                                                                                                   |                                                                                                                                      CCGATCAGGCCTACGCAGTCCTTGCAATATATTGAATTTGGACGATTTTGTAGGCCTGGTAAGTCGTTTACGCCGCATCAGGCATAGACAAAGAGCACTTTATCAACAAACTACATTTGTACCCGTAACGCTATGCTTCCGGCGGCAGACAGACGCCAATTCCGCCAATTCCACAGTAACCATACGGGTTTTTATGCAGATATTGCTGGTGGTCATCTTCGGCATAATAAAACGGTGTGGCGTTAGCGATTTCCGTGGTGATGT  >  NC_000913/4441408‑4441671

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