Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,094,275	C→T	*289* (TAG→TAA)	insF1 ←	IS3 transposase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,094,275	0	C	T	100.0%	130.8 / NA	44	*289* (TAG→TAA)	insF1	IS3 transposase B
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/25);  total (19/25)

TATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCC  >  NC_000913/1094139‑1094381                                                                                                                                         |                                                                                                           tATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACGATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTa                                                                                                          >  2:1321553/1‑139 (MQ=255)   tCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAg                                                                                                        >  1:1028002/1‑139 (MQ=255)     gACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCg                                                                                                      >  2:789712/1‑139 (MQ=255)          gCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGtt                                                                                                 >  2:104566/1‑139 (MQ=255)          gCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGtt                                                                                                 >  2:415297/1‑139 (MQ=255)            tGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCt                                                                                               <  1:162536/139‑1 (MQ=255)                      cTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCaaa                                                                                     <  1:1291352/139‑1 (MQ=255)                        ttACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAAtt                                                                                   >  1:473150/1‑139 (MQ=255)                                 gAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGAc                                                                          >  1:671111/1‑139 (MQ=255)                                    agagACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGa                                                                       <  1:626837/139‑1 (MQ=255)                                     gagaCCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAg                                                                      <  1:415297/139‑1 (MQ=255)                                                  ataaGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCaa                                                         >  1:78425/1‑139 (MQ=255)                                                       gtggtTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACtgtg                                                    >  1:73983/1‑139 (MQ=255)                                                           ttGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAg                                                                      >  1:1154273/1‑117 (MQ=255)                                                           ttGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAg                                                                      <  2:1154273/117‑1 (MQ=255)                                                             gtgAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAg                                                                                                        >  1:338758/1‑81 (MQ=255)                                                             gtgAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAg                                                                                                        <  2:338758/81‑1 (MQ=255)                                                                 gggATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCAc                                          <  1:104566/139‑1 (MQ=255)                                                                 gggATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCAc                                          <  1:1045455/139‑1 (MQ=255)                                                                 gggATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCAc                                          <  2:73983/139‑1 (MQ=255)                                                                  ggATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCAcc                                         <  1:1243051/139‑1 (MQ=255)                                                                    aTGTGGTGAATTTTTTTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGa                                       <  2:671111/139‑1 (MQ=255)                                                                    aTGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGa                                       <  2:473150/139‑1 (MQ=255)                                                                    aTGTGATGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGa                                       <  1:258059/139‑1 (MQ=255)                                                                     tgtgGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGAt                                      >  1:1036462/1‑139 (MQ=255)                                                                       tggtgAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTg                                    >  2:56893/1‑139 (MQ=255)                                                                        ggtgAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGa                                                                       <  1:863953/103‑1 (MQ=255)                                                                        ggtgAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGa                                                                       >  2:863953/1‑103 (MQ=255)                                                                         gtgAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTc                                                                                                <  2:1266339/77‑1 (MQ=255)                                                                         gtgAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTc                                                                                                >  1:1266339/1‑77 (MQ=255)                                                                           gAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAAt                                <  2:1028002/139‑1 (MQ=255)                                                                               ttttGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATcaca                            <  1:675007/139‑1 (MQ=255)                                                                                      aGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGtt                                                                                >  1:597554/1‑80 (MQ=255)                                                                                      aGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGtt                                                                                <  2:597554/80‑1 (MQ=255)                                                                                      aGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGata                     <  2:1171437/139‑1 (MQ=255)                                                                                            cATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTa               >  2:309822/1‑139 (MQ=255)                                                                                                  gCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCaaa                                                                                     <  1:782358/63‑1 (MQ=255)                                                                                                  gCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCaaa                                                                                     >  2:782358/1‑63 (MQ=255)                                                                                                     gCATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGtt      <  2:78425/139‑1 (MQ=18)                                                                                                      cATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATaa                  <  2:1395113/126‑1 (MQ=18)                                                                                                      cATTTTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATaa                  >  1:1395113/1‑126 (MQ=18)                                                                                                        ttttGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGcc   <  2:1357054/139‑1 (MQ=14)                                                                                                         tttGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGccc  <  1:309822/139‑1 (MQ=11)                                                                                                         tttGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGccc  <  1:14711/139‑1 (MQ=11)                                                                                                                                         |                                                                                                           TATCGACTGGCTGATAACGCACTTTACGAGGCGAAAGAGACCGGGCGTAATAAGGTGGTTGTGAGGGATGTGGTGAATTTTTGTGAGTCACCATAAAGCGGCATTTTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCC  >  NC_000913/1094139‑1094381

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