Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,587,525	Δ1 bp	coding (589/744 nt)	ugpQ ←	glycerophosphodiester phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,587,523	0	T	.	100.0%	218.0 / NA	51	coding (591/744 nt)	ugpQ	glycerophosphodiester phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (22/29);  total (22/29)

TAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTT  >  NC_000913/3587381‑3587653                                                                                                                                               |                                                                                                                                   tAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGC‑‑‑GCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCt                                                                                                                                    >  1:278686/1‑140 (MQ=255)    aGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCt                                                                                                                                    <  2:427384/140‑1 (MQ=255)           aCCAATCACGTCAATCGCATCGGTACAGATGCAATCCACGCCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAg                                                                                                                            <  2:979833/140‑1 (MQ=255)           aCCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAg                                                                                                                            <  2:97594/140‑1 (MQ=255)           aCCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAg                                                                                                                            <  1:340838/140‑1 (MQ=255)           aCCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAg                                                                                                                            <  1:390561/140‑1 (MQ=255)                tCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACt                                                                                                                       >  1:179892/1‑140 (MQ=255)                tCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACt                                                                                                                       >  1:585776/1‑140 (MQ=255)                  aCGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTa                                                                                                                     <  1:28476/140‑1 (MQ=255)                       aaTCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGAtt                                                                                                                <  2:496543/140‑1 (MQ=255)                       aaTCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGAtt                                                                                                                <  1:750573/140‑1 (MQ=255)                       aaTCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGAtt                                                                                                                <  2:115548/140‑1 (MQ=255)                         tCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTga                                                                                                              <  2:179892/140‑1 (MQ=255)                          cGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTgag                                                                                                             <  1:106189/140‑1 (MQ=255)                          cGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTgag                                                                                                             <  1:968550/140‑1 (MQ=255)                          cGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTgag                                                                                                             <  2:106606/140‑1 (MQ=255)                          cGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTgag                                                                                                             <  1:695514/140‑1 (MQ=255)                          cGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTgag                                                                                                             <  2:540282/140‑1 (MQ=255)                          cGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTgag                                                                                                             <  1:648629/140‑1 (MQ=255)                                  tACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATag                                                                                                     <  2:644819/140‑1 (MQ=255)                                       aTGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACg                                                                                                <  2:887937/140‑1 (MQ=255)                                                                      ctgcGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAgt                                                                 <  2:579932/140‑1 (MQ=255)                                                                           cgcTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTcgcg                                                                      <  2:552158/130‑1 (MQ=255)                                                                           cgcTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTcgcg                                                                      >  1:552158/1‑130 (MQ=255)                                                                           cgcTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCa                                                                   >  2:774009/1‑133 (MQ=255)                                                                           cgcTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCa                                                                   <  1:774009/133‑1 (MQ=255)                                                                             cTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTcgcg                                                          <  2:173776/140‑1 (MQ=255)                                                                                      aTTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATc                                                 <  2:278686/140‑1 (MQ=255)                                                                                      aTTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATc                                                 <  2:712743/140‑1 (MQ=255)                                                                                        tGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCAGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCa                                               >  2:185572/1‑140 (MQ=255)                                                                                            ggTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGcaa                                           >  1:279112/1‑140 (MQ=255)                                                                                                taAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACaaa                                       >  1:48547/1‑140 (MQ=255)                                                                                                   aCCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACcg                                    <  2:279112/140‑1 (MQ=255)                                                                                                        aaTCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCgg                               >  1:520413/1‑140 (MQ=255)                                                                                                        aaTCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCgg                               >  1:3978/1‑140 (MQ=255)                                                                                                        aaTCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCgg                               >  2:29035/1‑140 (MQ=255)                                                                                                        aaTCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCgg                               >  1:656833/1‑140 (MQ=255)                                                                                                            cGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGt                           >  1:565701/1‑140 (MQ=255)                                                                                                               aGTCCGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTcc                        >  1:432525/1‑140 (MQ=255)                                                                                                                  ccGGCGTCTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGt                     >  1:270961/1‑140 (MQ=255)                                                                                                                         cTTTCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCt              <  2:520413/140‑1 (MQ=255)                                                                                                                           ttCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCtgt            >  2:433486/1‑140 (MQ=255)                                                                                                                           ttCAACTGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCtgt            >  1:405445/1‑140 (MQ=255)                                                                                                                                 tGCATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCa      <  2:48547/140‑1 (MQ=255)                                                                                                                                   cATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGc       <  2:174597/137‑1 (MQ=255)                                                                                                                                   cATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGc       >  1:174597/1‑137 (MQ=255)                                                                                                                                   cATCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGc    >  2:463011/1‑140 (MQ=255)                                                                                                                                    aTCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGc       <  2:631628/136‑1 (MQ=255)                                                                                                                                    aTCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGc       >  1:631628/1‑136 (MQ=255)                                                                                                                                    aTCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCt   >  1:526176/1‑140 (MQ=255)                                                                                                                                     tCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTtgtg        <  2:769928/134‑1 (MQ=255)                                                                                                                                     tCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTtgtg        >  1:769928/1‑134 (MQ=255)                                                                                                                                     tCACTCGCGC‑TTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCtt  >  2:500904/1‑140 (MQ=255)                                                                                                                                               |                                                                                                                                   TAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTT  >  NC_000913/3587381‑3587653

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