Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,914,439	A→G	R126R (CGT→CGC)	prc ←	tail‑specific protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,914,439	0	A	G	100.0%	201.3 / NA	57	R126R (CGT→CGC)	prc	tail‑specific protease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (35/22);  total (35/22)

ACTTTACTGTCCCACAGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAACGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCTTTTTCGCGAACTGTTCAACATCGCTTGCCAGCAGCACGTTGT  >  NC_000913/1914308‑1914572                                                                                                                                    |                                                                                                                                      aCTTTACTGTCCCACAGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCaaa                                                                                                                               >  1:109924/1‑140 (MQ=255)     tACTGTCCCACAGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACgg                                                                                                                           <  1:35530/140‑1 (MQ=255)          tCCCACAGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCtt                                                                                                                      >  2:887027/1‑140 (MQ=255)          tCCCACAGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCtt                                                                                                                      >  2:772571/1‑140 (MQ=255)           cccACAGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCttt                                                                                                                     >  2:1062200/1‑140 (MQ=255)               caGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTgcg                                                                                                                 <  1:1124343/140‑1 (MQ=255)               caGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTgcg                                                                                                                 >  2:231349/1‑140 (MQ=255)               caGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTgcg                                                                                                                 <  1:843479/140‑1 (MQ=255)                 gcgcGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGcc                                                                                                               >  1:1122761/1‑140 (MQ=255)                  cgcgTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCa                                                                                                              <  1:772571/140‑1 (MQ=255)                    cgTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGa                                                                                                            >  2:1025506/1‑140 (MQ=255)                     gTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGAt                                                                                                           <  1:1025506/140‑1 (MQ=255)                      ttCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTgcg                                                                                                                 >  2:495943/1‑133 (MQ=255)                      ttCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTgcg                                                                                                                 <  1:495943/133‑1 (MQ=255)                          aCTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTag                                                                                                      <  2:900039/140‑1 (MQ=255)                            tCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTagag                                                                                                    >  1:64151/1‑140 (MQ=255)                              aGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGAt                                                                                                  >  1:246262/1‑140 (MQ=255)                                     tttttCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGaa                                                                                           >  2:403946/1‑140 (MQ=255)                                        ttCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCACTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAAc                                                                                        >  1:318172/1‑140 (MQ=255)                                                gggCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGtt                                                                                <  2:246262/140‑1 (MQ=255)                                                gggCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGtt                                                                                >  1:657940/1‑140 (MQ=255)                                                gggCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGtt                                                                                >  2:533771/1‑140 (MQ=255)                                                gggCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGtt                                                                                <  1:887027/140‑1 (MQ=255)                                                  gcgcTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCAACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTg                                                                              >  1:708107/1‑140 (MQ=255)                                                           tgcGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGc                                                                     >  2:639569/1‑140 (MQ=255)                                                             cGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCGGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAg                                                                   >  2:631623/1‑140 (MQ=255)                                                               gTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGtt                                                                 >  2:587485/1‑140 (MQ=255)                                                                 cAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCa                                                               <  2:639907/140‑1 (MQ=255)                                                                 cAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCa                                                               <  1:587485/140‑1 (MQ=255)                                                                  aaGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGc                                                                             >  1:433349/1‑125 (MQ=255)                                                                  aaGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGc                                                                             <  2:433349/125‑1 (MQ=255)                                                                         taAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTaa                                                       >  2:731639/1‑140 (MQ=255)                                                                            gtgtCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTc                                                    >  2:1018748/1‑140 (MQ=255)                                                                              gtCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCgg                                                  >  2:116640/1‑140 (MQ=255)                                                                               tCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGt                                                 <  2:64151/140‑1 (MQ=255)                                                                               tCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGt                                                 >  2:817052/1‑140 (MQ=255)                                                                               tCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGt                                                 >  2:87082/1‑140 (MQ=255)                                                                               tCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGt                                                 >  2:263162/1‑140 (MQ=255)                                                                                       ggTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCtttt                                         >  2:953820/1‑140 (MQ=255)                                                                                         tGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCTTTTTc                                       <  1:87082/140‑1 (MQ=255)                                                                                            aaTCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCTTTTTcgcg                                    <  1:403946/140‑1 (MQ=255)                                                                                            aaTCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCTTTTTcgcg                                    <  1:231349/140‑1 (MQ=255)                                                                                              tCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCTTTTTCGCGaa                                  <  1:731639/140‑1 (MQ=255)                                                                                                   cGGCTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCTTTTTCGCGAACTGtt                             >  1:552212/1‑140 (MQ=255)                                                                                                      cTTTTCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCTTTTTCGCGAACTGTTCaa                          <  1:1018748/140‑1 (MQ=255)                                                                                                         ttCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCTTTTTCGCGAACTGTTCAACAt                       <  2:109924/140‑1 (MQ=255)                                                                                                          tCCAGTACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCTTTTTCGCGAACTGTTCAACATc                      >  1:608355/1‑140 (MQ=255)                                                                                                               tACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCTTTTTCGCGAACTGTTCAACATCGCTTg                 >  2:339804/1‑140 (MQ=255)                                                                                                               tACCGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCTTTTTCGCGAACTGTTCAACATCGCTTg                 >  1:422930/1‑140 (MQ=255)                                                                                                                 ccGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCTTTTTCGCGAACTGTTCAACATCGCTTGcc               >  2:930403/1‑140 (MQ=255)                                                                                                                 ccGACAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTCTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCTTTTTCGCGAACTGTTCAACATCGCTTGcc               <  1:1062200/140‑1 (MQ=255)                                                                                                                    aCAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCTTTTTCGCGAACTGTTCAACATCGCTTGCcagc            >  2:498759/1‑140 (MQ=255)                                                                                                                    aCAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCTTTTTCGCGAACTGTTCAACATCGCTTGCcagc            <  2:628006/140‑1 (MQ=255)                                                                                                                    aCAAAGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCTTTTTCGCGAACTGTTCAACATCGCTTGCcagc            >  2:758559/1‑140 (MQ=255)                                                                                                                        aGCGTACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCTTTTTCGCGAACTGTTCAACATCGCTTGCCAgcagca        >  2:113399/1‑140 (MQ=255)                                                                                                                            tACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCTTTTTCGCGAACTGTTCAACATCGCTTGCCAGCAGCACGtt    >  1:489131/1‑140 (MQ=255)                                                                                                                            tACTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCTTTTTCGCGAACTGTTCAACATCGCTTGCCAGCAGCACGtt    <  1:953820/140‑1 (MQ=255)                                                                                                                              cTGGTAGCGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCTTTTTCGCGAACTGTTCAACATCGCTTGCCAGCAGCACGTtgt  >  1:1000364/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                                      ACTTTACTGTCCCACAGCGCGTTCAACTCAGCCTCGTTTTTCGGCCAGGGCGCTTTGCTGCGGTCAAGGTTATAAGTGTCGTTGCCGGTGAAATCCATCGGCTTTTCCAGTACCGACAAAGCGTACTGGTAACGCTCAAAACGGCGCTTTTGCGCCAGATTGTAGAGATCGTAGAAAACGTCGAGTTTGCCTGAACGCAGTTCATCGCCTAACTCGGTTTTCTTTTTCGCGAACTGTTCAACATCGCTTGCCAGCAGCACGTTGT  >  NC_000913/1914308‑1914572

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