Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	391,739	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	391,739	0	C	T	100.0%	141.5 / NA	45	*289* (TAG→TAA)	insF1	IS3 transposase B
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (25/20);  total (25/20)

GAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCC  >  NC_000913/391622‑391844                                                                                                                      |                                                                                                          gAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCa                                                                                      <  2:1078878/139‑1 (MQ=255)       ggCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGt                                                                                >  1:1167257/1‑139 (MQ=255)       ggCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGt                                                                                >  2:554947/1‑139 (MQ=255)                  tACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAg                                                                     >  1:51896/1‑139 (MQ=255)                   aCCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAg                                                                     >  1:1290989/1‑138 (MQ=255)                   aCCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAg                                                                     <  2:1290989/138‑1 (MQ=255)                     cGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgcc                                                                  <  2:600945/139‑1 (MQ=255)                       atcatcTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgccgc                                                                >  2:1119139/1‑139 (MQ=255)                       atcatcTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgccgc                                                                <  2:463075/139‑1 (MQ=255)                           tcTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGa                                                                      <  2:783643/129‑1 (MQ=255)                           tcTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGa                                                                      >  1:783643/1‑129 (MQ=255)                           tcTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCcaca                                                            <  2:51896/139‑1 (MQ=255)                           tcTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCcaca                                                            >  2:1283854/1‑139 (MQ=255)                                             cAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCa                                          >  2:503321/1‑139 (MQ=255)                                                  cTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGttt                                                                                         >  2:72204/1‑87 (MQ=255)                                                  cTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGttt                                                                                         <  1:72204/87‑1 (MQ=255)                                                     ttGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCcaca                                                            <  2:488442/113‑1 (MQ=255)                                                     ttGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCcaca                                                            >  1:488442/1‑113 (MQ=255)                                                     ttGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAAc                                                       >  1:958393/1‑118 (MQ=255)                                                     ttGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAAc                                                       <  2:958393/118‑1 (MQ=255)                                                     ttGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGt                                  >  1:1164036/1‑139 (MQ=255)                                                     ttGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGGTCTTGGTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGt                                  >  1:796554/1‑139 (MQ=255)                                                       gCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCAc                                         <  2:596098/130‑1 (MQ=255)                                                       gCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCAc                                         >  1:596098/1‑130 (MQ=255)                                                        cGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAAt                               >  1:437932/1‑139 (MQ=255)                                                             aCCGAATCTACTGGGT‑GGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgccgc                                                                >  1:371665/1‑100 (MQ=255)                                                             aCCGAATCTACTGGGT‑GGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgccgc                                                                <  2:371665/100‑1 (MQ=255)                                                              ccGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACAtt                         <  1:503321/139‑1 (MQ=255)                                                               cGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTc                        >  1:733938/1‑139 (MQ=255)                                                                   tCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGata                    <  2:130853/139‑1 (MQ=255)                                                                       cTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAAt                >  2:1265352/1‑139 (MQ=255)                                                                            tgggGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCcac                                                             <  2:179262/89‑1 (MQ=255)                                                                            tgggGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCcac                                                             >  1:179262/1‑89 (MQ=255)                                                                            tgggGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACAtt                         <  2:613515/125‑1 (MQ=255)                                                                            tgggGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACAtt                         >  1:613515/1‑125 (MQ=255)                                                                              gggggCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCaaa                                                                                    >  1:1033089/1‑64 (MQ=255)                                                                              gggggCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCaaa                                                                                    <  2:1033089/64‑1 (MQ=255)                                                                              gggggCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGa                      <  1:1275070/126‑1 (MQ=37)                                                                              gggggCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGa                      >  2:96001/1‑126 (MQ=25)                                                                              gggggCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGa                      <  1:96001/126‑1 (MQ=25)                                                                              gggggCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGa                      >  2:1275070/1‑126 (MQ=37)                                                                               ggggCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCATATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCa                                          <  2:915163/105‑1 (MQ=255)                                                                               ggggCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACAcc        >  1:612918/1‑139 (MQ=25)                                                                               ggggCTGATTGATCCAACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCa                                          >  1:915163/1‑105 (MQ=255)                                                                                     gattgatCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGcc  <  2:1167257/139‑1 (MQ=14)                                                                                                                      |                                                                                                          GAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCC  >  NC_000913/391622‑391844

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