Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,812,859	Δ1 bp	intergenic (‑105/+46)	ygaG ← / ← gshA	hypothetical protein/gamma‑glutamate‑cysteine ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,812,858	0	C	.	100.0%	282.9 / NA	62	intergenic (‑104/+47)	ygaG/gshA	hypothetical protein/gamma‑glutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (36/26);  total (36/26)

TTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACC  >  NC_000913/2812784‑2812929                                                                           |                                                                        ttttCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTca                                                                        <  2:2599272/76‑2 (MQ=255) ttttCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTca                                                                        <  2:2622342/76‑2 (MQ=255)  tttCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTcag                                                                       >  2:2060074/1‑74 (MQ=255)  tttCTGAACTCGTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATGTGTTATCATCATCCCTGAATTcag                                                                       >  2:2483143/1‑74 (MQ=255)   ttCTGAACTCGGTCTTCCCATGCGAGTCTGAGTATATGAAAGACGCGCACTTGTTAGCATCATCCCTGAATTcaga                                                                      >  2:1521794/1‑73 (MQ=255)     cTGAACTCTTTCTTCCCGGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAgaga                                                                   >  1:1015133/1‑76 (MQ=255)      tGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGAt                                                                  <  1:1564074/76‑1 (MQ=255)       gAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATg                                                                 <  2:2082618/76‑1 (MQ=255)        aaCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGa                                                                >  1:1385106/1‑76 (MQ=255)        aaCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGa                                                                >  1:610803/1‑76 (MQ=255)         aCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGaa                                                               >  1:1086919/1‑76 (MQ=255)          ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGaaa                                                              >  2:3409374/1‑76 (MQ=255)          ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGaaa                                                              <  1:2490181/76‑1 (MQ=255)           tctttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAAt                                                             <  1:2628198/76‑1 (MQ=255)           tctttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAAt                                                             <  1:3113180/76‑1 (MQ=255)            ctttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAAtt                                                            >  2:3553203/1‑76 (MQ=255)              ttcttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAAtttt                                                          >  1:1678561/1‑76 (MQ=255)               tcttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTg                                                         >  1:1353949/1‑76 (MQ=255)               tcttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTg                                                         >  1:4016463/1‑76 (MQ=255)               tcttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTg                                                         <  1:3684224/76‑1 (MQ=255)                cttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTgg                                                        >  1:235592/1‑76 (MQ=255)                cttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTgg                                                        >  2:286796/1‑76 (MQ=255)                 ttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGc                                                       >  2:1925041/1‑76 (MQ=255)                   cccAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCa                                                     >  1:1788023/1‑76 (MQ=255)                   cccAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCa                                                     <  1:1871006/76‑1 (MQ=255)                   cccAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCACCATCCCTGAATT‑CAGAGATGAAATTTTGGCCa                                                     >  2:3842981/1‑76 (MQ=255)                     cAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACt                                                   >  1:182204/1‑76 (MQ=255)                      aGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTc                                                  >  2:2762848/1‑76 (MQ=255)                      aGGCGAGTCTGAGGATATGAAAGACGCGCAGTTGTTATCATCAACCCTGAATT‑CACAGATGAAATTTTGGGCAcac                                                  >  2:2401611/1‑74 (MQ=255)                        gCGAGTCTGAGTATATGAAAGATGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCAc                                                <  2:4799092/76‑1 (MQ=255)                        gCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCAc                                                >  1:652207/1‑76 (MQ=255)                        gCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCAc                                                <  1:4780783/76‑1 (MQ=255)                         cGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACg                                               <  1:1644426/76‑1 (MQ=255)                          gAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGa                                              >  1:4035367/1‑76 (MQ=255)                           aGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAg                                             <  1:4713467/76‑1 (MQ=255)                            gTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGACATTTTGGCCACTCACGAGt                                            >  2:1120681/1‑76 (MQ=255)                            gTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGt                                            >  2:1082111/1‑76 (MQ=255)                                gAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGcc                                        >  1:63976/1‑76 (MQ=255)                                 aGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCt                                       >  2:3749049/1‑76 (MQ=255)                                  gTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCtt                                      <  2:1831606/76‑1 (MQ=255)                                    atatGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCtttt                                    <  1:2629343/76‑1 (MQ=255)                                     tatGAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAAATTTGGCCACTCACGAGTGGCCttttt                                   <  2:3460718/76‑1 (MQ=255)                                        gcaaGACGCGCCTTTGTTATCATCATCCCTGCATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctt                                >  2:2510080/3‑76 (MQ=255)                                        gAAAGACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctt                                >  2:2485393/1‑76 (MQ=255)                                            gACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTtttcttttct                            >  2:743148/1‑76 (MQ=255)                                            gACGCGCATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTtttcttttct                            >  2:138520/1‑76 (MQ=255)                                                cgcATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCa                        <  2:1136541/76‑1 (MQ=255)                                                cgcATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCa                        <  2:431506/76‑1 (MQ=255)                                                cgcATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCa                        <  1:2060075/76‑1 (MQ=255)                                                 gcATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAg                       >  1:1276182/1‑76 (MQ=255)                                                 gcATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAg                       <  1:1782421/76‑1 (MQ=255)                                                 gcATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAg                       >  2:3546545/1‑76 (MQ=255)                                                  cATTTGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAgg                      >  1:236360/1‑76 (MQ=255)                                                    tttGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCg                    >  1:3319401/1‑76 (MQ=255)                                                     ttGTTATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCgt                   >  1:4062960/1‑76 (MQ=255)                                                        ttATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGtt                <  2:3138735/76‑1 (MQ=255)                                                         tATCATCATCCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGttt               <  2:2484674/76‑1 (MQ=255)                                                          atcatcatcCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGtttt              <  2:2595007/76‑1 (MQ=255)                                                              tcatcCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttct                                 <  1:3502097/53‑1 (MQ=255)                                                              tcatcCCTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttct                                 >  2:3502097/1‑53 (MQ=255)                                                                 tcCCTGTATT‑CAGAGAGTAAATTTTGGCGACTCGCGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGcc       <  1:3742296/76‑1 (MQ=255)                                                                   ccTGTATT‑CAGAGATGAAATTTTGGCCACTCACGGGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCcac     <  2:2465282/76‑1 (MQ=255)                                                                   ccTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCcac     >  2:4810858/1‑76 (MQ=255)                                                                    cTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCcacc    >  2:200368/1‑75 (MQ=255)                                                                    cTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCcaca    <  2:3149739/76‑1 (MQ=255)                                                                    cTGAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCcaca    <  2:3140388/76‑1 (MQ=255)                                                                      gAATT‑CAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACAcc  >  2:2766972/1‑76 (MQ=255)                                                                           |                                                                        TTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACC  >  NC_000913/2812784‑2812929

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