Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	922,366	C→T	*75* (TAG→TAA)	cspD ←	inhibitor of DNA replication, cold shock protein homolog

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	922,366	0	C	T	100.0%	239.6 / NA	73	*75* (TAG→TAA)	cspD	inhibitor of DNA replication, cold shock protein homolog
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (38/35);  total (38/35)

TACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAATGGTGGAAT  >  NC_000913/922238‑922499                                                                                                                                 |                                                                                                                                      tACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTg                                                                                                                             >  1:33805/1‑139 (MQ=255)                cGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGa                                                                                                              <  2:14240/139‑1 (MQ=255)                  gCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACg                                                                                                            <  1:1006197/139‑1 (MQ=255)                                       gCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACggg                                                                                                          >  1:119053/1‑120 (MQ=255)                                       gCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACggg                                                                                                          <  2:119053/120‑1 (MQ=255)                                            ggCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCg                                                                                                               <  2:1376161/110‑1 (MQ=255)                                            ggCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCg                                                                                                               >  1:1376161/1‑110 (MQ=255)                                             gCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGcc                                                                                                                           <  2:864900/97‑1 (MQ=255)                                             gCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGcc                                                                                                                           >  1:864900/1‑97 (MQ=255)                                                     gTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGccc                                                                         <  2:33805/139‑1 (MQ=255)                                                             tGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGa                                                                 >  1:1393656/1‑139 (MQ=255)                                                              gagaTAAAAATGCCAGCCAATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGAc                                                                >  1:1508505/1‑139 (MQ=255)                                                                gaTAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACg                                                                                                            >  1:532101/1‑93 (MQ=255)                                                                gaTAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACg                                                                                                            <  2:532101/93‑1 (MQ=255)                                                                        tGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGa                                                      <  1:969672/139‑1 (MQ=255)                                                                          ccAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACgg                                                                                                           >  2:850829/1‑84 (MQ=255)                                                                          ccAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACgg                                                                                                           <  1:850829/84‑1 (MQ=255)                                                                                gATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTg                                              >  1:22343/1‑139 (MQ=255)                                                                                 aTCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGAt                                                                                      <  2:926528/98‑1 (MQ=255)                                                                                 aTCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGAt                                                                                      >  1:926528/1‑98 (MQ=255)                                                                                        tGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGc                                                                                                         >  2:1038404/1‑72 (MQ=255)                                                                                        tGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGc                                                                                                         <  1:1038404/72‑1 (MQ=255)                                                                                        tGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTg                                                                                             >  1:1063859/1‑84 (MQ=255)                                                                                        tGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTg                                                                                             <  2:1063859/84‑1 (MQ=255)                                                                                            aTTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAg                                  <  1:1817252/139‑1 (MQ=255)                                                                                                tGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggtgg                                                                  >  2:1851762/1‑103 (MQ=255)                                                                                                tGCCTTTAGGATGTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggtgg                                                                  <  1:1851762/103‑1 (MQ=255)                                                                                                    tttAGGATGTAC‑CAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTa                         >  2:1493777/1‑139 (MQ=255)                                                                                                            gTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAg                                  <  2:202987/123‑1 (MQ=255)                                                                                                            gTACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAg                                  >  1:202987/1‑123 (MQ=255)                                                                                                             tACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACg                                                                                                            >  1:490352/1‑48 (MQ=255)                                                                                                             tACACAATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACg                                                                                                            <  2:490352/48‑1 (MQ=255)                                                                                                               cacaATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGc                                                                           <  1:456108/79‑1 (MQ=255)                                                                                                               cacaATGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGc                                                                           >  2:456108/1‑79 (MQ=255)                                                                                                                  aaTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACg                                                                                                            >  1:514669/1‑43 (MQ=255)                                                                                                                  aaTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACg                                                                                                            <  2:514669/43‑1 (MQ=255)                                                                                                                   aTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACggg                                                                                                          >  1:501361/1‑44 (MQ=255)                                                                                                                   aTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACggg                                                                                                          <  2:501361/44‑1 (MQ=255)                                                                                                                   aTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTg                                                                             <  1:1876930/73‑1 (MQ=255)                                                                                                                   aTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTg                                                                             >  2:1876930/1‑73 (MQ=255)                                                                                                                   aTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa           <  1:1493777/139‑1 (MQ=255)                                                                                                                   aTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa           >  2:497088/1‑139 (MQ=255)                                                                                                                   aTGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGaa           <  2:1508505/139‑1 (MQ=255)                                                                                                                   aTGAGACAGAAGAGTTATCCGACTGCCGCTTCTACTTGGACGGGCACAATAACACTTGCCTGATTGCCTTTTGGCCGCTGATGGACATTAAACTGACCGTATTGTCCAGCTTTTAGCGTCCTGTAACCATCCATCTGaa           >  1:43381/1‑139 (MQ=255)                                                                                                                    tGAGCCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATACCACTGGCGTGATTGCCTTTTGGCCCCtggtg                                                                   <  2:806225/82‑1 (MQ=255)                                                                                                                    tGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACggg                                                                                                          <  1:865830/43‑1 (MQ=255)                                                                                                                    tGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACggg                                                                                                          >  2:865830/1‑43 (MQ=255)                                                                                                                    tGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAac                                                                                                 >  2:1209476/1‑52 (MQ=255)                                                                                                                    tGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAac                                                                                                 <  1:1209476/52‑1 (MQ=255)                                                                                                                    tGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtggtg                                                                   >  1:806225/1‑82 (MQ=255)                                                                                                                    tGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGa                                                 <  1:1668583/100‑1 (MQ=255)                                                                                                                    tGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGa                                                 >  2:1668583/1‑100 (MQ=255)                                                                                                                    tGAGACAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAAt          >  2:1804561/1‑139 (MQ=255)                                                                                                                     gagaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAAtg         <  2:537140/139‑1 (MQ=255)                                                                                                                       gaCAGAAGAGTTATGCGTCTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTg                                                                             >  2:1175816/1‑69 (MQ=255)                                                                                                                       gaCAGAAGAGTTATGCGCCTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAAtggt       >  2:1177419/1‑139 (MQ=255)                                                                                                                       gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTg                                                                             >  2:1608391/1‑69 (MQ=255)                                                                                                                       gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTg                                                                             <  1:1175816/69‑1 (MQ=255)                                                                                                                       gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTg                                                                             <  1:1608391/69‑1 (MQ=255)                                                                                                                       gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGc                                                                           <  2:988842/71‑1 (MQ=255)                                                                                                                       gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGc                                                                           >  1:988842/1‑71 (MQ=255)                                                                                                                       gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGa                                                 >  2:1195375/1‑97 (MQ=255)                                                                                                                       gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGa                                                 <  1:1195375/97‑1 (MQ=255)                                                                                                                       gaCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGG‑‑CAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAAtggtgg     >  2:1642703/1‑139 (MQ=255)                                                                                                                        aCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGa                                                 >  1:904426/1‑96 (MQ=255)                                                                                                                        aCAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGa                                                 <  2:904426/96‑1 (MQ=255)                                                                                                                         cAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTg                                                                             <  2:718717/67‑1 (MQ=255)                                                                                                                         cAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTg                                                                             >  1:718717/1‑67 (MQ=255)                                                                                                                         cAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGcccc                                                                        >  2:1692431/1‑72 (MQ=255)                                                                                                                         cAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGcccc                                                                        <  1:1692431/72‑1 (MQ=255)                                                                                                                         cAGAAGAGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAAtggtgg     >  1:753000/1‑139 (MQ=255)                                                                                                                          agaagaGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAATGGTGGa    <  1:983937/139‑1 (MQ=255)                                                                                                                            aagaGTTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAATGGTGGaat  >  2:1334572/1‑139 (MQ=255)                                                                                                                                 |                                                                                                                                      TACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTTTTAGCGTTCTGTAACCATCCATCTGAATGGTGGAAT  >  NC_000913/922238‑922499

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