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%	179.4 / NA	60	*289* (TAG→TAA)	insF1	IS3 transposase B
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/28);  total (32/28)

CTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCC  >  NC_000913/391605‑391845                                                                                                                                       |                                                                                                           cttcGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAg                                                                                                        <  1:1149740/139‑1 (MQ=255)         aaCATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTc                                                                                                >  1:127930/1‑139 (MQ=255)              tGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGtt                                                                                           >  1:1486514/1‑139 (MQ=255)                        ggCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAg                                                                                                        >  1:459867/1‑116 (MQ=255)                        ggCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAg                                                                                                        <  2:459867/116‑1 (MQ=255)                             cGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCaaa                                                                                     <  1:1301578/130‑1 (MQ=255)                             cGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCaaa                                                                                     >  2:1301578/1‑130 (MQ=255)                                   tACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAg                                                                      >  1:1052674/1‑139 (MQ=255)                                      cGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgcc                                                                   <  2:622539/139‑1 (MQ=255)                                       gATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgccg                                                                  >  1:1221552/1‑139 (MQ=255)                                        atcatcTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgccgc                                                                 >  1:1429485/1‑139 (MQ=255)                                                  tATTAAGTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACt                                                       <  1:111468/139‑1 (MQ=255)                                                       aCTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTATTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTgcc                                                  >  1:1046973/1‑139 (MQ=255)                                                        cTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTgccg                                                 <  1:931405/139‑1 (MQ=255)                                                        cTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTgccg                                                 <  2:871172/139‑1 (MQ=255)                                                            tACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTgccgccgc                                             >  1:767656/1‑139 (MQ=255)                                                            tACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTgccgccgc                                             >  2:909777/1‑139 (MQ=255)                                                              cAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCa                                           >  2:759830/1‑139 (MQ=255)                                                              cAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCa                                           <  2:767656/139‑1 (MQ=255)                                                               aGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCAc                                          <  2:1052674/139‑1 (MQ=255)                                                                      ttGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgcc                                                                   <  2:771483/107‑1 (MQ=255)                                                                      ttGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgcc                                                                   >  1:771483/1‑107 (MQ=255)                                                                      ttGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGt                                   >  2:48070/1‑139 (MQ=255)                                                                      ttGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGt                                   >  1:1208771/1‑139 (MQ=255)                                                                           ttGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCa                                           >  1:1034095/1‑126 (MQ=255)                                                                           ttGACCGAATCTACTGGGTGGAGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCa                                           <  2:1034095/126‑1 (MQ=255)                                                                              aCCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGa                                                                           <  1:192706/91‑1 (MQ=255)                                                                              aCCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGa                                                                           >  2:192706/1‑91 (MQ=255)                                                                              aCCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAAt                                >  1:1478178/1‑134 (MQ=255)                                                                              aCCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAAt                                <  2:1478178/134‑1 (MQ=255)                                                                              aCCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACAt                           >  2:71371/1‑139 (MQ=255)                                                                                   aTCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgcc                                                                   <  1:1332454/94‑1 (MQ=255)                                                                                   aTCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGgcc                                                                   >  2:1332454/1‑94 (MQ=255)                                                                                   aTCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTgccgcc                                               >  1:297116/1‑114 (MQ=255)                                                                                   aTCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTgccgcc                                               <  2:297116/114‑1 (MQ=255)                                                                                    tCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGata                     <  2:1429485/139‑1 (MQ=255)                                                                                    tCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGata                     <  2:1208771/139‑1 (MQ=255)                                                                                    tCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGata                     <  1:71371/139‑1 (MQ=255)                                                                                    tCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGata                     <  1:759830/139‑1 (MQ=255)                                                                                             tgggGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGata                     >  2:98814/1‑130 (MQ=38)                                                                                             tgggGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGata                     <  1:98814/130‑1 (MQ=38)                                                                                             tgggGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAac            >  2:873721/1‑139 (MQ=38)                                                                                             tgggGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAac            >  1:1308764/1‑139 (MQ=255)                                                                                              ggggggCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAaca           >  1:627211/1‑139 (MQ=37)                                                                                               gggggCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAacac          >  2:1611738/1‑139 (MQ=37)                                                                                                ggggCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAg                                                                      >  2:1492423/1‑78 (MQ=255)                                                                                                ggggCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAg                                                                      <  1:1492423/78‑1 (MQ=255)                                                                                                ggggCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAg                                                                      >  1:988265/1‑78 (MQ=255)                                                                                                ggggCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAAt                                <  1:743781/116‑1 (MQ=25)                                                                                                ggggCTGATTGATCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAAt                                >  2:743781/1‑116 (MQ=37)                                                                                                gcggCTGATTGATCCTACCCACGTAACATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAg                                                                      <  2:988265/76‑1 (MQ=11)                                                                                                     tgattgatCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTgc                                                   <  1:215844/92‑1 (MQ=18)                                                                                                     tgattgatCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTgc                                                   >  2:215844/1‑92 (MQ=18)                                                                                                     tgattgatCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAAt                                <  2:406774/111‑1 (MQ=17)                                                                                                     tgattgatCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAAt                                >  1:406774/1‑111 (MQ=17)                                                                                                     tgattgatCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATcaca                            <  1:1398832/115‑1 (MQ=17)                                                                                                     tgattgatCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATcaca                            >  2:1398832/1‑115 (MQ=17)                                                                                                       attgatCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGccc  <  2:1486514/139‑1 (MQ=11)                                                                                                       attgatCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGccc  >  2:1331/1‑139 (MQ=11)                                                                                                       attgatCCTACCCACGTAATATGGACACAGGCTTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGccc  <  2:1308764/139‑1 (MQ=11)                                                                                                                                       |                                                                                                           CTTCGCTGAACATTGGTGAAGATGGCTACGTTGATACCGATCATCTGACTATTAACTCCTACAGTACTGTTGCGTTGACCGAATCTACTGGGTGGGGGGCTGATTGATCCTACCCACGTAATATGGACACAGGCCTAAGCGAGGTTCTTGTTTTCAAATTGTTCCGGACTGAGGCCGCCACACCAACTGTGCCGCCGCCACCGATTGTAATCACATTCGATATAATTAAACACCGTTGCCC  >  NC_000913/391605‑391845

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