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_0009132,812,8580C.98.9% 386.7 / ‑4.3 87intergenic (‑104/+47)ygaG/gshAhypothetical protein/gamma‑glutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (51/35);  minor base G (0/1);  total (51/36)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACC  >  NC_000913/2812783‑2812929
                                                                           |                                                                       
aTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTc                                                                         >  1:3562205/1‑76 (MQ=255)
aTTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCACCCCTGAATTc                                                                         >  2:904303/1‑76 (MQ=255)
 ttttCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTca                                                                        <  2:917322/76‑2 (MQ=255)
  tttCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTcag                                                                       >  1:3871720/1‑74 (MQ=255)
  tttCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTcag                                                                       <  1:4282874/76‑3 (MQ=255)
  tttCTGAACTCTGTCGTCCCAGGCGAGTCGTAGTATATGAAAGACGCTCATTTGTTATCATCATCCCTGAATTcag                                                                       <  1:4291708/76‑3 (MQ=255)
   ttCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTcaga                                                                      >  1:4505270/1‑73 (MQ=255)
     cTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAgaga                                                                   >  2:5631675/1‑76 (MQ=255)
     cTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAgaga                                                                   >  2:1994337/1‑76 (MQ=255)
      tGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGAt                                                                  >  2:4649942/1‑76 (MQ=255)
        caCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGa                                                                >  2:4654683/2‑76 (MQ=255)
        aaCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGa                                                                >  2:5114022/1‑76 (MQ=255)
        aaCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGa                                                                <  1:564424/76‑1 (MQ=255)
         aCTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGaa                                                               >  2:4806203/1‑76 (MQ=255)
          ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGaaa                                                              >  2:5438652/1‑76 (MQ=255)
          ctctTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACACGCATTTGTTATCATCATCCCTGAATTCAGAGATGaaa                                                              >  1:2660556/1‑76 (MQ=255)
           tctttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAAt                                                             <  1:329561/76‑1 (MQ=255)
            ctttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGa                                                                >  2:4552344/1‑72 (MQ=255)
            ctttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGa                                                                <  1:4552345/72‑1 (MQ=255)
            ctttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAAtt                                                            >  1:5725125/1‑76 (MQ=255)
             tttcttCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAAttt                                                           >  2:4790114/1‑76 (MQ=255)
               tcttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTg                                                         >  1:5305715/1‑76 (MQ=255)
               tcttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTg                                                         >  2:4537475/1‑76 (MQ=255)
                cttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTgg                                                        >  1:3185346/1‑76 (MQ=255)
                 ttcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGc                                                       >  1:1395406/1‑76 (MQ=255)
                  tcCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGcc                                                      <  1:126624/76‑1 (MQ=255)
                  tcCCAGGCGAGTCTGAGTATATGAAAGACGAGCATTTGGTATCATCATCCCTGAATTCAGAGATGAAATTGTGGcc                                                      >  1:3403268/1‑76 (MQ=255)
                   cccAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCa                                                     <  2:3157708/76‑1 (MQ=255)
                   cccAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCa                                                     <  2:3153532/76‑1 (MQ=255)
                    ccAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCAc                                                    >  2:5028695/1‑76 (MQ=255)
                    ccAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCAc                                                    <  2:4010210/76‑1 (MQ=255)
                     cAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACt                                                   <  2:5211171/76‑1 (MQ=255)
                      aGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGa                                                                <  2:4778293/62‑1 (MQ=255)
                      aGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGa                                                                >  1:4778294/1‑62 (MQ=255)
                      aGGCGAGTCGGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTc                                                  <  2:3159788/76‑1 (MQ=255)
                       ggCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCa                                                 >  2:1622891/1‑76 (MQ=255)
                        gCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGCGATTAAACTTTGGCCCCTCAc                                                <  1:1366838/76‑1 (MQ=255)
                        gCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                                >  1:1253071/1‑76 (MQ=255)
                        gCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                                <  2:4076339/76‑1 (MQ=255)
                        gCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCAc                                                <  1:594075/76‑1 (MQ=255)
                         cGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACg                                               <  1:5500092/76‑1 (MQ=255)
                          gAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGa                                              >  2:2675173/1‑76 (MQ=255)
                            gTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGt                                            >  1:2804214/1‑76 (MQ=255)
                            gTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGt                                            <  1:2695185/76‑1 (MQ=255)
                             tCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTg                                           <  2:162742/76‑1 (MQ=255)
                               tGAGAATATGAAAGCCGCGCATTTGCTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGCGGc                                         >  2:3899549/1‑76 (MQ=255)
                                gAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGcc                                        >  2:5238029/1‑76 (MQ=255)
                                 aGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTg                                                         <  1:1417069/58‑1 (MQ=255)
                                 aGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTg                                                         >  2:1417069/1‑58 (MQ=255)
                                   tataTGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttt                                     >  2:4937059/1‑76 (MQ=255)
                                   tataTGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttt                                     >  1:5143366/1‑76 (MQ=255)
                                   tataTGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttt                                     >  2:4989591/1‑76 (MQ=255)
                                   tataTGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttt                                     >  1:4384157/1‑76 (MQ=255)
                                    atatGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCtttt                                    >  2:2868969/1‑76 (MQ=255)
                                     tatGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCttttt                                   <  1:1622891/76‑1 (MQ=255)
                                        gAAATACGCGCATTTGTTATCATCATCCCTGAATTCAGAGAGGAAAGTTTGGACACGCACGAGTGGCCTttttctt                                <  1:2669151/76‑1 (MQ=255)
                                        gAAAGACGCGCATTTGTTATCATCATCCCTGAATTCGGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttctt                                <  1:2631046/76‑1 (MQ=255)
                                         aaaGACGCGCATTTGTTATCATCATCCCTGAATTCAGGTGTGAAATTTTGGCCACTCGGGCGTGGCCTttttcttt                               <  2:2430657/76‑1 (MQ=255)
                                         aaaGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttt                               <  2:1555088/76‑1 (MQ=255)
                                         aaaGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttt                               >  1:1555088/1‑76 (MQ=255)
                                         aaaGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttt                               >  1:4072096/1‑76 (MQ=255)
                                           aGACGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTttttcttttc                             <  2:2921853/76‑1 (MQ=255)
                                             aCGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTg                           <  2:5109902/76‑1 (MQ=255)
                                             aCGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTg                           >  2:4890833/1‑76 (MQ=255)
                                             aCGCGCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTCTTCTTTTCTg                           >  1:2452892/1‑76 (MQ=255)
                                              cgcgCATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGt                          >  1:2891078/1‑76 (MQ=255)
                                                cgcATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCa                        >  2:179462/1‑76 (MQ=255)
                                                cgcATTTGTTATCATCAGCCCTGAATTCAGAGATGAAATTTTGGCAAATCAAGAGTCACATTGTTCTTTTCTGTCa                        >  1:3711621/1‑76 (MQ=255)
                                                 gcATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAg                       <  1:4743959/76‑1 (MQ=255)
                                                  cATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGt                          <  1:1332547/72‑1 (MQ=255)
                                                  cATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGt                          >  2:1332547/1‑72 (MQ=255)
                                                  cATTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAgg                      >  2:4877751/1‑76 (MQ=255)
                                                   aTTTGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGc                     >  1:1410479/1‑76 (MQ=255)
                                                    tttGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCg                    >  2:221954/1‑76 (MQ=255)
                                                    tttGTTATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCg                    >  2:3765461/1‑76 (MQ=255)
                                                        ttATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGtt                >  2:2924316/1‑76 (MQ=255)
                                                         tATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGttt               <  2:4072095/76‑1 (MQ=255)
                                                         tATCATCATCCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGttt               >  2:4791725/1‑76 (MQ=255)
                                                           tcatcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGttttt             <  1:2973815/76‑1 (MQ=255)
                                                            catcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTc            <  1:3484629/76‑1 (MQ=255)
                                                            catcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTc            <  2:3415810/76‑1 (MQ=255)
                                                            catcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTc            <  1:1968307/76‑1 (MQ=255)
                                                            catcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTc            >  1:935670/1‑76 (MQ=255)
                                                             ctcatcCCTGAAAGCAGAGATGACATTTTGGCCGCTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTcc           <  2:816391/75‑1 (MQ=255)
                                                              tcatcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGTGTGTATTTCCa          >  2:4229122/1‑76 (MQ=255)
                                                               catcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAg         <  1:5438653/76‑1 (MQ=255)
                                                                atcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGc        >  2:301617/1‑76 (MQ=255)
                                                                atcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGc        >  2:3899393/1‑76 (MQ=255)
                                                                 tcCCTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCCGcc       >  2:5757247/1‑76 (MQ=255)
                                                                  cccTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCATGCGTGTTTTTCCAGCca      >  2:4421122/1‑76 (MQ=255)
                                                                   ccTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCcac     <  2:2891078/76‑1 (MQ=255)
                                                                   ccTGAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCcac     <  2:5725123/76‑1 (MQ=255)
                                                                     tGAATTCAGAGTTGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCacac   >  1:3283932/1‑76 (MQ=255)
                                                                      gAATTCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACAcc  <  1:179462/76‑1 (MQ=255)
                                                                           |                                                                       
ATTTTCTGAACTCTTTCTTCCCAGGCGAGTCTGAGTATATGAAAGACGCGCATTTGTTATCATCATCCCTGAATTCCAGAGATGAAATTTTGGCCACTCACGAGTGGCCTTTTTCTTTTCTGTCAGGCGTGTTTTTCCAGCCACACC  >  NC_000913/2812783‑2812929

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