Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,618,281	Δ1 bp	coding (628/666 nt)	marC ←	UPF0056 family inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,618,281	0	T	.	100.0%	252.0 / NA	52	coding (628/666 nt)	marC	UPF0056 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (33/19);  total (33/19)

CGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATC  >  NC_000913/1618136‑1618423                                                                                                                                                  |                                                                                                                                               cgGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGttaa‑aataa                                                                                                                                           >  2:508392/1‑150 (MQ=255)           ttttGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACAccc                                                                                                                                >  2:779520/1‑151 (MQ=255)             ttGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCcat                                                                                                                              >  2:659736/1‑151 (MQ=255)               gCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCcatac                                                                                                                            <  2:527571/151‑1 (MQ=255)                 cgcgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCcatacat                                                                                                                          <  1:628658/151‑1 (MQ=255)                   cgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCatacatac                                                                                                                        >  2:717929/1‑151 (MQ=255)                       aaTTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACcagc                                                                                                                    <  1:776080/151‑1 (MQ=255)                       aaTTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACcag                                                                                                                     <  2:757113/150‑1 (MQ=255)                            ggTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGa                                                                                                                >  1:823475/1‑150 (MQ=255)                              tCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATc                                                                                                             >  2:538066/1‑151 (MQ=255)                               cAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATcc                                                                                                            <  2:447252/151‑1 (MQ=255)                                   cATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCAt                                                                                                         <  2:688892/150‑1 (MQ=255)                                             cgccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACggg                                                                                                 >  2:976707/1‑148 (MQ=255)                                               ccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACggg                                                                                                 >  2:190587/1‑146 (MQ=255)                                               ccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATc                                                                                            <  1:367341/151‑1 (MQ=255)                                                 gcTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCAATACAAACAAGAAGAATTCCCACAAGACgg                                                                                                  >  2:261958/1‑143 (MQ=255)                                                    ggCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTc                                                                                       <  1:663909/151‑1 (MQ=255)                                                    ggCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTc                                                                                       <  2:970465/151‑1 (MQ=255)                                                    ggCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTc                                                                                       >  1:610146/1‑151 (MQ=255)                                                      ccAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCa                                                                                      >  2:670687/1‑150 (MQ=255)                                                       cAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAAt                                                                                    <  1:810460/151‑1 (MQ=255)                                                       cAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAAt                                                                                    <  1:606374/151‑1 (MQ=255)                                                          tGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGcc                                                                                 >  1:162780/1‑151 (MQ=255)                                                                    agaaCAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGccc                                                                       >  1:715578/1‑151 (MQ=255)                                                                       aCAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCccacc                                                                    >  2:65890/1‑151 (MQ=255)                                                                        cAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCccacc                                                                    >  2:284956/1‑150 (MQ=255)                                                                            cGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCg                                                               >  2:403070/1‑151 (MQ=255)                                                                             gCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGc                                                              >  1:455755/1‑151 (MQ=255)                                                                               aaTAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCAt                                                            >  2:333919/1‑151 (MQ=255)                                                                                aTAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCAt                                                            >  2:80969/1‑150 (MQ=255)                                                                                  aGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAAt                                                         >  1:571045/1‑151 (MQ=255)                                                                                      gAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTgcg                                                     <  2:113621/151‑1 (MQ=255)                                                                                       aaaGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTgcgc                                                    <  2:78426/151‑1 (MQ=255)                                                                                              cATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCt                                             >  2:897906/1‑151 (MQ=255)                                                                                              cATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCt                                             <  2:680974/151‑1 (MQ=255)                                                                                                ttCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGc                                           >  1:859414/1‑151 (MQ=255)                                                                                                  cGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGt                                         <  1:6183/151‑1 (MQ=255)                                                                                                   gggCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTa                                        <  1:141229/151‑1 (MQ=255)                                                                                                      ccTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGc                                     <  1:738694/151‑1 (MQ=255)                                                                                                           tttAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTccc                                >  2:313325/1‑151 (MQ=255)                                                                                                              aaTGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCcac                             >  1:493510/1‑151 (MQ=255)                                                                                                                 ggTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACaa                           >  1:254574/1‑150 (MQ=255)                                                                                                                       ttttAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAAcc                    >  2:385421/1‑151 (MQ=255)                                                                                                                          tAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGcc                 >  2:486315/1‑151 (MQ=255)                                                                                                                          tAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGcc                 >  2:96652/1‑151 (MQ=255)                                                                                                                               aTTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAg            <  2:346391/151‑1 (MQ=255)                                                                                                                                  tCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGa          >  1:93388/1‑150 (MQ=255)                                                                                                                                  tCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGa          >  1:460559/1‑150 (MQ=255)                                                                                                                                   ccAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACTAGGa          >  2:887253/1‑149 (MQ=255)                                                                                                                                   ccAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGaa         >  2:177984/1‑150 (MQ=255)                                                                                                                                     aGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGaaaaa      <  1:508452/151‑1 (MQ=255)                                                                                                                                         tGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGaaaa       >  1:291259/1‑146 (MQ=255)                                                                                                                                           ccGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATc  >  2:213598/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                               CGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATC  >  NC_000913/1618136‑1618423

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