Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,650,904	T→C	100%	pseudogene (60/693 nt)	insD1 →	IS2 transposase TnpB;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,650,904	0	T	C	100.0%	104.6 / ‑6.7	33	pseudogene (60/693 nt)	insD1	IS2 transposase TnpB;IS, phage, Tn; Transposon‑related functions; extrachromosomal; transposon related
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (16/16);  total (16/17)

AGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGTTGCTTG  >  NC_000913/1650808‑1650941                                                                                                 |                                      aGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACgcgt                                     <  2:2124881/99‑1 (MQ=255)  gCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGttt                                   >  1:1093469/1‑100 (MQ=255)  gCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGttt                                   >  1:1084710/1‑100 (MQ=255)      aCATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCg                               <  1:347835/100‑1 (MQ=255)       cATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCgg                              >  2:1081295/1‑100 (MQ=255)         tGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGAt                            >  2:2339123/1‑100 (MQ=255)         tGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGAt                            <  1:2131709/100‑1 (MQ=255)         tGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGAt                            <  1:2131978/100‑1 (MQ=255)         tGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGAt                            <  2:1919504/100‑1 (MQ=255)         tGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGAt                            >  1:2292372/1‑100 (MQ=255)         tGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGAt                            >  1:580073/1‑100 (MQ=255)          ggAAAACACTGGATGACGAACTGGCCTATGGGCGCTGTTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATc                           >  2:1372400/1‑100 (MQ=255)          ggAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATc                           >  1:2027008/1‑100 (MQ=255)           gAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCa                          <  2:569744/100‑1 (MQ=255)           gAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCa                          <  2:569591/100‑1 (MQ=255)            aaaaCACTGGATGACGAACTGGCCTATGGGCGCTGCTTCACAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCAt                         <  2:547068/100‑1 (MQ=255)             aaaCACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATg                        <  1:615087/100‑1 (MQ=255)              aaCACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATgc                       <  2:52700/100‑1 (MQ=255)               acacTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATgcg                      <  1:417371/100‑1 (MQ=255)                    ggATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGa                 >  2:2080758/1‑100 (MQ=255)                    ggATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGa                 >  1:136766/1‑100 (MQ=255)                    ggATGACGAACTGGCCTATGAGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGa                 >  2:283595/1‑100 (MQ=255)                      aTGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAAt               >  1:2264703/1‑100 (MQ=255)                         aCGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATgcg            <  2:136764/100‑1 (MQ=255)                         aCGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATgcg            >  1:1287595/1‑100 (MQ=255)                          cGAACTGGCCTATGGGCGCTGCTTCGCAGAGAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATgcgc           <  2:238451/100‑1 (MQ=255)                          cGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATgcgc           <  2:1940505/100‑1 (MQ=255)                           gAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATGCGCt          <  1:856428/100‑1 (MQ=255)                            aaCTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATGCGCTg         >  2:1010731/1‑100 (MQ=255)                              cTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATGCGCTGtt       >  2:656450/1‑100 (MQ=25)                                  ccTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGTTg      <  2:2478185/97‑1 (MQ=31)                                   cTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATGCGCTGTTGCTTg  <  1:713990/100‑1 (MQ=11)                                   cTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATGCGCTGTTGCTTg  >  1:1551656/1‑100 (MQ=11)                                   cTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATGCGCTGTTGCTTg  <  1:1288514/100‑1 (MQ=11)                                                                                                 |                                      AGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGTTGCTTG  >  NC_000913/1650808‑1650941

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑