Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,650,904	T→C	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	97.2%	108.4 / ‑4.0	36	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 (17/18);  total (17/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.25e-01

AGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGTTGCTTG  >  NC_000913/1650808‑1650941                                                                                                 |                                      aGCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGtt                                    >  1:1246328/1‑100 (MQ=255)  gCAGACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGttt                                   <  1:1670459/100‑1 (MQ=255)     gACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTAcc                                <  1:116461/100‑1 (MQ=255)     gACATGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTGTAcc                                >  1:855245/1‑100 (MQ=255)        aTGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGa                             >  1:1696172/1‑100 (MQ=255)         tGGAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGAt                            <  2:923375/100‑1 (MQ=255)          ggAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATc                           >  2:103118/1‑100 (MQ=255)           gAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCa                          >  1:1202186/1‑100 (MQ=255)           gAAAACACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCa                          <  1:774416/100‑1 (MQ=255)            aaaaCACTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCAt                         <  1:1073217/100‑1 (MQ=255)                 acTGGATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGcc                    <  1:897008/100‑1 (MQ=255)                    ggATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGa                 >  2:1248238/1‑100 (MQ=255)                    ggATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGa                 >  1:1540644/1‑100 (MQ=255)                    ggATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGa                 >  1:766623/1‑100 (MQ=255)                    ggATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGa                 >  1:1899951/1‑100 (MQ=255)                    ggATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACCTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGa                 >  1:298327/1‑100 (MQ=255)                     gATGACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGaa                >  1:1297873/1‑100 (MQ=255)                        gACGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATgc             >  1:670458/1‑100 (MQ=255)                         aCGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATgcg            <  2:1389458/100‑1 (MQ=255)                          cGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATgcgc           <  2:615696/100‑1 (MQ=255)                          cGAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATgcgc           <  2:1686603/100‑1 (MQ=255)                          cGAACGGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATgcgc           <  1:693038/100‑1 (MQ=255)                           gAACTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATGCGCt          >  2:419797/1‑100 (MQ=255)                            aaCTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATGCGCt          <  2:255530/99‑1 (MQ=255)                            aaCTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATGCGCTg         >  2:412501/1‑100 (MQ=255)                             aCTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATGCGCTGt        <  2:64547/100‑1 (MQ=255)                             aCTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAAAATGCGCTGt        >  1:64245/1‑100 (MQ=255)                              cTGGCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGTGTTTACCGGATCATGCGCCAGAATGCGCTGtt       <  2:1503230/100‑1 (MQ=35)                                ggCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATGCGCTGTTGc     <  2:1905331/100‑1 (MQ=18)                                ggCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATGCGCTGTTGc     >  2:72496/1‑100 (MQ=18)                                 gCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATGCGCTGTTGCt    <  1:660757/100‑1 (MQ=17)                                 gCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATGCGCTGTTGCt    <  2:1077393/100‑1 (MQ=17)                                 gCCTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATGCGCTGTTGCt    <  1:980112/100‑1 (MQ=17)                                  ccTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATTCCAAACGCGTTTACCGGATCATGCGCCAGAATGCGCTGTTGCtt   >  2:1593385/1‑100 (MQ=9)                                  ccTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATGCGCTGTTGCtt   <  2:1502194/100‑1 (MQ=14)                                   cTATGGGCGCTGCTTCGCAGACAGGCAGAACTTGATGGTATGCCTGCGATCAATGCCAAACGCGTTTACCGGATCATGCGCCAGAATGCGCTGTTGCTTg  <  2:666329/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: ‑