Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,729,487	G→A	intergenic (‑43/+129)	gltW ← / ← rrsG	tRNA‑Glu/16S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,729,487	0	G	A	100.0%	258.1 / NA	80	intergenic (‑43/+129)	gltW/rrsG	tRNA‑Glu/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (35/45);  total (35/45)

ACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTG  >  NC_000913/2729461‑2729584                           |                                                                                                  aCGCGTTGCGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTa                         <  2:4349912/101‑1 (MQ=11) aCGCGTTGCGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTa                         <  1:4350025/101‑1 (MQ=11) aCGCGTTGCGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTa                         <  2:4350124/101‑1 (MQ=11) aCGCGTTGCGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTa                         <  1:4350230/101‑1 (MQ=11)   gcgTTGCGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTAtt                       >  1:4350556/1‑101 (MQ=17)   gcgTTGCGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTAtt                       <  1:4350033/101‑1 (MQ=17)   gcgTTGCGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTAtt                       >  1:4350555/1‑101 (MQ=17)   gcgTTGCGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTAtt                       >  1:4350554/1‑101 (MQ=17)    cgTTGCGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTc                      >  2:4350557/1‑101 (MQ=17)     gTTGCGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCa                     >  1:4350558/1‑101 (MQ=17)     gTTGCGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCa                     >  2:4350560/1‑101 (MQ=17)     gTTGCGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCa                     <  1:4350265/101‑1 (MQ=17)     gTTGCGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCa                     >  2:4350559/1‑101 (MQ=17)     gTTGCGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCa                     <  2:4350049/101‑1 (MQ=17)      ttGCGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCAc                    >  1:4350561/1‑101 (MQ=17)      ttGCGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCAc                    <  2:4350051/101‑1 (MQ=17)       tGCGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACt                   >  1:4350563/1‑101 (MQ=17)       tGCGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACt                   >  2:4350562/1‑101 (MQ=17)        gCGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtt                  >  2:4350564/1‑101 (MQ=17)        gCGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtt                  <  2:4350136/101‑1 (MQ=17)        gCGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtt                  <  1:4350043/101‑1 (MQ=17)         cGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTATACGCCTTGCTATTCACttt                 >  2:4350566/1‑101 (MQ=12)         cGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACttt                 <  2:4350066/101‑1 (MQ=17)         cGTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACttt                 >  2:4350565/1‑101 (MQ=17)          gTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtttt                >  2:4350568/1‑101 (MQ=18)          gTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtttt                <  1:4350054/101‑1 (MQ=18)          gTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtttt                <  1:4350281/101‑1 (MQ=18)          gTGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtttt                >  1:4350567/1‑101 (MQ=18)           tGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtc               >  1:4350570/1‑101 (MQ=25)           tGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtc               <  2:4350208/101‑1 (MQ=25)           tGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtc               <  1:4350233/101‑1 (MQ=25)           tGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtc               <  1:4350191/101‑1 (MQ=25)           tGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtc               <  1:4350185/101‑1 (MQ=25)           tGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtc               >  2:4350569/1‑101 (MQ=25)           tGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtc               >  2:4350571/1‑101 (MQ=25)           tGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtc               >  2:4350572/1‑101 (MQ=25)           tGATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtc               <  2:4350195/101‑1 (MQ=25)            gATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtca              >  1:4350574/1‑101 (MQ=25)            gATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtca              >  1:4350576/1‑101 (MQ=25)            gATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtca              <  2:4350100/101‑1 (MQ=25)            gATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtca              >  2:4350573/1‑101 (MQ=25)            gATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtca              >  2:4350575/1‑101 (MQ=25)            gATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtca              <  1:4350015/101‑1 (MQ=25)            gATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtca              <  2:4350016/101‑1 (MQ=25)            gATATTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtca              <  2:4350013/101‑1 (MQ=25)             atatTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACtttt                <  2:4350232/98‑1 (MQ=21)             atatTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcat             >  1:4350579/1‑101 (MQ=255)             atatTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcat             >  1:4350578/1‑101 (MQ=38)             atatTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcat             >  1:4350577/1‑101 (MQ=38)             atatTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcat             <  1:4350439/101‑1 (MQ=38)             atatTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcat             <  1:4350197/101‑1 (MQ=38)              tatTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcatc            >  1:4350581/1‑101 (MQ=255)              tatTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcatc            >  2:4350580/1‑101 (MQ=255)               atTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcatca           >  1:4350582/1‑101 (MQ=255)               atTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcatca           <  1:4350170/101‑1 (MQ=255)               atTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcatca           <  1:4350089/101‑1 (MQ=255)               atTTTCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtcatca           >  2:4350583/1‑101 (MQ=255)                ttttCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAg          <  2:4349981/101‑1 (MQ=255)                ttttCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAg          <  1:4350446/101‑1 (MQ=255)                 tttCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGa         <  1:4349992/101‑1 (MQ=255)                  ttCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACTGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGAc        <  2:4350105/101‑1 (MQ=255)                  ttCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTtc               <  1:4350065/94‑1 (MQ=25)                  ttCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGAc        <  2:4350201/101‑1 (MQ=255)                  ttCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGAc        >  1:4350584/1‑101 (MQ=255)                  ttCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGAc        <  1:4350164/101‑1 (MQ=255)                   tCGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACa       <  1:4349881/101‑1 (MQ=255)                    cGTGTAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACaa      <  2:4350367/101‑1 (MQ=255)                      tgtAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATc    >  2:4350585/1‑101 (MQ=255)                      tgtAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATc    <  1:4350103/101‑1 (MQ=255)                      tgtAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATc    <  2:4350095/101‑1 (MQ=255)                       gtAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCt   <  1:4350177/101‑1 (MQ=255)                       gtAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCt   <  1:4350110/101‑1 (MQ=255)                       gtAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCt   <  1:4350107/101‑1 (MQ=255)                       gtAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCt   <  2:4350132/101‑1 (MQ=255)                       gtAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCt   <  1:4350021/101‑1 (MQ=255)                        tAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtg  >  2:4350586/1‑101 (MQ=255)                        tAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtg  >  2:4350587/1‑101 (MQ=255)                        tAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtg  <  1:4350113/101‑1 (MQ=255)                        tAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtg  >  1:4350588/1‑101 (MQ=255)                        tAGAGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCtg  <  2:4350508/101‑1 (MQ=255)                           |                                                                                                  ACGCGTTGCGTGATATTTTCGTGTAGGGTGAGCTTTCATTAATAGAAAGCGAACGGCCTTATTCTCTTCAGCCTCACTCCCAACGCGTAAACGCCTTGCTATTCACTTTTCATCAGACAATCTG  >  NC_000913/2729461‑2729584

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