Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,593,254	2 bp→GC	coding (59‑60/807 nt)	yjjM ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,593,254	0	C	G	100.0%	262.8 / NA	98	A20G (GCG→GGC)	yjjM	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base G (50/48);  total (50/48)
*	NC_000913	4,593,255	0	G	C	100.0%	262.8 / NA	98	A20G (GCG→GGC)	yjjM	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (50/48);  total (50/48)

CAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACCGCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTT  >  NC_000913/4593186‑4593330                                                                     ||                                                                            cAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCAtt                                                                       <  1:2969890/76‑1 (MQ=255)  aGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTc                                                                      <  1:1173718/76‑1 (MQ=255)    gCGAACCCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTcgc                                                                    <  1:2009200/76‑1 (MQ=255)    gCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTcgc                                                                    >  1:1948075/1‑76 (MQ=255)     cGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCCTTcgcg                                                                   >  2:4579149/1‑76 (MQ=255)     cGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTcgcg                                                                   >  1:3383496/1‑76 (MQ=255)      gAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGt                                                                  >  2:3549421/1‑76 (MQ=255)       aagccGTCATCATGGTCAGGTTTACGCGCAATACCATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTa                                                                 <  2:4017749/73‑1 (MQ=255)        aaCCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTaa                                                                <  2:1773677/76‑1 (MQ=255)         aCCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAg                                                               <  1:869130/76‑1 (MQ=255)          ccGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTacgg                                                              >  1:584114/1‑73 (MQ=255)          ccGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAgg                                                              >  1:579238/1‑75 (MQ=255)          ccGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCATCAAGACGCCGCATTCGCGTAAGt                                                              <  2:2879956/76‑1 (MQ=255)            gTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTgg                                                            >  1:4479950/1‑76 (MQ=255)              catcatGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCt                                                          >  2:2232785/1‑76 (MQ=255)               atcatGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTg                                                         >  1:4797391/1‑76 (MQ=255)               agcatGGTCAGGTTGACGCGCAGTGACAGCGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTg                                                         <  2:4230443/74‑1 (MQ=255)                tcatGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTga                                                        <  1:3281793/76‑1 (MQ=255)                 catGGTCAGGTTTACGCTCAAGAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTgag                                                       <  1:2115621/76‑1 (MQ=255)                  atGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTgaga                                                      >  2:2779795/1‑76 (MQ=255)                   tGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGaa                                                     >  1:2438678/1‑76 (MQ=255)                     gTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAata                                                   <  2:3864409/76‑1 (MQ=255)                     gTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAata                                                   <  1:3624414/76‑1 (MQ=255)                      tCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAatat                                                  <  1:414506/76‑1 (MQ=255)                      tCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAatat                                                  <  2:901043/76‑1 (MQ=255)                      tCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAatat                                                  >  1:3176411/1‑76 (MQ=255)                      tCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAatat                                                  <  2:4500933/76‑1 (MQ=255)                      tCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAatat                                                  <  2:1830778/76‑1 (MQ=255)                       cAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATg                                                 <  1:869828/76‑1 (MQ=255)                        aGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGa                                                >  1:1520401/1‑76 (MQ=255)                        aGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGa                                                <  1:4805967/76‑1 (MQ=255)                         ggTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGAc                                               <  1:738983/76‑1 (MQ=255)                         ggTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAACTCGCTGAGAACATTAc                                               >  2:1006107/1‑76 (MQ=255)                          gTTTGCGAGAAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACg                                              <  2:1347021/76‑1 (MQ=255)                          gTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACg                                              >  1:4727930/1‑76 (MQ=255)                            ttACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCa                                            >  1:4855949/1‑76 (MQ=255)                             tACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCAc                                           <  1:1273255/76‑1 (MQ=255)                             tACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCAc                                           >  1:906626/1‑76 (MQ=255)                              aCGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACt                                          >  1:1907061/1‑76 (MQ=255)                              aCGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACt                                          <  1:3705514/76‑1 (MQ=255)                              aCGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACt                                          >  2:2635596/1‑76 (MQ=255)                               cgcgCAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACtg                                         >  2:1662460/1‑76 (MQ=255)                                gcgcAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACtgt                                        >  1:1584993/1‑76 (MQ=255)                                gcgcAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACtgg                                        >  2:4361118/1‑75 (MQ=255)                                gcgcAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACtgg                                        >  1:9279/1‑75 (MQ=255)                                 cgcAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACtgtg                                       <  2:2819627/76‑1 (MQ=255)                                  gcAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTgtgg                                      >  2:4087919/1‑76 (MQ=255)                                  gcAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTgtgg                                      >  2:4408952/1‑76 (MQ=255)                                  gcAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTgtgg                                      <  2:1312751/76‑1 (MQ=255)                                  gcAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTgtgg                                      >  2:3095421/1‑76 (MQ=255)                                   cAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTgtggt                                     <  1:4087919/76‑1 (MQ=255)                                   cAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTgtggt                                     <  2:4679313/76‑1 (MQ=255)                                   cAATAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTgtggg                                     >  2:2284941/1‑75 (MQ=255)                                   cAATAACATAGCCGTTGCCAACCGGAGCCGAGACGCCGCAGTCGCGTAAGTGGCTGAGAGTATGGCGCACTgtggt                                     >  2:3446832/1‑76 (MQ=255)                                    aaTAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTgtggtg                                    <  1:3930098/76‑1 (MQ=255)                                    aaTAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTgtggtg                                    >  2:2809032/1‑76 (MQ=255)                                     aTAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGc                                   <  2:3253827/76‑1 (MQ=255)                                     aTAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGc                                   >  2:2442227/1‑76 (MQ=255)                                      tAACATAGTCGTTTCCAGCCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCg                                  <  2:4180238/76‑1 (MQ=255)                                      tAACATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCg                                  >  2:642191/1‑76 (MQ=255)                                       aaCATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCgg                                 >  1:3415546/1‑76 (MQ=255)                                        aCATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGc                                >  1:3308022/1‑76 (MQ=255)                                         cATAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGcg                               >  1:338721/1‑75 (MQ=255)                                          aTAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTg                              >  1:4743764/1‑76 (MQ=255)                                           tAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGa                             >  1:856340/1‑76 (MQ=255)                                           tAGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGa                             <  2:3997640/76‑1 (MQ=255)                                           tAGTCGTGTCCAACCTGCGGCAAGACGCCGCAGGCGCGTAAGTGGCTGAGAATATGACGCACTGTGGGGCGGCTGa                             <  2:4033741/76‑1 (MQ=255)                                            aGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAt                            >  1:3663380/1‑76 (MQ=255)                                            aGTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAt                            >  1:3687741/1‑76 (MQ=255)                                             gTCGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtg                           >  2:3446170/1‑76 (MQ=255)                                              tCGTTTACAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgt                          <  1:24169/76‑1 (MQ=255)                                               cGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgtt                         >  2:2151959/1‑76 (MQ=255)                                               cGTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgtt                         <  2:121265/76‑1 (MQ=255)                                                gTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgttg                        >  2:2636238/1‑76 (MQ=255)                                                gTTTCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgttg                        <  2:2332017/76‑1 (MQ=255)                                                gTTGCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgttg                        <  2:2673418/76‑1 (MQ=255)                                                 tttCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGAtgctgt                       >  2:2640354/1‑76 (MQ=255)                                                   tCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTAc                     >  2:2139111/1‑76 (MQ=255)                                                   tCCAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTAc                     >  2:1703138/1‑76 (MQ=255)                                                    ccAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACa                    <  2:1760192/76‑1 (MQ=255)                                                    ccAACCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACa                    <  1:2779795/76‑1 (MQ=255)                                                     cAACCTGCGTCAAGCCGCCGCATTCGCGTAAGTGGCTGCGAATATGACGCACTGTGGTGCGGCTGATGTTGTACAt                   <  1:3503298/76‑1 (MQ=255)                                                     cAACCCGCGGCAAGACGCCGCATTCGCGTAAGTGGGTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACAt                   <  2:299313/76‑1 (MQ=255)                                                      aaCCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACAtt                  >  1:4337123/1‑76 (MQ=255)                                                      aaCCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACAtt                  <  2:4657052/76‑1 (MQ=255)                                                       aCCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACAttt                 <  1:4387923/76‑1 (MQ=255)                                                       aCCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACAttt                 <  1:2767610/76‑1 (MQ=255)                                                       aCCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACAttt                 <  2:36861/76‑1 (MQ=255)                                                       aCCTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACAttt                 >  2:1715885/1‑76 (MQ=255)                                                         cTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTcc               >  1:720868/1‑76 (MQ=255)                                                         cTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTcc               >  1:3340948/1‑76 (MQ=255)                                                         cTGCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTcc               >  1:3345334/1‑76 (MQ=255)                                                           gCGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGc             <  1:344656/76‑1 (MQ=255)                                                            cGTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGcc            <  1:1225035/76‑1 (MQ=255)                                                             gTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCa           <  1:2635596/76‑1 (MQ=255)                                                             gTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCa           <  2:3237951/76‑1 (MQ=255)                                                             gTCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCa           <  2:1850915/76‑1 (MQ=255)                                                              tCAAGACGCCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAg          >  2:2764207/1‑76 (MQ=255)                                                                      gcgcATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACttt  <  1:1537273/76‑1 (MQ=255)                                                                     ||                                                                            CAGGCGAAACCGTCATCATGGTCAGGTTTACGCGCAATAACATAGTCGTTTCCAACCTGCGTCAAGACCGCGCATTCGCGTAAGTGGCTGAGAATATGACGCACTGTGGTGCGGCTGATGTTGTACATTTCCGCCAGCGCACTTT  >  NC_000913/4593186‑4593330

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