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%	247.8 / NA	51	coding (628/666 nt)	marC	UPF0056 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (24/27);  total (24/27)

CGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGC  >  NC_000913/1618136‑1618426                                                                                                                                                  |                                                                                                                                                  cgGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATaaa                                                                                                                                             <  1:568364/151‑1 (MQ=255)   gTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACt                                                                                                                                           >  1:562478/1‑151 (MQ=255)        tGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGcac                                                                                                                                       >  2:781287/1‑150 (MQ=255)         gCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGcac                                                                                                                                       <  2:298458/149‑1 (MQ=255)              tGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCcat                                                                                                                                 >  2:575789/1‑150 (MQ=255)               gCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCcata                                                                                                                                <  1:483869/150‑1 (MQ=255)                ccGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCcatac                                                                                                                               >  1:286364/1‑150 (MQ=255)                ccGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCcatac                                                                                                                               >  2:399015/1‑150 (MQ=255)                 cgcgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCtatacat                                                                                                                             >  1:905002/1‑151 (MQ=255)                 cgcgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCcatacat                                                                                                                             >  2:40152/1‑151 (MQ=255)                 cgcgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCcatacat                                                                                                                             >  2:204710/1‑151 (MQ=255)                 cgcgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCcatacat                                                                                                                             <  2:372919/151‑1 (MQ=255)                  gcgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑ACTAAACTGCACACCcatacata                                                                                                                            <  1:754043/151‑1 (MQ=255)                  gcgTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCcatacata                                                                                                                            <  1:521586/151‑1 (MQ=255)                        aTTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACcagca                                                                                                                      <  2:713621/151‑1 (MQ=255)                            ggTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGaa                                                                                                                  <  1:272460/151‑1 (MQ=255)                            ggTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGa                                                                                                                   <  1:573086/150‑1 (MQ=255)                                aaTCATTATATTTCGCCGCTGGCCAGTGATACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATccc                                                                                                              <  2:67146/151‑1 (MQ=255)                                       atatTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGAc                                                                                                       <  2:539709/151‑1 (MQ=255)                                            tCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACggg                                                                                                    >  2:365269/1‑149 (MQ=255)                                              gccgcTGGCCAGTGACACTCGAAGAACAGACGCAATTGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGAc                                                                                                       >  1:56644/1‑144 (MQ=255)                                               ccgcTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACggg                                                                                                    >  2:59813/1‑146 (MQ=255)                                                cgcTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCg                                                                                              <  2:436121/151‑1 (MQ=255)                                                   tGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCtt                                                                                           >  2:282028/1‑151 (MQ=255)                                                         gTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGc                                                                                     <  2:927530/151‑1 (MQ=255)                                                                         aGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCccacc                                                                       >  2:233602/1‑149 (MQ=255)                                                                           aCGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGcc                                                                   <  1:549065/151‑1 (MQ=255)                                                                            cGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCg                                                                  <  2:453484/151‑1 (MQ=255)                                                                                aTAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCAt                                                               >  2:419549/1‑150 (MQ=255)                                                                                aTAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCAt                                                               >  2:466400/1‑150 (MQ=255)                                                                                    ttGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTg                                                          >  2:892553/1‑151 (MQ=255)                                                                                       aaaGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTgcgc                                                       <  1:204701/151‑1 (MQ=255)                                                                                          ggCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCgg                                                    <  2:556150/151‑1 (MQ=255)                                                                                            cccATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCgg                                                    >  2:89865/1‑149 (MQ=255)                                                                                              cATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCt                                                <  2:745239/151‑1 (MQ=255)                                                                                                         tttttAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCtt                                      <  1:513385/150‑1 (MQ=255)                                                                                                               aTGGTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCcac                                >  2:340171/1‑150 (MQ=255)                                                                                                                 ggTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACaa                              <  2:44678/150‑1 (MQ=255)                                                                                                                  gTACGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACaaaa                            <  2:286391/151‑1 (MQ=255)                                                                                                                     cGTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATa                          >  2:451127/1‑150 (MQ=255)                                                                                                                      gTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAAc                        >  2:866684/1‑151 (MQ=255)                                                                                                                      gTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAAc                        >  1:152419/1‑151 (MQ=255)                                                                                                                      gTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAAc                        >  2:65389/1‑151 (MQ=255)                                                                                                                      gTTTTAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAAc                        <  2:638951/151‑1 (MQ=255)                                                                                                                         ttAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCg                      <  1:398975/150‑1 (MQ=255)                                                                                                                         ttAATGATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCg                      <  2:271478/150‑1 (MQ=255)                                                                                                                              gATTTCCAGGATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGccacca                <  2:867298/151‑1 (MQ=255)                                                                                                                                      ggATGCCGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGaaaa          >  1:856197/1‑149 (MQ=255)                                                                                                                                           ccGTTAA‑AATAAACTGCACACCTATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAg   <  2:905101/151‑1 (MQ=255)                                                                                                                                           ccGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAg   <  1:365233/151‑1 (MQ=255)                                                                                                                                            cGTTAA‑AATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAgc  >  1:195869/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                  CGGTGCCTGCTTTTGCCGCGTTAATTTGGTCAATCATTATATTTCGCCGCTGGCCAGTGACACTCGAAGAACAGACGCAATAGTTGAAAGGCCCATTCGGGCCTTTTTTAATGGTACGTTTTAATGATTTCCAGGATGCCGTTAATAATAAACTGCACACCCATACATACCAGCAGGAATCCCATCAGACGGGAGATCGCTTCAATGCCACCCTTGCCCACCAGCCGCATAATTGCGCCGGAGCTGCGTAGGCTTCCCCACAAAATAACCGCCACCAGGAAAAAGATCAGC  >  NC_000913/1618136‑1618426

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: ‑