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%	201.1 / NA	63	*302* (TAG→TAA)	insD1	IS2 transposase TnpB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/31);  total (32/31)

TAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCT  >  NC_000913/2996241‑2996493                                                                                                                                    |                                                                                                                          taaataaTCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATAttt                                                                                                                    >  2:1553818/1‑139 (MQ=255) taaataaTCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATAttt                                                                                                                    >  1:1033367/1‑139 (MQ=255)    ataatCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGTCAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCa                                                                                                                 <  1:792846/139‑1 (MQ=255)              aGTACAGCTGAAAACTTCTAAATAATCCCGCGTTCTTCAAACCTTTCTCGATAAAAGTTACAACCACTGAATAAACAAATCAATGGAAGTACATTCCAAGTGTTTATTGAAATTGCCCTTATATCTCCAGACATCGGtt                                                                                                       <  1:1264285/139‑1 (MQ=255)                  cATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCa                                                                                                   >  1:394367/1‑139 (MQ=255)                     cAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACa                                                                                                             >  2:927982/1‑126 (MQ=255)                     cAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACtt                                                                                                >  1:1126710/1‑139 (MQ=255)                     cAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACtt                                                                                                >  1:46617/1‑139 (MQ=255)                     cAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTACCCTTATATTTCCAGACa                                                                                                             <  1:927982/126‑1 (MQ=255)                         aaCCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCTATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                            <  2:1126710/139‑1 (MQ=255)                         aaCCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                            <  1:1553818/139‑1 (MQ=255)                          aCCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAccc                                                                                           >  2:1654175/1‑139 (MQ=255)                               ataAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTtat                                                                                                                        >  2:1380479/1‑105 (MQ=255)                               ataAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTtat                                                                                                                        <  1:1380479/105‑1 (MQ=255)                                      aaCCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTc                                                                                                                   >  2:898914/1‑103 (MQ=255)                                      aaCCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTc                                                                                                                   <  1:898914/103‑1 (MQ=255)                                      aaCCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                               <  2:394367/139‑1 (MQ=255)                                                           tctcGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCg                                                             <  1:1105030/136‑1 (MQ=255)                                                           tctcGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCg                                                             >  2:1105030/1‑136 (MQ=255)                                                                  aaaaGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                            >  2:738521/1‑98 (MQ=255)                                                                  aaaaGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                            <  1:738521/98‑1 (MQ=255)                                                                       tGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGtt                                                                                                       <  1:349146/82‑1 (MQ=255)                                                                       tGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGtt                                                                                                       >  2:349146/1‑82 (MQ=255)                                                                         aCAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCa                                            <  1:1654175/139‑1 (MQ=255)                                                                           aaCCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                            <  2:1126881/89‑1 (MQ=255)                                                                           aaCCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                            >  1:1126881/1‑89 (MQ=255)                                                                             ccACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                               >  2:604646/1‑100 (MQ=255)                                                                             ccACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                               <  1:604646/100‑1 (MQ=255)                                                                                        aCAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                            <  1:1449364/76‑1 (MQ=255)                                                                                        aCAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAAcc                                                                                            >  2:1449364/1‑76 (MQ=255)                                                                                        aCAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAg                                           <  1:870541/125‑1 (MQ=255)                                                                                        aCAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAg                                           >  2:870541/1‑125 (MQ=255)                                                                                        aCAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGAt                             >  1:513924/1‑139 (MQ=255)                                                                                         cAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                                 >  1:753125/1‑70 (MQ=255)                                                                                         cAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                                 <  2:753125/70‑1 (MQ=255)                                                                                           gATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgcg                                         <  1:833186/124‑1 (MQ=255)                                                                                           gATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgcg                                         >  2:833186/1‑124 (MQ=255)                                                                                           gATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGcc                          <  2:1033367/139‑1 (MQ=255)                                                                                            aTCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCa                         >  1:310898/1‑139 (MQ=255)                                                                                            aTCAATGGAAGTAAATTCCAATTGTTTAGTGGATTTGCCCTTATATTTCCAGACA‑CTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCAt                        <  2:1280576/139‑1 (MQ=255)                                                                                             tCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGc                                                                           <  2:738292/88‑1 (MQ=255)                                                                                             tCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGc                                                                           >  1:738292/1‑88 (MQ=255)                                                                                             tCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGata                                                                   <  1:500904/96‑1 (MQ=255)                                                                                             tCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGata                                                                   >  2:500904/1‑96 (MQ=255)                                                                                             tCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgc                                                         <  1:523602/106‑1 (MQ=255)                                                                                             tCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGgc                                                         >  2:523602/1‑106 (MQ=255)                                                                                              cAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCAtt                       >  1:52566/1‑139 (MQ=255)                                                                                               aaTGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTc                      >  2:476145/1‑139 (MQ=255)                                                                                                 tGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                      <  1:806094/57‑1 (MQ=255)                                                                                                 tGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTa                                                                                                      >  2:806094/1‑57 (MQ=255)                                                                                                      gTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                                 >  2:903068/1‑57 (MQ=255)                                                                                                      gTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACt                                                                                                 <  1:903068/57‑1 (MQ=255)                                                                                                      gTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                               <  2:724927/75‑1 (MQ=255)                                                                                                      gTAAATTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCg                                                                               >  1:724927/1‑75 (MQ=255)                                                                                                         aaTTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCa                                       <  1:776444/112‑1 (MQ=255)                                                                                                         aaTTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCa                                       >  2:776444/1‑112 (MQ=255)                                                                                                         aaTTCCAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGGTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGc            >  2:321273/1‑139 (MQ=255)                                                                                                              cAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATg                            <  2:1519041/118‑1 (MQ=255)                                                                                                              cAATTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATg                            >  1:1519041/1‑118 (MQ=255)                                                                                                               aaTTGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAAc      <  2:52566/139‑1 (MQ=37)                                                                                                                 ttGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCACTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACg     <  1:695260/138‑1 (MQ=21)                                                                                                                 ttGTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCACTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACg     >  2:695260/1‑138 (MQ=255)                                                                                                                   gTTTATTGGATTTGCCCTTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCct  <  2:513924/139‑1 (MQ=18)                                                                                                                                    |                                                                                                                          TAAATAATCCCTTAGTACATCAGGAACCTTATAAAAAAACCCGCGTTCTGCAAACCTTTCTCGATAAAAGTGACAACCACTGAGTAAACAGATCAATGGAAGTAAATTCCAATTGTTTATTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCT  >  NC_000913/2996241‑2996493

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