Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,581,598	A→G	R38R (CGT→CGC)	b1499 ←	putative ARAC‑type regulatory protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,581,598	0	A	G	100.0%	267.2 / NA	86	R38R (CGT→CGC)	b1499	putative ARAC‑type regulatory protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (37/49);  total (37/49)

TGGATGTTAGTGCATTGTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTACGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGTTAGCGTTGAAGGCATGATGAATAAA  >  NC_000913/1581509‑1581687                                                                                          |                                                                                          tGGATGTTAGTGCATTGTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGa                                                                                        >  1:3322265/1‑93 (MQ=255)  ggATGTTAGTGCATTGTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGaa                                                                                       <  2:2717771/93‑1 (MQ=255)   gATGTTAGTGCATTGTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAAc                                                                                      <  2:1568246/93‑1 (MQ=255)    aTGTTAGTGCATTGTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACa                                                                                     >  2:2450771/1‑93 (MQ=255)      gTTAGTGCATTGTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGc                                                                                   <  1:1871396/93‑1 (MQ=255)       ttAGTGCATTGTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCg                                                                                  >  2:2681571/1‑93 (MQ=255)          gTGCATTGTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTa                                                                               <  1:350636/93‑1 (MQ=255)          gTGCATTGTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTa                                                                               <  2:2306929/93‑1 (MQ=255)           tGCATTGTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTAt                                                                              <  2:4286952/93‑1 (MQ=255)           tGCATTGTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTAt                                                                              <  2:4089503/93‑1 (MQ=255)           tGCATTGTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTAt                                                                              <  1:4567958/93‑1 (MQ=255)           tGCATTGTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTAt                                                                              <  1:3566050/93‑1 (MQ=255)              aTTGTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTa                                                                           >  1:3100989/1‑93 (MQ=255)              aTTGTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTa                                                                           >  1:4164643/1‑93 (MQ=255)              aTTGTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTa                                                                           <  2:2284198/93‑1 (MQ=255)              aTTGTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTa                                                                           >  2:1950040/1‑93 (MQ=255)               ttGTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTAc                                                                          <  1:347364/93‑1 (MQ=255)                 gTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCa                                                                        >  2:2289917/1‑93 (MQ=255)                 gTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCa                                                                        <  1:768125/93‑1 (MQ=255)                 gTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCa                                                                        <  2:510705/93‑1 (MQ=255)                   aGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATc                                                                      >  2:2881812/1‑93 (MQ=255)                    gATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCa                                                                     >  2:1710240/1‑93 (MQ=255)                    gATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCa                                                                     >  2:966791/1‑93 (MQ=255)                     atatTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAg                                                                    <  2:608132/93‑1 (MQ=255)                       atTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGa                                                                  <  2:3056018/93‑1 (MQ=255)                        ttttttAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGAt                                                                 >  1:3099155/1‑93 (MQ=255)                        ttttttAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGAt                                                                 >  1:3671014/1‑93 (MQ=255)                         tttttAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGAtc                                                                <  2:4633891/93‑1 (MQ=255)                          ttttAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGAtct                                                               <  1:1149988/93‑1 (MQ=255)                          ttttAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGAtct                                                               >  2:4064783/1‑93 (MQ=255)                          ttttAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGAtct                                                               <  1:1878095/93‑1 (MQ=255)                             tAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGc                                                            <  1:1568170/93‑1 (MQ=255)                               aTGGTATTTTCAGTAAATCCGATCTCTTTAAAATGACGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCg                                                          >  1:2401402/1‑93 (MQ=255)                                tGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGt                                                         >  2:1326837/1‑93 (MQ=255)                                tGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGt                                                         >  2:2479236/1‑93 (MQ=255)                                 ggTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTc                                                        >  1:4285097/1‑93 (MQ=255)                                   tATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCag                                                      >  1:4416636/1‑93 (MQ=255)                                    aTTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCaga                                                     <  2:4678081/93‑1 (MQ=255)                                      tttCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAgaga                                                   <  2:3310761/93‑1 (MQ=255)                                       ttCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGaa                                                  <  2:2689614/93‑1 (MQ=255)                                         cAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAgg                                                <  2:2821444/93‑1 (MQ=255)                                          aGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGc                                               >  2:2443105/1‑93 (MQ=255)                                             aaaTCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTa                                            >  1:1693646/1‑93 (MQ=255)                                              aaTCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTaa                                           <  1:272206/93‑1 (MQ=255)                                              aaTCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTaa                                           <  1:482000/93‑1 (MQ=255)                                               aTCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAAt                                          <  1:1315417/93‑1 (MQ=255)                                               aTCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAAt                                          <  2:1604974/93‑1 (MQ=255)                                                tCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATc                                         >  1:1944473/1‑93 (MQ=255)                                                tCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATc                                         <  1:3896060/93‑1 (MQ=255)                                                  cGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCtt                                       <  2:2863031/93‑1 (MQ=255)                                                  cGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCtt                                       <  1:283538/93‑1 (MQ=255)                                                  cGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCtt                                       >  2:1485562/1‑93 (MQ=255)                                                   gATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCttt                                      <  1:3131849/93‑1 (MQ=255)                                                    aTCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTa                                     >  1:4127009/1‑93 (MQ=255)                                                    aTCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTa                                     >  2:2055197/1‑93 (MQ=255)                                                    aTCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTa                                     >  2:2968121/1‑93 (MQ=255)                                                    aTCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTa                                     >  1:209347/1‑93 (MQ=255)                                                       tctTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCt                                  >  1:1435649/1‑93 (MQ=255)                                                        ctTTAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTa                                 >  2:3630896/1‑93 (MQ=255)                                                         tttAAAATGGCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTaa                                >  2:3068962/1‑93 (MQ=255)                                                          ttAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTaaa                               <  2:90396/93‑1 (MQ=255)                                                           tAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTaaaa                              <  2:4083901/93‑1 (MQ=255)                                                           tAAAATGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTaaaa                              <  2:3657897/93‑1 (MQ=255)                                                            aaaaTGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAAt                             <  2:135264/93‑1 (MQ=255)                                                             aaaTGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATg                            <  1:797521/93‑1 (MQ=255)                                                             aaaTGCCGCTCATAAGGTTCAAAATGCGCGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATg                            >  1:590192/1‑93 (MQ=255)                                                              aaTGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGt                           <  1:135359/93‑1 (MQ=255)                                                               aTGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGtt                          >  1:1651536/1‑93 (MQ=255)                                                                tGCCGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCAGTCAGAGAAGGCGTAATCTTTATCTAAAATGTTa                         <  1:2644195/93‑1 (MQ=255)                                                                   cGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGTTAGCg                      <  1:1036755/93‑1 (MQ=255)                                                                   cGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGTTAGCg                      >  2:1499073/1‑93 (MQ=255)                                                                   cGCTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGTTAGCg                      >  1:2796392/1‑93 (MQ=255)                                                                     cTCATGAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGTTAGCGtt                    <  2:1988520/93‑1 (MQ=255)                                                                     cTCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGTTAGCGtt                    <  2:2172197/93‑1 (MQ=255)                                                                      tCATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGTTAGCGTTg                   >  2:453515/1‑93 (MQ=255)                                                                       cATAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGTTAGCGTTGa                  <  1:186885/93‑1 (MQ=255)                                                                        aTAAGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGTTAGCGTTGaa                 <  1:2473408/93‑1 (MQ=255)                                                                         aaaGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGTTAGCGTTGAAg                <  2:3228764/92‑1 (MQ=255)                                                                           aGGTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGTTAGCGTTGAAGGc              <  1:77880/93‑1 (MQ=255)                                                                             gTTCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGTTAGCGTTGAAGGCat            <  1:2768146/93‑1 (MQ=255)                                                                              ttCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGTTAGCGTTGAAGGCatg           >  1:2557679/1‑93 (MQ=255)                                                                              ttCAAAATGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGTTAGCGTTGAAGGCatg           >  2:1844552/1‑93 (MQ=255)                                                                                 aaaaTGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGTTAGCGTTGAAGGCatgatg        <  1:2194813/93‑1 (MQ=255)                                                                                  aaaTGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGTTAGCGTTGAAGGCAtgatga       <  1:4407809/93‑1 (MQ=255)                                                                                   aaTGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGTTAGCGTTGAAGGCATGATGaa      <  1:2068708/93‑1 (MQ=255)                                                                                   aaTGCGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGTTAGCGTTGAAGGCATGATGaa      <  2:1792204/93‑1 (MQ=255)                                                                                       cGTGCGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGTTAGCGTTGAAGGCATGATGaataaa  >  1:738591/1‑93 (MQ=255)                                                                                          |                                                                                          TGGATGTTAGTGCATTGTAGATATTTTTTAATGGTATTTTCAGTAAATCCGATCTCTTTAAAATGCCGCTCATAAGGTTCAAAATGCGTACGAACAGCGTTATCTACCATCAGGATCTCGCCGTCAGAGAAGGCGTAATCTTTATCTAAAATGTTAGCGTTGAAGGCATGATGAATAAA  >  NC_000913/1581509‑1581687

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