Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,996,372	C→T	*302* (TAG→TAA)	insD1 ←	IS2 transposase TnpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,996,372	0	C	T	100.0%	122.9 / NA	39	*302* (TAG→TAA)	insD1	IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/20);  total (19/20)

TAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGC  >  NC_000913/2996241‑2996491                                                                                                                                    |                                                                                                                        taaataaTCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATAttt                                                                                                                  >  2:547957/1‑139 (MQ=255)  aaataaTCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTta                                                                                                                       <  2:287616/133‑1 (MQ=255)  aaataaTCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTta                                                                                                                       >  1:287616/1‑133 (MQ=255)  aaataaTCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTc                                                                                                                 >  2:416842/1‑139 (MQ=255)       atCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGAc                                                                                                            >  1:304108/1‑139 (MQ=255)         cccTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACAt                                                                                                          <  2:146744/139‑1 (MQ=255)              aGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGtt                                                                                                     <  1:416842/139‑1 (MQ=255)               gTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTtat                                                                                                                      >  2:329439/1‑121 (MQ=255)               gTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTtat                                                                                                                      <  1:329439/121‑1 (MQ=255)                   aTCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAATTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACa                                                                                                           <  2:787358/128‑1 (MQ=255)                   aTCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACa                                                                                                           >  1:787358/1‑128 (MQ=255)                     cAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACtt                                                                                              >  1:307784/1‑139 (MQ=255)                        gAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAc                                                                                           >  2:182178/1‑139 (MQ=255)                         aaCCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                          <  1:561769/139‑1 (MQ=255)                         aaCCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                          <  2:307784/139‑1 (MQ=255)                             ttATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCAtt                                                                                      >  2:58451/1‑139 (MQ=255)                                      aaCCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                             <  1:547957/139‑1 (MQ=255)                                      aaCCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                             <  2:304108/139‑1 (MQ=255)                                           gcgTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                          >  1:346374/1‑121 (MQ=255)                                           gcgTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                          <  2:346374/121‑1 (MQ=255)                                           gcgTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGcc                                                                        <  2:946718/139‑1 (MQ=255)                                           gcgTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGcc                                                                        <  1:182178/139‑1 (MQ=255)                                                cTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGa                                                                   <  2:1004200/139‑1 (MQ=255)                                                           tctcGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTgg                                                        >  2:344659/1‑139 (MQ=255)                                                                         aCAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGc                                                                      <  1:621477/111‑1 (MQ=255)                                                                         aCAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGc                                                                      >  2:621477/1‑111 (MQ=255)                                                                           aaCCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgc                                        >  2:361495/1‑139 (MQ=255)                                                                            aCCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgcg                                       <  1:58451/139‑1 (MQ=255)                                                                                 tGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGa                                                                   >  2:332092/1‑106 (MQ=255)                                                                                 tGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGa                                                                   <  1:332092/106‑1 (MQ=255)                                                                                        aCAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGc                                                                            <  2:282980/90‑1 (MQ=255)                                                                                        aCAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGc                                                                            >  1:282980/1‑90 (MQ=255)                                                                                             tCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                          <  2:413122/71‑1 (MQ=255)                                                                                             tCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                          >  1:413122/1‑71 (MQ=255)                                                                                             tCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCAt                      >  1:477870/1‑139 (MQ=255)                                                                                                  ggAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGtt                 >  2:817038/1‑139 (MQ=255)                                                                                                      gTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATaa             >  1:1097774/1‑139 (MQ=255)                                                                                                              cAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGaa     <  2:1097774/139‑1 (MQ=37)                                                                                                                 ttGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGc  <  2:574644/139‑1 (MQ=21)                                                                                                                                    |                                                                                                                        TAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGC  >  NC_000913/2996241‑2996491

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