Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,503,643	T→C	intergenic (‑212/‑415)	insG ← / → yjhB	IS4 transposase/putative MFS transporter, membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,503,643	0	T	C	100.0%	207.5 / NA	54	intergenic (‑212/‑415)	insG/yjhB	IS4 transposase/putative MFS transporter, membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (34/20);  total (34/20)

ATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTC  >  NC_000913/4503495‑4503784                                                                                                                                                     |                                                                                                                                              aTCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGcat                                                                                                                                             >  1:626597/1‑151 (MQ=255)   ccTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGGTTTTTGATCTGGTATAACAGGTATAAAGGCATAc                                                                                                                                           >  1:351318/1‑151 (MQ=255)      tAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATAc                                                                                                                                           >  2:414617/1‑148 (MQ=255)      tAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAg                                                                                                                                         <  2:442779/150‑1 (MQ=255)          tGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGaaa                                                                                                                                      >  2:543804/1‑149 (MQ=255)           gATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGaaaa                                                                                                                                     <  2:362558/149‑1 (MQ=255)           gATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGc                                                                                                                                   >  1:348019/1‑151 (MQ=255)            aTCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCa                                                                                                                                  <  2:354692/151‑1 (MQ=255)              cGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCaaa                                                                                                                                <  2:416065/151‑1 (MQ=255)               ggCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCaa                                                                                                                                 >  2:108087/1‑149 (MQ=255)                   ttAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAAt                                                                                                                           >  1:271400/1‑151 (MQ=255)                       tCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTg                                                                                                                       >  2:483777/1‑151 (MQ=255)                       tCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTg                                                                                                                       >  2:301226/1‑151 (MQ=255)                        cTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGc                                                                                                                      >  2:319286/1‑151 (MQ=255)                           ggTTCTGGGGTTTGTAACACAGTATCAGCTACAAAACTATGCTCAATTTGTGACATCAGTGACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGAATAAAGGCATACAGAACAGCAAAGAAATACTGCa                                                                                                                     >  1:60149/1‑149 (MQ=255)                                 ggTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGaaa                                                                                                              >  2:7023/1‑150 (MQ=255)                                  gtgtTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAAc                                                                                                            <  1:168132/151‑1 (MQ=255)                                   tgtTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAAc                                                                                                            >  1:315142/1‑150 (MQ=255)                                    gtTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAAc                                                                                                            >  1:121852/1‑149 (MQ=255)                                    gtTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAg                                                                                                          <  1:322276/151‑1 (MQ=255)                                           aCAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAAc                                                                                                   <  2:640390/151‑1 (MQ=255)                                           aCAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAAc                                                                                                   >  1:96267/1‑151 (MQ=255)                                                  aTCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGc                                                                                             <  1:567337/150‑1 (MQ=255)                                                         aCAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGta                                                                                      <  1:543799/150‑1 (MQ=255)                                                          cAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGtaat                                                                                    >  1:340058/1‑151 (MQ=255)                                                                tatGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACctc                                                                              >  2:465793/1‑151 (MQ=255)                                                                  tGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCtctc                                                                            >  2:499111/1‑151 (MQ=255)                                                                   gCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCa                                                                           <  2:335007/151‑1 (MQ=255)                                                                       aaTTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCa                                                                           <  1:113374/147‑1 (MQ=255)                                                                       aaTTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCa                                                                           >  2:113377/1‑147 (MQ=255)                                                                          ttGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTc                                                                    >  1:541439/1‑151 (MQ=255)                                                                           tgtgACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTcc                                                                   >  2:176280/1‑151 (MQ=255)                                                                           tgtgACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTc                                                                    <  1:492574/150‑1 (MQ=255)                                                                                cATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCAtt                                                              <  2:485787/151‑1 (MQ=255)                                                                                 aTCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCAttt                                                             >  1:233030/1‑151 (MQ=255)                                                                                       aaCAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGa                                                        >  1:9075/1‑150 (MQ=255)                                                                                         cAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAgg                                                     >  2:155832/1‑151 (MQ=255)                                                                                           aaaCGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAggg                                                    >  1:231363/1‑150 (MQ=255)                                                                                           aaaCGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGa                                                   <  2:569860/151‑1 (MQ=255)                                                                                               gcgTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGttt                                               <  1:517875/151‑1 (MQ=255)                                                                                               gcgTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGttt                                               >  2:417572/1‑151 (MQ=255)                                                                                                 gTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGtttt                                              >  2:481454/1‑150 (MQ=255)                                                                                                      gTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCaa                                        <  2:315144/151‑1 (MQ=255)                                                                                                           gtgGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTg                                    >  1:536703/1‑150 (MQ=255)                                                                                                               aTTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACgg                               <  1:646910/151‑1 (MQ=255)                                                                                                                    ttGTTTTTTGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTAc                          >  2:557846/1‑151 (MQ=255)                                                                                                                       ttttttGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGa                        >  2:523382/1‑150 (MQ=255)                                                                                                                         ttttGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGaaa                      >  2:477758/1‑150 (MQ=255)                                                                                                                          tttGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAAt                    >  2:457185/1‑151 (MQ=255)                                                                                                                            tGATCTGGTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATg                   >  1:667831/1‑150 (MQ=255)                                                                                                                                   gTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAAt           >  1:634568/1‑151 (MQ=255)                                                                                                                                   gTATAACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAAt           <  1:176279/151‑1 (MQ=255)                                                                                                                                     ataACAGGTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATg          <  2:340061/150‑1 (MQ=255)                                                                                                                                           ggTATAAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGt   <  2:666577/151‑1 (MQ=255)                                                                                                                                             tataAAGGCATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTc  >  2:13085/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                              ATCCTTAACTGATCGGCATTAATCTTGGTTCTGGTGTTTGTAACAAACTATCAGCTACAAAAATATGCTCAATTTGTGACATCAGTAACAAAACGCGTTTTGTTATGTGGATTGCTTGTTTTTTGATCTGGTATAACAGGTATAAAGGTATACAGAAAAGCAAAGAAATACTGCAAAGGAAAACAGCTATAACGTAAGCTAAAGTAATAACCTCTCAGTCTTTCCTCATTTGACGAAGGGAGTTTTATTCAACCTGAACGGACTACGAAAATGAGCACAATGAATAAGTC  >  NC_000913/4503495‑4503784

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