Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	973,216	C→T	L62L (CTG→CTA)	elyC ←	envelope biogenesis factor; DUF218 superfamily protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	973,216	0	C	T	100.0%	202.0 / NA	52	L62L (CTG→CTA)	elyC	envelope biogenesis factor; DUF218 superfamily protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/25);  total (27/25)

CAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTG  >  NC_000913/973075‑973358                                                                                                                                              |                                                                                                                                               caGGCGTGGCAGACTGTTATTGATTAAATTCGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCg                                                                                                                                       >  1:370752/1‑151 (MQ=255)     cGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCt                                                                                                                                    >  1:87902/1‑150 (MQ=255)     cGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCt                                                                                                                                    >  1:703109/1‑150 (MQ=255)         gCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCgg                                                                                                                               >  1:375336/1‑151 (MQ=255)            gACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTg                                                                                                                            <  2:469874/151‑1 (MQ=255)             aCTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGc                                                                                                                           >  1:33465/1‑151 (MQ=255)                gTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCaa                                                                                                                         <  2:245637/150‑1 (MQ=255)                  tATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGc                                                                                                                                            <  2:260649/129‑1 (MQ=255)                  tATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGc                                                                                                                                            >  1:260616/1‑129 (MQ=255)                   aTTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACt                                                                                                                      <  2:639113/150‑1 (MQ=255)                   aTTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTc                                                                                                                     <  1:510445/151‑1 (MQ=255)                      gATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCaa                                                                                                                   >  2:603906/1‑150 (MQ=255)                               agagCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAgc                                                                                                         <  1:498161/151‑1 (MQ=255)                                  gCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCgcc                                                                                                      >  1:237895/1‑151 (MQ=255)                                          cccACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCAccc                                                                                              >  1:512181/1‑151 (MQ=255)                                             aCTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAAt                                                                                           <  1:441576/151‑1 (MQ=255)                                             aCTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAAt                                                                                           <  2:313728/151‑1 (MQ=255)                                                gCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCt                                                                                        <  2:126898/151‑1 (MQ=255)                                                 cGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCAccc                                                                                              <  1:3229/144‑1 (MQ=255)                                                 cGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCAccc                                                                                              >  2:3229/1‑144 (MQ=255)                                                       ccAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATaaa                                                                                  >  2:415208/1‑150 (MQ=255)                                                           gTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCt                                                                             >  2:193029/1‑151 (MQ=255)                                                            tataACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCtt                                                                            >  1:738076/1‑151 (MQ=255)                                                             ataACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTc                                                                           >  2:468975/1‑151 (MQ=255)                                                             ataACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTc                                                                           <  1:6118/151‑1 (MQ=255)                                                                  caccTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCgg                                                                       >  2:499413/1‑150 (MQ=255)                                                                   accTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGtt                                                                     <  2:375383/151‑1 (MQ=255)                                                                         gAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTg                                                                >  1:350725/1‑150 (MQ=255)                                                                          aGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGaa                                                              <  2:87933/151‑1 (MQ=255)                                                                              ccacGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACg                                                          <  1:24853/151‑1 (MQ=255)                                                                               cacGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACgg                                                         >  2:368997/1‑151 (MQ=255)                                                                                   atgtaGTCCACTTTTTGCGAATTATTCCATGTCGGATTGGTTCTTTCGATGGGACGCCTTACGCGGTCGTCCGCGGGTTGCAAACTCAATAACACGAGCGCCAGCCACCCAATGCTAATAAACATCTTCCCGGATTTGTGCAAACGGCTa                                                      >  1:514590/4‑150 (MQ=255)                                                                                    tataGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCTGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGGGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTaa                                                     >  1:246965/1‑150 (MQ=255)                                                                                       aGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAAcc                                                  >  2:31481/1‑150 (MQ=255)                                                                                       aGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACca                                                 <  2:629440/151‑1 (MQ=255)                                                                                            aCTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCAcagca                                            <  2:350760/151‑1 (MQ=255)                                                                                              tttttGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACagcagc                                          >  2:650078/1‑151 (MQ=255)                                                                                               ttttGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACagcagc                                          <  2:624855/150‑1 (MQ=255)                                                                                                   gCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCgccag                                     <  1:268068/151‑1 (MQ=255)                                                                                                    cGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCgccagg                                    <  2:604713/151‑1 (MQ=255)                                                                                                         tattCCATGTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAgc                               >  2:86177/1‑151 (MQ=255)                                                                                                                 gTCGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCaataa                       <  2:32769/151‑1 (MQ=255)                                                                                                                   cGGATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCaataat                      <  2:709107/150‑1 (MQ=255)                                                                                                                    ggATAGGTGCTTTCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATc                     >  1:254005/1‑150 (MQ=255)                                                                                                                               ttCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCa           <  2:23084/149‑1 (MQ=255)                                                                                                                               ttCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCAGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGc         >  1:352773/1‑151 (MQ=255)                                                                                                                                tCGATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATTAGCg        >  2:364386/1‑151 (MQ=255)                                                                                                                                  gATGGGACGCAGTAGGTGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCgg       >  2:39321/1‑150 (MQ=255)                                                                                                                                  gATGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCgg       <  2:594073/150‑1 (MQ=255)                                                                                                                                    tGGGACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGa      >  1:32308/1‑149 (MQ=255)                                                                                                                                      ggACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTg  <  1:528153/151‑1 (MQ=255)                                                                                                                                      ggACGCAGTAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTg  <  2:639936/151‑1 (MQ=255)                                                                                                                                              |                                                                                                                                               CAGGCGTGGCAGACTGTTATTGATTAAATTAGAGCTCGGTGCCCACTGCGGGTTCCAGGTATAACCACCTCCGAGCACCACGATATAGTCCACTTTTTGCGAATTATTCCATGTCGGATAGGTGCTTTCGATGGGACGCAGCAGGCGGTCGGCTACCGGTTGCAAACTCAATAACAAGAGCGCCAGCCACCCAATGCTAATAAAAATCTTCCCGGTTTTTTGAAAACGGCTAAACCACAGCAGCGCCAGGCCAGCGCCAATAATCAACAACATCAGCGGAAGTG  >  NC_000913/973075‑973358

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